| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 1,921,182 | (C)5→6 | coding (960/1440 nt) | rsmF → | 16S rRNA m(5)C1407 methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,921,177 | 1 | . | C | 93.1% | 94.5 / 0.2 | 29 | coding (955/1440 nt) | rsmF | 16S rRNA m(5)C1407 methyltransferase |
| Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (12/15); total (14/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.24e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCG‑C‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAA > NC_000913/1921041‑1921311 | aCTTGGCGATCTCTTCCCTGGTTCAAACAAAACGCCGACTGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCG‑c‑cc < 2:327705/140‑1 (MQ=255)aCTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACGGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTAACGACTGCAAAGCCTTCTTCGTTGCTCGTCTGCGTAAACCTCAGGCGATCCCCGCCTTACCCG‑c‑cc < 2:971466/140‑1 (MQ=255) ttGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGGCCGAAGAAGGCTTTTTGCATGTTTCCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCG‑c‑cc‑c‑c < 2:882789/140‑1 (MQ=255) ttGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCG‑c‑cc‑c‑c > 2:89962/1‑140 (MQ=255) ggCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCG‑c‑cc‑c‑cca > 2:922529/1‑138 (MQ=255) tGGTGCAAACAAAGGGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTGAAACTCAGACGATTCCAGCCTTACCCG‑C‑CCGC‑CAAATACAAAGTCGGt > 2:381124/1‑140 (MQ=255) tGCAAACAAAGCGCTGAGCGTAGAAGTCTTTTTGCATGTTTTCCCACACATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAAt > 2:745357/1‑140 (MQ=255) tGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGGTTCACGCGTTACCCG‑CGCC‑C‑CTAATACAAAGTCGGTAAt > 2:293230/1‑140 (MQ=255) caaacaaaGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAAttt > 1:72005/1‑140 (MQ=255) aaGGGCTGACCGAAGTAGGCTTCTTTCATTTTTTCCCGCAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGtt < 1:381124/140‑1 (MQ=255) aGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTc > 1:604446/1‑140 (MQ=255) cTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGcc < 1:816478/140‑1 (MQ=255) cTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGcc < 2:107838/140‑1 (MQ=255) cTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGcc < 2:92446/140‑1 (MQ=255) cTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACGGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGcc < 2:689942/140‑1 (MQ=255) cGAAGAAGGCTTTTTGCATGTTTTCCCACAATTTTCCGATTGCGAAGGGTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTg < 2:393794/140‑1 (MQ=255) tttGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGaa > 2:729828/1‑140 (MQ=255) cccACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAAtt > 1:911480/1‑140 (MQ=255) acaAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGt < 1:729828/140‑1 (MQ=255) aaTTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAg < 2:808036/140‑1 (MQ=255) tttACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAggc < 2:911480/140‑1 (MQ=255) gACTGCGAAGGCTTCTTCGTTGCTCCTCTTCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTa < 1:353304/140‑1 (MQ=255) aCTGCGAAGGCTTCTACGTTGCCCGACTGCGTAAAACACAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTAc < 1:300040/140‑1 (MQ=255) cGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGGCCGGTGAAGGATTGCGGAGCTGGACGAATTCGTCAGGCGGCTACAggg > 2:300040/1‑139 (MQ=255) aaGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgtgt < 2:604446/140‑1 (MQ=255) aGGCTTCTTCGTTGCTTCTCTGAGTAAAACTCAAGCGATCCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGtt < 2:445892/140‑1 (MQ=255) ttcttcGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCt < 2:72866/140‑1 (MQ=255) tCTGCGTAAAACTCGGGCGATTCCCGCCTTATCCG‑C‑CC‑CGCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGaa < 2:433445/140‑1 (MQ=255) tCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGTCTAGAGGTGTTGGATTAAACTGGGATGaa > 1:433445/1‑140 (MQ=255) tCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATCCGTCAGGCGGCTACAGCTGTTGGCTTAAACTGGGATGaa > 2:403008/1‑140 (MQ=255) cTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGaaa > 1:34136/1‑140 (MQ=255) cGTAAAACTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAAcc > 1:744150/1‑140 (MQ=255) atAAAGCTCAGGCGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCt < 1:745357/139‑1 (MQ=255) cGATTCCCGCCTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCa < 2:744150/140‑1 (MQ=255) ccTTACCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACaa > 1:21187/1‑140 (MQ=255) aCCCGCC‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGaa > 1:166770/1‑140 (MQ=255) | ACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCG‑C‑CC‑C‑CAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATGAAAACCTGCGCCTCTGGCAGCGTGACAAAGAA > NC_000913/1921041‑1921311 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |