| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,423,287 | 0 | G | C | 51.3% | 27.7 / 53.5 | 37 | intergenic (+93/+136) | yhdZ/rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
| Reads supporting (aligned to +/- strand): ref base G (3/15); new base C (19/0); total (22/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.71e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-05 | |||||||||||
TCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATG > NC_000913/3423146‑3423425 | tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCTT‑cccccccccccc > 2:92004/1‑148 (MQ=255)tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCCG‑cccccccccccc > 1:63393/1‑148 (MQ=255)tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑cccccccccccc > 2:557162/1‑148 (MQ=255)tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTG‑cccccccccccc > 2:215684/1‑148 (MQ=255)tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTG‑CCCGcccccccc > 2:403790/1‑148 (MQ=255)tCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTC‑cccccccccccc > 2:210034/1‑148 (MQ=255) cTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCT‑Ccccccccccccccg > 2:572896/1‑149 (MQ=255) aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGc < 1:134030/151‑1 (MQ=255) aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTT‑CCCCCCCCCCCCCggg > 2:496781/1‑150 (MQ=255) gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCCGGCCGttt > 1:464048/1‑150 (MQ=37) gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCG‑CCCCCCCCCCCCCCGGCCGGTGttt > 2:758514/1‑151 (MQ=255) gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTT‑CCCCCCCCCCCCCCGGGCGTTGCTTc > 2:180741/1‑150 (MQ=37) tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCG‑Cccccccccccccccggcccttgtttg > 1:779514/1‑136 (MQ=255) tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCG‑CCCCCCCCCCCCCCGGCCGTTGCtttg > 2:477231/1‑148 (MQ=255) aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCG‑CCCCCCACCCCCCCCGCCGCTGCTTTTC‑AAc > 1:756537/1‑150 (MQ=255) gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCG‑CCCCCCCCCCCCCCCGCCGTTgtctctcaca > 1:86280/1‑140 (MQ=255) ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTT‑CCCCCCCCCCTCCCGGCCCGTGCTTCTC‑AACGTTCAAAt > 1:587777/1‑151 (MQ=25) tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATc < 2:393774/151‑1 (MQ=255) tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATc < 2:792122/151‑1 (MQ=255) ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATcc < 2:602174/151‑1 (MQ=255) gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCG‑CCCCCCCCCCCCCCGGCCGGTGCTTTGC‑AACCTTTCAATcccccccccg > 2:408849/1‑148 (MQ=21) ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGAtt < 1:27425/151‑1 (MQ=255) ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCCCCCGGCCGCTGCTTTGC‑CAAGTTTAAAATCGCCCCccgcggggt > 2:772102/1‑148 (MQ=18) aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTcc < 2:647097/151‑1 (MQ=255) aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATTGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCACCCCCCCGGGCGGTTCTTCGC‑AAAGATTAAAACCCCCCCCGGCGGGGTTGTc > 2:538078/1‑149 (MQ=18) gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTT‑CCCGCCACCCCCCCGGGCGTTTCTTTGC‑AAACTTATAATCCCCTCCCCGCGGGGTTGTCCtt > 2:651382/1‑148 (MQ=255) aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCC‑CCCCCCACCCCCCCGGGCGCTGTTTTCT‑AAAAATCAAAATCGCTCCccgcgggattttgcttctc > 2:612013/1‑137 (MQ=12) ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGGCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCCCCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCggg < 1:424333/150‑1 (MQ=255) tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGgag < 1:496543/151‑1 (MQ=255) cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCCCCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCAcc < 2:634788/149‑1 (MQ=255) cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCAcc < 2:392186/149‑1 (MQ=255) ccgtcaggATGGCCTTCTGCTTAATTTGATGCCTGGCAGTTTATGGCGGGCTGCCTT‑ACCCCCCCCCCCCCTGCCCTTGCTTTTC‑AATCGTCAAATTCTCTTCCCGCGGTCTTTTTTTA‑CCCCCGAGAGTGTTCTCCGCCCTacaacac > 1:199758/7‑148 (MQ=0) ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTT‑CCCCCCCCCCCCCCGGCCGGTGCTTCCC‑CACCTTCAAAACCCCCCCCCGCGGGTTTTTCCTA‑CTCCGGAGGGGGTTTACCCacaaaaaaaaaaa > 2:552307/1‑146 (MQ=9) gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCC‑CCCCCCCCCCCCCCGGCCGGTGCTTTGC‑AAAGGTCAAATCCCCCCCCCGCGGGGTTTTTCTT‑CTCCGGGGGGGGGTCACCCCCaaaaaaaaaaaaaaaaaaaaag > 2:352581/1‑149 (MQ=9) ccTTTCGCTTGATTTGATGTCGGGCAGTTTATGGCGGGCGCCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGc < 1:403578/151‑1 (MQ=255) gatttgatGTCTGGCAGTTTATGGCGGGCGCCCTG‑CCCCCCACCCCCCGGGCCGCTGCTTTGC‑AACGTTCAAAACCCCCCCCCGCGGGTTTTTCCTT‑CTCCGGAGAGGGTTTACCCACAAAAAAAAGAAAAAAAAAAAAGCCCCGGCTTc > 2:87560/1‑150 (MQ=255) atttgatGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCCCCCCCCCGGCCGGTGCTTCGC‑CAACTTTAAATCCCCTACCCGCGGGGTTGTTCTC‑CTCCGGCGGGTGTTCACCCACAAAACAAAGATAAAAAAAAAAGCCCCGttcttcg > 2:148155/1‑145 (MQ=2) atttgatGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCCCCCGGCCGGTGCTTCTCAAAGGTTCAAACCGCCCCCCGGCGG‑TTTGTCCTACCCCGGGGGGGGGTTCCCGAACAAACAACAAAAAAAACAAAAGACCCCGCCTTc > 1:317514/1‑150 (MQ=2) atttgatGTCTGGCAGTTTATGGCGGGCGGCCCG‑CCCCCCCCCCCCCCGGCCCGTGCTTTGC‑AAAGGTTAAATTCCCTCCCCGCGGGTTTTTCCTT‑CTCCGGGGAGGGTTTACCCCCAAACAACAGAAAAAAAAAAAAAGCCAGGCTTcc > 2:523413/1‑150 (MQ=255) ttgatGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc < 1:792889/151‑1 (MQ=38) gatGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCCCCCCCCCCCCGGCCGGTGCTTCGC‑CACCCTCCAATCCGCCCCCCGCCGGTTTTTCCTA‑CTCCGGGGGGTGGTCACCCACCAACAACAAATAAAAAAAAAAGCCCACTCCTCCggccg > 1:237480/1‑147 (MQ=2) tttATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt < 2:393855/149‑1 (MQ=38) gcgggcgTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt < 1:500565/150‑1 (MQ=38) gcgTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg < 2:566356/150‑1 (MQ=38) | TCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTG‑CCCGCCACCCTCCGGGCCGTTGCTTCGC‑AACGTTCAAATCCGCTCCCGGCGGATTTGTCCTA‑CTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATG > NC_000913/3423146‑3423425 |
| 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: ‑ |