| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,369,516 | 0 | T | A | 44.8% | 29.1 / 10.5 | 30 | S42T (TCA→ACA) | pspD | peripheral inner membrane phage‑shock protein |
| Reads supporting (aligned to +/- strand): ref base T (1/15); new base A (13/0); total (14/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.06e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAAT > NC_000913/1369377‑1369664 | aaCTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGGGGCGGGCTGGGGGAT‑AAAAAAAGTTTCt > 2:51892/1‑148 (MQ=255) gtgAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCg < 1:577134/149‑1 (MQ=255) ggAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGG‑AAAAACAGTTTGCCGCCGAcc > 2:207130/1‑148 (MQ=255) aaaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGT‑AAAAACAGGTGGTCGCCCAcc > 1:372902/1‑146 (MQ=255) aaaGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGa < 2:139895/151‑1 (MQ=255) aGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGGT‑AAAAATAGGTGCTCGcc > 1:642198/1‑140 (MQ=255) aaCAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGTAAAAATAAGTT‑CTCGCCGGCCGCCGGAAAAGGTTGTtggtggtgg > 1:98655/1‑148 (MQ=255) aGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGT‑AAAAAAAGTTGCTCCCCCACCG‑CTGGAAAAGTTGCTGGGTGTGGGGCTgg > 1:53235/1‑150 (MQ=255) aaaaGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGtt < 1:531946/151‑1 (MQ=255) tAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAC‑AAAAACAAGTGCTCGGCCCCCC‑CTGGAAAAGTTTGGGGGTGGGGGGCTGGGAACGCCGttta > 2:116934/1‑145 (MQ=255) aaaGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCgg < 1:597300/151‑1 (MQ=255) ggTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCta < 2:661342/151‑1 (MQ=255) gTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGT‑AAAAAAAGTTGCTCGCCGACCG‑CCGAAAATGTTGCTGGGTGTGGGGCTGGGAACCCTGTTAAATCggg > 1:47835/1‑144 (MQ=255) tttCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGCG‑AAAAAAAGTTGGCCCCCGGCCG‑CCTGAAATGTTGGTGGCTGTGGGGCTGGAAACCCTGtttaaacggggcg > 2:198106/1‑135 (MQ=255) tttCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGGT‑AAAAACAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGGTGTGGCGCTGGGACCGCCGTTTAGTCGGGCTGCtaa > 1:128762/1‑150 (MQ=255) aaTTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGGGGT‑AAAAACAGATGCTCGCCGACCC‑CTGGAAATGTTTCTGGGTGTGGCGCTGGAAACCCTGTTTAATCGGGGTGGTAATaaa > 1:80413/1‑149 (MQ=255) gCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGcac < 1:65691/151‑1 (MQ=255) cAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGTT‑GAAAACATGTGGTCGCCCGCCG‑CCGGAAAAGTTTGTGGGTGGGGCGGTGGGAACGCCGGTTAGGACGGCTGCTCATaaaaaggg > 2:536218/1‑143 (MQ=255) ggCAAGCTGGTCCTTCTCACCTCACTGCGCTACGGCCCGGCGGGTTTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGTACAGc < 2:175611/151‑1 (MQ=255) gTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGGTGGGGGGT‑AAAAACAGTTGGTCCCCCCCCC‑CCGGAAAAGTTGGTGGGTGTGGGGGTGGGACCGCCGTTTAATTGGGCTGGTTAAAAACATGGCCCAggtttttaa > 2:210067/1‑142 (MQ=255) gTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATaaaa < 1:616175/151‑1 (MQ=255) cTGCGCTATGGCCCGGCGGGTGTGGCGGGCCGGGGGGT‑AAAAACAGTTGCTCGCCCGCCG‑CCGGAAAATTTGTTGGGTGGGGGGCTGGGACCGCTGTTTAAGCGGGGTGCTAATAAAA‑TGGGACAGCGGTTTAAAAAgggggggggggtt > 1:378663/1‑145 (MQ=255) gcacTATGGCCCGTCGGGTGTGGCTGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTc < 2:295998/146‑1 (MQ=255) cgcTATGGCCCGGCGGGTGTGGCTGGCTGGGCGAG‑AAAAGCGGTTGCTCAGCGGCCG‑CTGAAAATGTTATTGGGTTTGGCGCTGTAACCGCTGTTAAGTCGTGCTGCTAACGATC‑TGGCACAGCTTTATAAAAGGTGAGGGGAATTGCGCa < 2:357232/151‑1 (MQ=255) cgcTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCa < 1:559185/151‑1 (MQ=255) ggCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAAtt < 2:609819/151‑1 (MQ=255) gCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTg < 2:398812/151‑1 (MQ=255) gtgGCGGGCTGGGCGGT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCa < 2:651335/150‑1 (MQ=255) gtgGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAg < 1:102652/151‑1 (MQ=255) aaTCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAAt < 1:548835/151‑1 (MQ=255) | AACTGTGAGGAAAGTTATGAATACTCGCTGGCAACAGGCCGGGCAAAAGGTAAAGCCTGGTTTCAAATTAGCAGGCAAGCTGGTACTTCTTACCGCACTGCGCTATGGCCCGGCGGGTGTGGCGGGCTGGGCGAT‑AAAATCAGTTGCTCGCCGACCG‑CTGAAAATGTTGCTGGCTGTGGCGCTGGAACCGCTGTTAAGTCGGGCTGCTAATAAAC‑TGGCACAGCGTTATAAAAGGTGAGGGGAGTTTCGCAAAAATTGTTAAATCTCAGGCGTATAATGGATGGCAAT > NC_000913/1369377‑1369664 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |