| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,423,287 | 0 | G | C | 41.2% | 47.8 / 17.3 | 34 | 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 (6/14); new base C (14/0); total (20/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.05e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01 | |||||||||||
AATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTC > NC_000913/3423151‑3423436 | aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCT‑G‑Ccccccccccc‑cc > 1:609702/1‑146 (MQ=255)aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCT‑G‑CCCCCCACCCC‑CCcgggc > 1:91893/1‑149 (MQ=255) cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCTCCCT‑CCGGGCCGt < 1:607627/150‑1 (MQ=255) gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCt < 1:768487/150‑1 (MQ=255) gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGGCC‑T‑CCCCCCCCCCT‑CCCGGCCGGTTCtt > 1:715490/1‑149 (MQ=255) gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCg < 1:747367/151‑1 (MQ=255) aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑C‑CCCCCCCCCCC‑CCGGGCCGTTGCTTTGCAACCTTCAAATcccc > 1:362452/1‑149 (MQ=21) aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCGGCCCT‑G‑CCCCCCACCCC‑CCCCGCCGCTGCTTTCCAAccatctaatac > 1:606495/1‑139 (MQ=255) gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGGCCC‑G‑CCCGCCCCCCC‑CCCGGCCCGTTCCTCGCAAAGTTCAAATCCCCTcc > 1:590890/1‑151 (MQ=255) tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCT‑G‑CCCCCCCCCCC‑CCCGGCCGGTTCTTTGCCACGTTTAAATCCCCCCCccgc‑ggg > 2:173246/1‑149 (MQ=21) aaTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTc < 2:546154/150‑1 (MQ=255) ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCCCCCCCCC‑CCCGGCCGCTGCTTTGCAAAGTTTAAAACCCCCCCCCGC‑GGGTTTGt > 2:419759/1‑147 (MQ=18) gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACt < 1:192537/151‑1 (MQ=255) aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCC‑C‑CCCCCCCCCCC‑CCCGGCCGGTTCTTTTCAACCTTTACAACCCCCCCCCGC‑CGGTTTTTCCTAc > 1:726710/1‑149 (MQ=16) aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATAGCGGGAGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAAAGTTCAAAACCGCTCCCGGC‑GGAGTTTTCCTACTc > 1:127586/1‑151 (MQ=18) aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCT‑C‑TCCCCCCCCCC‑CCCGGGCGCTGCCTTGCAACGTTCCAAACCCCCCCCCGC‑GGGTTTTTccctctcc > 2:102113/1‑145 (MQ=255) tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCCCCCC‑CCCGGCCGGTGCTTTGCAACGTTTAAAACCCCCCCCCGC‑GGGTTTTTTCTACTCCGGGGAgt‑g > 2:601617/1‑151 (MQ=17) tCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGAGCCCC‑T‑CCCCCCCCCCTCCCGGCCCGTTG‑TTCCCAAACTTCAAATCCACCCCCCGC‑CGGTTTTTTctcctcccgagg > 1:255419/1‑142 (MQ=255) aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑Gtt < 2:235502/150‑1 (MQ=255) aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCCCCACCCC‑CCGGGCCGGTGCTTCTCCACGTTTCAATTCCCCCCCCGC‑GGGGTTGTTCTACTcctgcgggtgtt > 1:100980/1‑138 (MQ=17) cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGa < 2:3276/151‑1 (MQ=255) ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGAt < 1:125211/151‑1 (MQ=255) tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGG‑CCCCCCCCCCC‑CCGGGCCGTT‑TTTTGCAACGGTTCAATACCCCCCCCGC‑GGGTTTTTTCTACTCGGGGGGGT‑GGTGACCCACAAAAAACAGaaaaaaaaaa > 2:14616/1‑150 (MQ=255) gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCCTCCC‑GCCCCCCCCCCCC‑CCGGCCCGTTCTTTTCTAACTTTCAATTCCCCCCCCGCC‑GGGTTT‑TTCTCCTCCGGGGAGTGGTTCTCC‑CCAAACAACAGCGAaaaaaaaaa > 1:98856/1‑149 (MQ=255) ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑T‑CCCCCCCCCCT‑CCCGGCCGTTTCTTCCCCACGTTTAAATACACCCCCCGC‑GGATTTGTTCTACCCTGGAGCAG‑GTTCACCCACAAAAAACAaaaaaaaaaaaaa > 2:312451/1‑148 (MQ=9) tttgatGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg < 1:298485/150‑1 (MQ=38) atGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa < 2:437057/151‑1 (MQ=38) tGTCTGGCAGTTTATGGCGGGCGGCCT‑G‑CCCCCCCCCCC‑CCCGGCCGCTGCTTTTCAACCTTTAAAACCGCCCCCCGCGGGGTTTTTCCTCCCCGGGGGGGG‑GTTC‑CCCACCAAAAAAAGAAAAAAAAAAAAGGCCCGTCTTTCGGCTggg > 1:143315/1‑149 (MQ=2) gCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt < 2:598485/151‑1 (MQ=37) tGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt > 1:207602/1‑151 (MQ=38) ggcgggcgTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt < 1:19960/151‑1 (MQ=38) ggcgggcgTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt < 1:250417/150‑1 (MQ=38) cgTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc > 2:153364/1‑151 (MQ=39) cT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc > 2:414563/1‑151 (MQ=40) cGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTc < 1:419708/151‑1 (MQ=40) | AATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCT‑G‑CCCGCCACCCT‑CCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGT‑GTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTC > NC_000913/3423151‑3423436 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |