| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 48.7% | 39.2 / 16.8 | 39 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (16/4); new base A (0/19); total (16/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTG > NC_000913/436367‑436650 | tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc > 2:624662/1‑150 (MQ=255) tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCaa < 1:717862/150‑1 (MQ=255) cATGTTGAACCGTACCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑Accg > 2:431030/1‑151 (MQ=255) gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑Accgtcg < 2:328928/151‑1 (MQ=255) ggAACGGGCGCTGCGCTGGGCGCTCTCGGGCGGTGAAGTTTACGAGTGGTTTTTCACTGTGCCGAAACTAAACCGGTGCGCGTGGTATGTGGCTCCCGGCGCCCGGGGCGTACCGTTTCCCTGTTGCGGCAAAATG‑ACCGCCGATATCGaa < 2:202709/151‑1 (MQ=255) aaCGGGCGCTGCGCTGGGCGCTCTCTGGGGGTGAAGATTACGATTTTTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGGTGGTTGTGGCTCTCGGACCCGTGGGCGTCCCTTTTCCCTGTATGGGGAAAATG‑ACCGCCGATATCGAAgg < 1:27967/151‑1 (MQ=255) tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGttt > 2:427444/1‑151 (MQ=255) tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGttt > 2:684317/1‑151 (MQ=255) gcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGtttt > 2:680744/1‑151 (MQ=255) gcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGtttt > 1:582347/1‑151 (MQ=255) gctgGCCCCTCTCTGGAGGTGAGATTTAGGAGTTGTTTTTCGCTGTCACGGAGTTGAACCGGGACGCGTTGGCTGTGGCTCGGGGCCACGTGGGCGTACCTTTGCCTTGTATGGGGAAAGGG‑ACCGCCGATATCGAAGGGCTTTg < 2:708464/145‑1 (MQ=255) gcgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg > 2:49472/1‑151 (MQ=255) cgcTCTCTGGCGGTGACGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg > 2:293390/1‑150 (MQ=255) tctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg > 2:210316/1‑151 (MQ=255) tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa > 2:547117/1‑150 (MQ=255) tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa > 1:11048/1‑150 (MQ=255) gggggTGAAGATTAGGGTTTGTGTTTCACTGTCCCGAAACTAAACGGGGGCGGGCTGGATGGGGCTCCCGGACCCCGGGGTGTACTTTTTACTTGTAGGGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa < 2:7747/140‑1 (MQ=255) aaGATTAGGATTTGTTTTCCACTGTGCCGGAACTGACCCGGGGCGGGGGGGAGGGGGCTCTCGGACCCTGGGGCTTCCCGTTTCCCGGTATGGGAAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa < 1:427670/150‑1 (MQ=255) cgggtttttttttCCCTGGGCGGGAACTGAACCGGGGCGCGGTGTTTGTGCCTCTCGCCCACCGGGCCGTCCTGTTTCCCTGTTTGGGGCAAATG‑GCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa < 1:252865/141‑1 (MQ=255) gggttgttTTTCACTGGGCCGGAACGGAACCGGGGCGGGCGGGATGTGGCTCTGGGACACGGGGGCGTCCCTTTTCCCTGTATCGGGAAAAGG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa < 1:154039/147‑1 (MQ=255) tttCACTGTGCGGGAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCGTTTCCCTGTATCGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAGCCTGTTACATTAGACTGGaaa < 1:256300/148‑1 (MQ=255) tttCACGGTGCGGGAACTGACCCGTGGCGCGTTGGATGTGGCTCTGGGACCCCTGGGCTTCCCGTTTCCCTTTATCGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg < 2:687060/149‑1 (MQ=255) cACTGGGCCGAAACTGACCCGTGGCGCGCTGAATGTGCCTCTCGGACACTTGGGCGTCCGTTTTTCTTTTAGCGGGAAAAGT‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata < 2:605547/150‑1 (MQ=255) tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa > 2:54087/1‑150 (MQ=255) ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt > 1:626283/1‑151 (MQ=255) ggCGCGTGGGTTGGGGCTCTCGGACCCCGGGGCGTACCTTTTCCTTTTATCGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa < 1:454618/150‑1 (MQ=255) cgtgggttgtgGCTCCGGGACACGGGGGCGTCCTTTTTCCTTGTTTCGGGAAATTG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt < 1:199862/143‑1 (MQ=255) gggggttgtgGCTCTGGGACCCTGGGGCGTCCGTTTTCCCTGT‑TGGGGCAAATGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGAt < 2:317609/143‑1 (MQ=255) ggATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgc > 1:656882/1‑151 (MQ=255) cggaccctggggCGTCCCTTTTCCCTGTATCGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg < 1:355762/139‑1 (MQ=255) ggACCCGTGGGCGTCTCTTTTCCCTTTATGGGAAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg < 2:640842/147‑1 (MQ=255) cgccggggCGTCCTGTTTCCCTGTTTGGGAAAAAAG‑CCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa < 1:178773/145‑1 (MQ=255) cacGTGGGCGTCCCGTTTCCCTGTCGCGGGAAAATA‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaaaa < 2:532180/140‑1 (MQ=255) tGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg > 2:52269/1‑150 (MQ=255) gggCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg > 1:173001/1‑149 (MQ=255) gCGTCCCTTTTCCTTGTATGGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa < 1:360633/151‑1 (MQ=255) aCCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc > 1:95834/1‑151 (MQ=255) tttcctttttcggggAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg < 1:720232/138‑1 (MQ=255) tttACTTGTACCGGGAAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTa < 1:401033/143‑1 (MQ=255) aCCTGTGACGGGAAAAGG‑CCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCg < 1:236704/145‑1 (MQ=255) | TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATG‑ACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTG > NC_000913/436367‑436650 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |