| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,880 | 0 | T | C | 58.1% | 8.4 / 30.7 | 31 | K271E (AAA→GAA) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base T (11/2); new base C (0/18); total (11/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.21e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-05 | |||||||||||
TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT > NC_000913/357732‑358027 | tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc > 2:340490/1‑150 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:720012/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg > 2:64535/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg > 2:413159/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTcc < 2:606931/148‑1 (MQ=255) atcccgcgCCGCAGACCTGTGCCGGCGGCTCCCGTGGTTCTTTCCCGCCAGCACCCGCGCTTTTACCCACGGCCAGGTGAAAAAGTCCCTCGGCGGCGGCTGTCCGCCATCTATTTTCCCCCGGCAACAAAACCGCCGTTCTCTCCagtag < 1:533056/144‑1 (MQ=255) aaTTCCAGCCGAAGACCTGTGCCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag < 2:726786/151‑1 (MQ=21) ccGCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCAAACACCCGCGCATTTACCCAACGCCTTGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg < 2:41479/151‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:688422/1‑151 (MQ=255) cccTGTGCCAGCGGCTCCCGTATTTCTTTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTACAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g < 1:114210/150‑1 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑Ggaaa > 1:582513/1‑150 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAA‑Gaaaa > 2:508485/1‑150 (MQ=255) ggtCTGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑Ggaaaa > 2:85445/2‑149 (MQ=255) gCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTCAAAAAAGGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCtt < 2:70546/150‑1 (MQ=25) ttCATTCCCGCCAACAACCGCGCATTTATCCACCGCCATGTGCAAAAATGCCTCCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgt < 2:250875/151‑1 (MQ=255) cATTCCCGCCAACACCCGCGCTTTTACCAAACGCAATGTAAAAAAAGCCCTCCGGGGCGGTTGTCTGCCGGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgttg < 1:320783/151‑1 (MQ=255) ccAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATCCTTCCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa < 2:691201/151‑1 (MQ=255) caaacaccCGCGCATTTCCCAAACCCCATGCGAAAAAACTCCTTGGCGGCGCCTGTAGCCCCGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa < 1:138579/149‑1 (MQ=255) caaacaccCGCGATTTTCCCAAACCCCTTGTAAAAAAATGCCTCCGGGGGGGCTTTCTGCGAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa < 1:426750/149‑1 (MQ=255) caacacCCGCGCTTTTACCCACCCCCATGGGAAAAAATGCCTCCGCGGCGGCTGCCTGCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 1:88149/151‑1 (MQ=255) cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGAGGAAAAACCTTTAAACCC‑CCCTGTTTTGTGGGGATTGCCGCT‑GGtt > 2:602988/1‑150 (MQ=25) tttACCCAACCCCTTGTGAAAAAACGCCTTCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAAc < 2:295957/151‑1 (MQ=255) gtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GGAAAAAATTTAAACCCCCCCTGTTGTTGGGGGATCTCTGCT‑GGTGACAAATTGGGAAgggga > 1:403562/1‑149 (MQ=255) gtaaaaaaaTGCCTTCGGGGCGCCTTTCGGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa < 1:544741/146‑1 (MQ=255) tgaaaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAg < 2:199253/148‑1 (MQ=255) tgCAAAAATGCCTGGGGGGCGGCGTTCTGCAAGTTTTTTTTCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAg < 1:708438/151‑1 (MQ=255) cggcggCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAACTTTTAACCCCCCCTGTTTTGTGGGAATCGCCGCTGGGGAAAAAATGGGAAAGGGAAggggggggatcaacc > 2:287571/1‑136 (MQ=255) gcggcTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc < 2:387515/150‑1 (MQ=255) ttGTCAGCAAGTTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:65304/150‑1 (MQ=255) cTGTCTGCAATTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:327374/151‑1 (MQ=255) cTGTCTCCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATATTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:255432/151‑1 (MQ=255) tttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:413572/151‑1 (MQ=255) tttcCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 1:211235/151‑1 (MQ=255) | TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT > NC_000913/357732‑358027 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |