| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 4,171,438 | T→A | 68.6% | noncoding (2798/2904 nt) | rrlB → | 23S ribosomal RNA of rrnB operon |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,171,438 | 0 | T | A | 68.6% | 85.4 / 46.6 | 57 | noncoding (2798/2904 nt) | rrlB | 23S ribosomal RNA of rrnB operon |
| Reads supporting (aligned to +/- strand): ref base T (7/11); new base A (12/27); total (19/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.60e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01 | |||||||||||
GTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAG > NC_000913/4171295‑4171586 | gTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGG > 3:779432‑M1/1‑149 (MQ=255)gTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGG < 4:1314497‑M1/149‑1 (MQ=255)gTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGG < 3:1074696‑M1/149‑1 (MQ=255)gTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGG < 2:993044‑M1/149‑1 (MQ=255)gTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGG > 2:329754‑M1/1‑149 (MQ=255)gTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGG > 2:1539497‑M1/1‑149 (MQ=255) tACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGT < 3:684646‑M1/149‑1 (MQ=255) tACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGT > 4:1425856‑M1/1‑149 (MQ=255) tACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGT < 1:357323‑M1/149‑1 (MQ=255) aCGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTC < 3:1425856‑M1/149‑1 (MQ=255) aCGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTC < 1:126801‑M1/149‑1 (MQ=255) aCGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTC < 1:1455357‑M1/149‑1 (MQ=255) aCGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTC < 2:1259692‑M1/149‑1 (MQ=255) gagaGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCag < 3:91938‑M1/149‑3 (MQ=255) gagaGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCag < 3:577179‑M1/149‑3 (MQ=255) gagaGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCag < 1:126383‑M1/149‑3 (MQ=255) agagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 1:1249379‑M1/149‑4 (MQ=255) agagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 3:234842‑M1/149‑4 (MQ=255) agagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 4:1278918‑M1/149‑4 (MQ=255) agagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 2:495230‑M1/149‑4 (MQ=255) agagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 4:831978‑M1/149‑4 (MQ=255) agagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 1:1386921‑M1/149‑4 (MQ=255) gagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg > 1:1246119‑M1/1‑145 (MQ=255) gagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg < 4:1115901‑M1/149‑5 (MQ=255) gagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg < 2:802540‑M1/149‑5 (MQ=255) gagGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg < 4:143176‑M1/149‑5 (MQ=255) agGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggga < 4:866527‑M1/149‑6 (MQ=255) agGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggga < 1:1319134‑M1/149‑6 (MQ=255) ggACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagggag < 2:1246119‑M1/149‑7 (MQ=255) gACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagggaga < 1:586662‑M1/149‑8 (MQ=255) ccGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 2:738378‑M1/143‑4 (MQ=255) ccGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg > 1:738378‑M1/1‑140 (MQ=255) aTGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg < 4:783678‑M1/110‑5 (MQ=255) aTGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg > 3:783678‑M1/1‑106 (MQ=255) aTGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg > 3:1473777‑M1/1‑106 (MQ=255) aTGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 4:1473777‑M1/109‑4 (MQ=255) ccAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 3:440014‑M1/106‑4 (MQ=255) ccAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg > 4:440014‑M1/1‑103 (MQ=255) aaTGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTC < 2:797860‑M1/101‑1 (MQ=255) aaTGGCACTGCCCGGTAGCTAAAAGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTC > 1:797860‑M1/1‑101 (MQ=255) cTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg > 4:1463682‑M1/1‑94 (MQ=255) cTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg < 1:1139113‑M1/98‑5 (MQ=255) cTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg < 3:1463682‑M1/98‑5 (MQ=255) cTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCaggg > 2:1139113‑M1/1‑94 (MQ=255) tGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg > 1:1311772‑M1/1‑93 (MQ=255) tGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 2:1311772‑M1/96‑4 (MQ=255) ggAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCag < 3:204921‑M1/77‑3 (MQ=255) ggAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCag > 4:204921‑M1/1‑75 (MQ=255) gAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagggagaa > 3:1570734‑M1/1‑74 (MQ=255) gAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagggagaa < 4:1570734‑M1/82‑9 (MQ=255) tAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagggagaa < 2:1255697‑M1/75‑9 (MQ=255) tAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAATTCTCCCTGACCCTTAAAGGGTCagggagaa > 1:1255697‑M1/1‑67 (MQ=255) aTCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg > 3:1209663‑M1/1‑53 (MQ=255) aTCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTAAAGGGTCagg < 4:1209663‑M1/56‑4 (MQ=255) ttCTCCCTGACTCCTTGAGAGTCagggagaactcatctcggggcaagtttcgtgcttagatgctttcagcacttatctcttccgcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggt > 1:1098263‑M1/1‑23 (MQ=255) tctcCCTGACTCCTTGAGAGTCagggagaactcatctcggggcaagtttcgtgcttagatgctttcagcacttatctcttccgcatttagctaccgggcagtgccattggcat > 3:809361‑M1/1‑22 (MQ=255) tctcCCTGACTCCTTGAGAGTCagggagaactcatctcggggcaagtttcgtgcttagatgctttcagcacttatctcttccgcatttagctaccgggcagtgccattggcat < 4:809361‑M1/113‑92 (MQ=255) tttAAGGGTCagggagaactcatctcggggcaagtttcgtgcttagatgctttcagcacttatctcttccgcatttagctaccgggcagtgccattggcatgacaacccgaaccccagtgatgcgtccactccggtcctctcgtactag < 4:577179‑M1/149‑140 (MQ=255) tAAGGGTCagggagaactcatctcggggcaagtttcgtgcttagatgctttcagcacttatctcttccgcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactagga < 3:1278918‑M1/149‑142 (MQ=255) | GTACGAGAGGACCGGAGTGGACGCATCACTGGTGTTCGGGTTGTCATGCCAATGGCACTGCCCGGTAGCTAAATGCGGAAGAGATAAGTGCTGAAAGCATCTAAGCACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAG > NC_000913/4171295‑4171586 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |