| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,604,937 | 0 | T | G | 51.4% | 39.2 / 59.1 | 35 | V182G (GTG→GGG) | rsmD | 16S rRNA m(2)G966 methyltransferase, SAM‑dependent |
| Reads supporting (aligned to +/- strand): ref base T (0/17); new base G (18/0); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.20e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-08 | |||||||||||
TGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑TTTGGGGATTTTTAATCCTCAACCTGGTGCATCCCTTCCCACGCCCG > NC_000913/3604793‑3605062 | tGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCt < 2:190082/149‑1 (MQ=255) gaagaGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGggggggtattgg > 1:263541/1‑142 (MQ=255) agaCGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGCGGGGTCAGGGGGtttt > 1:745444/1‑142 (MQ=255) aaGATAACGGCTGGCTGGCTGACGACGCCCTGATTTATGTCGAAAGCGAAGTCGAAACCGGTCTGCCCACTGTCCCAGCAAACTGGTCATTCCATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTTTCGACGCGAAG‑TACAAgg < 1:664486/151‑1 (MQ=255) gATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGgag < 2:526740/151‑1 (MQ=255) gATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGgag < 1:617496/151‑1 (MQ=255) tAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTg < 2:616435/127‑1 (MQ=255) tggctggctgACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTTATCGGCGGTTTCAAAACGAAA‑GACAAAGAGaaa > 2:625453/1‑145 (MQ=255) gctgACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGAtt < 2:23226/151‑1 (MQ=255) gATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCGGGGGGGTTTTTTGCTGTATTAACCCGAAG‑CACAAGGGGAAAGAGATGGTGGTTTATA‑TTTGTCGttt > 2:619562/1‑151 (MQ=255) aTTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTg < 2:245039/151‑1 (MQ=255) aTGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTa < 1:2342/150‑1 (MQ=255) gTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAAt < 2:567529/150‑1 (MQ=255) gTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATg < 1:32913/151‑1 (MQ=255) cGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTGAGGGGGCTTTTTGGGTTTTTTAACACCAAG‑GACAAGGAGAAAGAGGTGTTGATTAATA‑TTTGTGGTTTGTTTATGct > 2:734639/1‑151 (MQ=255) cGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTTATTGGCTGGATTAACGCGAAG‑CACAAAGAGGAAAGGGTGGTGGTTTATA‑TTTGTTGTTTGTTTATGct > 2:271884/1‑151 (MQ=255) cGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑t < 1:252861/151‑1 (MQ=255) cGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑t < 1:38707/151‑1 (MQ=255) aaGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑ttt < 1:614413/151‑1 (MQ=255) aGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTCGTTGGGTGTATTAAACCGGAG‑GACAAGGAGAAAAGGGTGCTGGTTTAAA‑TTTGTTGTTTGTTTATggttttcggtt > 1:776216/1‑140 (MQ=255) aaaCGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTTTTTGGCTGTATCAAAGCGAAG‑CACAAAGAGAAAGGGAGGGTGGTTAATA‑TTTGTCGGTTTTTTATGGTCTGC‑G‑TTTtgggggt > 2:192481/1‑148 (MQ=255) aaaCGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGGTATCGGCTGGTTTCACACGAAG‑GACAAAGAGAAAAGGAAGGTGTTTTATA‑TATTTCGTTTGTTTATATTCTGT‑GCTTTgggg > 1:705145/1‑149 (MQ=255) ggTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑TTTGGGGATTTTTa < 1:560901/151‑1 (MQ=255) gCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCCGGGGGCTTTTTGGCGGGTTTAAAGCGAAG‑GAAAAGGGGGAAAGGATGCTCGTTTTTT‑TTGGGTGTTTGTTTTTTGTCTTC‑G‑TTTGGGGGtttttttt > 1:443641/1‑145 (MQ=255) cccACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCCAGGGGGGTATCTGGCGTGTCAACGCGAAG‑GAAAAAGAGGAAAGGATGGTGATTTATT‑TTTTTTGGTTTTTAATTCTCTTC‑G‑TTTTGGGTTTTTTTATTCt > 2:571590/1‑150 (MQ=255) tGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTTTTTCGGTGTATTAACACGGAA‑CACCAAGGGAAAAGGGGGGTGGTTTATA‑TTTTTTGGTTTTTTATTCTCTTC‑GTTTTGGGGGTTTTTat > 2:609642/1‑141 (MQ=255) tCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCATGGGGCGTATTTGGTGTTTTCACGCGAAGAAAAAAGGGGAAAGGGA‑GGTGGTTTTTA‑TTTGTGGTTTGTTTATTATTTGC‑G‑GTTGGGGGgttttttt > 1:504983/1‑136 (MQ=255) aGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑TTTGGGGATTTTTAATCCTCAACCTGGTGCa < 1:117286/151‑1 (MQ=255) aGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGGTAGGGGGGCTATTCGCTGGGTTAAAGCGGAG‑AAAAAAGGGAAAAAGATGGTTTTTTATA‑TTTGTTGGTTGTTTATTGTTTGCGG‑TTTGGGGTTTTTTATTTTTCACCCTg > 2:19529/1‑147 (MQ=255) cAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑TTTGGGGATTTTTAATCCTCAACCTGGTGCATc < 2:328175/151‑1 (MQ=255) aaaCTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑TTTGGGGATTTTTAATCCTCAACCTGGTGCATcc < 2:673678/151‑1 (MQ=255) aaaCTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTTTTTGGGTGTATCAACGCGAAG‑CACAAGGGGAAAAAGATGGTGGTTTATA‑TTGGGCGTTTTTTTATTCTCTGCGG‑TTTGGGGGTTTTTATTCTCCACCCTgtggatcc > 2:154939/1‑144 (MQ=255) aCTGGTCATTACATCGGGAAAAAGTGGC‑GGGGCTGGGGGGGTATCGGGCGTTTCAAAGCGAAA‑AACAAAGGGAAAGGGAAGGTGGTTTAAA‑TTTGGCGTTTGTTTATTTTTTTC‑G‑TTTTGGGGGTTTTTATTTTCAACCTGGTGGCTccc > 1:291890/1‑150 (MQ=255) aTTACATCGGGAAAAAGTGGC‑GGGTCAGGGGGGCTATCGGGTGGATTAAAGCGGAG‑GACCAAGGGAAAATGAGGGTGGTTAAAATTTGGTGGTTGTTTAAT‑CTTTGC‑G‑GTTGGGGGTTTTTTATTCTCAACCCTGGGGATCCCt > 2:313928/1‑144 (MQ=255) cATCGGGAAAAAGTGGC‑GGGTCAGGGGGGGTTTTGGCTGTATCAACGCGGAA‑CAAAAGGAGGAAGAGGTGGTGATTTATA‑TTTGTGGTTTGTTTATTCTCTGCGG‑TTTGGGGGTTTTT‑ATTCTCAAACTGGTGCATTCCTTCCCACccccc > 1:510814/1‑150 (MQ=255) cATCGGGAAAAAGTGGC‑GGGGCAGGGGGGTTTTTGGCTGTGTTAAAGCGAAG‑CAAAAAGGGAAAAAGGTGGTGGTTTATA‑TTTGTTGTTTTTTTATGCTCTGC‑G‑TTTGGGGGGTTTTTATTCTCAACCTGGTGCATCCCTTCcccccccc > 2:272543/1‑150 (MQ=255) | TGTTAGAAGAGACGATAAATTTACTGGAAGATAACGGCTGGCTGGCTGACGAAGCCCTGATTTATGTCGAAAGCGAAGTCGAAAACGGTCTGCCCACTGTTCCAGCAAACTGGTCATTACATCGGGAAAAAGTGGC‑GGGTCAGGTGGCTTATCGGCTGTATCAACGCGAAG‑CACAAGGAGAAAGTGATGCTGATTAATA‑TTGGTCGTTTGTTAATGCTCTGC‑G‑TTTGGGGATTTTTAATCCTCAACCTGGTGCATCCCTTCCCACGCCCG > NC_000913/3604793‑3605062 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |