| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,423,290 | 0 | A | C | 42.9% | 31.3 / 10.9 | 35 | intergenic (+96/+133) | yhdZ/rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
| Reads supporting (aligned to +/- strand): ref base A (5/15); new base C (14/1); total (19/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.87e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGG > NC_000913/3423145‑3423430 | aTCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTc < 1:492664/151‑1 (MQ=255) tAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCTcccg > 1:82168/1‑148 (MQ=255) cAGCGCGTACTAGGGCATTTTGACCGTAGTTAATCCATTAATTTAATGTTAGTTCGTAAATCAAAAAGGCCATCTTTTCGGCTGCCTTTACGCTTGATTTGATGCCTGGCAGTTAATGGCGGGCGTCCGG‑CCCGCCACCCTCCGGGCCGtt < 2:131557/151‑1 (MQ=18) agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCCTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGAAGTTTATGTCGGGCGTCCTG‑CCCGCCCCCCCCCGGGCCGTTg < 2:122965/151‑1 (MQ=255) cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTGTGGCGGGGGTCCCG‑CCCCCCACCCCCCCGGCCGGTG‑Ct > 1:560783/1‑149 (MQ=255) aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTT‑TCCCCCCCCCCCCCGGCCGTTT‑TTTCGc > 1:31926/1‑150 (MQ=37) aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGGCTGGCAGTTTATGGCGGGCGGCCCG‑CCCGCCACCCCCCggcgcccttgt > 2:212427/1‑137 (MQ=255) gAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGGccctggtttgcaa > 1:210934/1‑139 (MQ=255) gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACg < 2:129219/150‑1 (MQ=255) tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTT‑CCCCCCCCCCCCCGGGCCGTTG‑CTTTGCAACGTTCaa > 1:471293/1‑150 (MQ=255) tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTGGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATcc < 2:268410/150‑1 (MQ=255) cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTT‑CCCCCCACCCCCCCCGCCCTTGCCTTCCCAAcctcaaattctccccc > 2:697971/1‑136 (MQ=255) aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTG‑CCCCCCCCCCCCCCGGCCCCTG‑CTTCTCCACCTTCAAATCCCCTCCccgcg > 2:281631/1‑150 (MQ=255) aTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGa < 2:667005/151‑1 (MQ=255) ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTT‑CCCCCCCCCCCCCCCGCCCCTG‑CTTTTCCACGTTCAAATCCCCCCCCCGCGGGtt > 2:348836/1‑151 (MQ=18) ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGAtt < 2:129789/151‑1 (MQ=255) aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCGCCG‑CCCCCCCCCCCCCCGGCCGCTT‑CTTCTCAACGTTCAAATCCCCTCCCCGCGGGttt > 2:352222/1‑150 (MQ=255) ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCTCCCGGCCGTTG‑CTTCGCAACGTTCAAAACCCCTCCCGGCGGGtt > 1:729059/1‑148 (MQ=18) ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCCCCCGGGCGGTG‑CTTTGCAACGTTTAAATCCGCTCCCCGCGGGttt > 1:508541/1‑149 (MQ=18) ttAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTG‑CCCGCCCCCCCCCCGGCCGGTG‑CTTCGCAACGTTCAAATCCGCCCCCCGCGGGttt > 1:785948/1‑149 (MQ=18) gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCggg < 1:372436/151‑1 (MQ=255) ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCG‑CCCCCCCCCCCCCCGGCCGTTG‑CTTTGCAACGTTCAAATCCCCTCCCCGCGGGTTTTTCCTCCTCcgg > 1:125005/1‑150 (MQ=17) cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCCCCCGGGCGTTG‑CTTTGCCACCCTCAAATCCCCTCCCCGCGGGGTTTTCCTTCTCtgggggggg > 2:153473/1‑144 (MQ=17) gAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCG‑CCCCCCCCCCTCCCCGCCGTTG‑CTTTGCAACGTTTAAATCCCCTCCCGGCGGGTTTTTTCTTCCTGGgaga > 1:417897/1‑146 (MQ=255) aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAc < 1:727471/151‑1 (MQ=255) aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAcc < 1:371476/151‑1 (MQ=255) ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGAt < 2:341738/151‑1 (MQ=255) ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATaa < 1:646508/150‑1 (MQ=255) ggCTTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAg < 1:174418/151‑1 (MQ=255) ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTCCTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGc < 2:130896/151‑1 (MQ=255) cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg < 1:51996/150‑1 (MQ=255) atttgatGTCTGGCAGTTTATGGCGGGCGTCCCG‑CCCCCCCCCCTCCCGGCCGTTG‑CTTTGCAACGGTCAAATTCCCTCCCCGCGGGTTTTTCCTACTCCGGGGAGGGTTCACCCACCAACAAAAGATAAAAAAAAAAGACCAGTCtttc > 1:45316/1‑148 (MQ=2) tGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa < 2:119677/150‑1 (MQ=38) gTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt > 1:429661/1‑151 (MQ=40) ccTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTTACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg > 2:515555/1‑151 (MQ=255) | ATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTG‑CTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGG > NC_000913/3423145‑3423430 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |