| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 436,516 | 0 | C | A | 43.9% | 60.2 / 17.7 | 41 | Q295K (CAA→AAA) | thiL | thiamine monophosphate kinase |
| Reads supporting (aligned to +/- strand): ref base C (18/5); new base A (0/18); total (18/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.69e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT > NC_000913/436367‑436660 | tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc > 1:98161/1‑150 (MQ=255) ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCaa < 2:493510/149‑1 (MQ=255) cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt < 2:657213/149‑1 (MQ=255) gCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAc > 2:724935/1‑151 (MQ=255) aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg < 2:666661/150‑1 (MQ=255) aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc > 1:330408/1‑149 (MQ=255) ggCGCGGGGCGGGGCGCTCTCGGGGGGGAAAGTTTAGGATTTGTTTTTCCCTGTCCCGAACCTGAC‑CGGGGGCGCGCTGGATGTGGCTCTCGGACCCGGGGGCGTCCCGTTTACCTTTATCGGGAAAATGACCGCCGATATCGAAGGGCtt < 2:589602/151‑1 (MQ=255) cgctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTg > 2:506746/1‑151 (MQ=255) ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt > 1:449357/1‑151 (MQ=255) ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt > 1:496496/1‑151 (MQ=255) ctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt > 2:80140/1‑151 (MQ=255) tcgggCGGTAAAGATTAGGGGTTGTTTTTCTCGGTGCCGAAACGGAC‑CCGTGGCGCGCGGGATGTGGCTCTCGGACCCCGGGGCGTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg < 2:146746/145‑1 (MQ=255) ggTGAAGTTTACGAGTTGTGTTTCCCTGGGCCGGAACGAAC‑CCGGGGCGCGTTGGATGTGGCTCTCGGACCCGGGGGCGTACCTTTTCCCTGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg < 2:322604/142‑1 (MQ=255) ttACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt > 1:734840/1‑150 (MQ=255) ttACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt > 1:245145/1‑149 (MQ=255) tacggtttttttttCCCTGGGCGGAACCTGAC‑CCGTGGCGCGGTGGATGTGTCTCTCGGCCCCCGGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa < 2:685417/139‑1 (MQ=255) ttGTTTTTCCCTGTGCCGGAACGGAC‑CGGTGCCGCGCTGGGTGGGGCTCTCGGACCCCGGGGCTTACCTTTTCCCTGTTTCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt < 2:422275/148‑1 (MQ=255) tttCCCTGTGCCGGACCTGAC‑CGGTGGCGCGTTGGATGTGGCTCTCGGCCCCCTGGGCGTCCTTTTTCCCTGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg < 1:267200/149‑1 (MQ=255) ttCACTGTCCCGGAACTGAC‑CGGTGGCGCGCGGGATGTGGCTCCGGGCCCCCGGGGGGTACCGTTTCCCTTTATGGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa < 2:382978/150‑1 (MQ=255) tCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat > 2:600711/1‑150 (MQ=255) cACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat > 2:181331/1‑151 (MQ=255) cACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat > 1:265996/1‑151 (MQ=255) ctgggcagaaaCTGAC‑CCGTGGCCCGTGGGGTGTGGCTCTCGGACACGTGGGCGTCCCTTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata < 2:618814/139‑1 (MQ=255) gtgcGGAACTTAACCCGTGGCGCGCGTGATTGGCCTCTCGGACCCCTGGGCGTCCTTTTTCCTTGTTTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa < 2:636130/146‑1 (MQ=255) ggAATTGAA‑GCTTGGCGCGCGGGTTGTGGCTCTGGGACACGGGGGCGTACCGTTTTCTTGTCGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt < 2:591385/150‑1 (MQ=255) aaCTGAC‑CGGGGGCGCGCTGGATGTGGCTCTCGGACACGGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg < 2:269839/150‑1 (MQ=255) cctgggCGCGTTGGATGGGCCTCTCGGCCCCCGGGGCGTCCCTTTTCCCTTTCTGGGCAAAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc < 2:730215/143‑1 (MQ=255) cgcgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa > 1:632059/1‑148 (MQ=255) gggggtgggggCTCCCGCACCCGTGGGCGTCCCTTTTACCTGTATGGGGAAAGTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGa < 1:109864/139‑1 (MQ=255) tgtgGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGa > 2:630119/1‑150 (MQ=255) gcctctGGGACCGCGGGGCGTACTTTTTCCTTGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa < 1:436338/146‑1 (MQ=255) tcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc > 2:184621/1‑151 (MQ=255) ggACCCCGGGGCGTCCCTTTTCCTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg < 2:640846/147‑1 (MQ=255) gcccccgtGGGCGTCCCTTTTCCCTGTATGGGGAAAAGAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg < 1:438129/144‑1 (MQ=255) ggggCGTCCTGTTTACCTGTAGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa < 1:599539/150‑1 (MQ=255) ggggttactttttCCCTTTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGa < 2:508717/140‑1 (MQ=255) aCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc > 1:444593/1‑151 (MQ=255) tttCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg < 2:651647/150‑1 (MQ=255) ttACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc > 2:153058/1‑151 (MQ=255) tgaacggggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt < 2:473798/144‑1 (MQ=255) tGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt > 1:394491/1‑150 (MQ=255) aTCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt > 1:669217/1‑151 (MQ=255) | TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAA‑CCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTT > NC_000913/436367‑436660 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |