| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,465,787 | 0 | G | A | 42.9% | 36.5 / 15.8 | 35 | R132C (CGC→TGC) | treB | trehalose‑specific PTS enzyme: IIB and IIC component |
| Reads supporting (aligned to +/- strand): ref base G (14/6); new base A (0/15); total (14/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.76e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTG > NC_000913/4465638‑4465936 | tttCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCg > 1:40796/1‑150 (MQ=255) cGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑g < 2:416024/151‑1 (MQ=255) gACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAg < 2:484070/149‑1 (MQ=255) gCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCaaa < 2:284160/151‑1 (MQ=255) ccgGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAAt > 1:255421/1‑151 (MQ=255) ggCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATc > 1:443801/1‑150 (MQ=255) aGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAAcg > 1:212297/1‑151 (MQ=255) accGCTTCCCCGATCAACCACAAAAAATCGTAGATCGTTGGAAGGGAGGGGTACTTTTGCGCCAGCGTCTGACCGTTGCTCGGGGGCAAATCGCCGACCAAATGGCGAAAACC‑GAGAATCAAACCGCCGCTAATCAACGCGGGCAGCAACg < 1:178023/149‑1 (MQ=255) ttCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGga > 2:146865/1‑149 (MQ=255) gcACCGACCAACCAAAGAAAATCGGAGATTGTTGGGAGGGAGGGGTACTTGTGCGCGGGCGTCTGCCCGTGGCTGGGGGGAAAAGCCCCGATCCCATAGCAAAAACG‑GGGCAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGgaag < 1:538945/149‑1 (MQ=255) cGATCAACCACAGAAAATCGTAGCTTGTTTGGAGGGAAGGGTACATTTGCGCCAGCTTCTGCCGGTTGCTCAGGGGCAAACGGCCGACCAGATGGCGAAAACG‑GAGGAAAAACCCGCCGCTAATCAACGCGGGCAGCAACGGGAagaaga < 1:63518/150‑1 (MQ=255) aaTCACCCACAGAAAATCTTAGTTCGTTGGAAGGGAGGGGTCCATTTCCGCCGGCGCCTGACCTTTCCTGAGGGGCAAATCGCCGGTCACATGGCAAAAACG‑GAGGAGCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATc < 1:444064/150‑1 (MQ=255) acaaaCCAAAAAAAACCGTAGTTCGTTTGGAGGGAGGGTTACTTTTGCCCCAGCTTCTGCCCGTTCCTCAGGGGAAACTCGCCAATCCCATGACAAAACCC‑GAGGAAAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAagaaga < 2:91618/145‑1 (MQ=255) tCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCAc > 2:115086/1‑151 (MQ=255) aaaaaagcttaactCGTTGGCAGGAAAGGTACCATTTGCCCCACCCTCTGACCGTTGCCCAGGGGAAAACCCCCACCCACATGCCAAAAACC‑GAGGATCAACCCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACcgc < 1:275814/133‑1 (MQ=255) aaaaaaTCGTAGTTCGTTTGAGGGGAAGGTTACTTTCGCCCCGGCGTCTGCCCTTTGCTCAGGGGAAAACCGCCAACCAAATTGGAAAACCC‑GAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGa < 2:249604/146‑1 (MQ=255) aaaaTCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATg > 2:421151/1‑150 (MQ=255) aaaTCGTAGTTCTTTTGCGGGAAAGGTTACATTTCCCCCAGCTTCTGCCCTTTGCTCATGGGCAAATCGCCGATCACATTGCAAAACCC‑GAGGAACAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATgag < 2:72994/151‑1 (MQ=255) aGATCGTTTGCGGGGAGGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCAGGGGCAAATCGCCGATCACATTGCAAAACCC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATc < 2:41668/149‑1 (MQ=255) gggAGGGGTACTTTTGCGCCAGCGTCTGCCCGTTGCTAATGGGAAAACCGCCAATCACATGGAAAAACCC‑GGGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCa < 2:190136/146‑1 (MQ=255) ttggccccAGCGTCTGCCCGTTGCTCAGGGGCAAACCCCCGATCAAATGGCAAAACCC‑GGGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCatat < 2:513140/141‑1 (MQ=255) ttGCGCCAGCGTCTGCCCGTTGCTCAGGGGAAAACCCCCGATCACATTGCAAAAACG‑GAGAATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCata < 2:400482/145‑1 (MQ=255) ccAGCGTCTGCCCGTTGCTCAGGGGAAAACCGCCGACCACATTGCAAAACCC‑GGGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggc < 2:415090/151‑1 (MQ=255) cAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGggcg > 1:174549/1‑151 (MQ=255) tGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttttt > 1:307813/1‑151 (MQ=255) ccGTTGCTCATGGGCAAATCGCCAATCACATTGCAAAACCC‑GAGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCttttt < 2:69879/148‑1 (MQ=255) ccGTTCCCCAGGGGCAAACCCCCGACCACA‑TGGGAAAACCGGGGGATAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCtttttt < 1:431347/149‑1 (MQ=255) cGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTAc > 1:109625/1‑151 (MQ=255) cATGGGAAAACCCCCGATCACATTGAAAAACCC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGctct < 2:340182/151‑1 (MQ=255) ggacAAATCGCCGATCACATTGCAAAACCC‑GAGGACAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGctct < 2:174609/145‑1 (MQ=255) tCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCt > 1:316900/1‑151 (MQ=255) cccccaaccacaTTGCAAAAACC‑GAAACCAAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCt < 1:115053/143‑1 (MQ=255) cGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTg > 2:153847/1‑151 (MQ=255) gATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGcc > 2:24808/1‑151 (MQ=255) tGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGt > 2:299295/1‑151 (MQ=255) cGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTg > 1:411369/1‑151 (MQ=255) | TTTCACCGCTGACCAGCAAATACCGACCGGCAGGTAGAAGAAGATCGCTTCACCGATCAACCACAGAAAATCGTAGATCGTTTGCAGGGAAGGGTACATTTGCGCCAGCGTCTGACCGTTGCTCATGGGCAAATCGCCGATCACATTGCGAAAACC‑GAGGATCAAACCGCCGCTAATCAACGCGGGCAGCAACGGGAAGAAGATCACCGCGAAATGAGAGATCAACTGCTCATGCCATTTCATATTATGCCGGGCGGCTTTTTTTACCTGCTCTTTATCAACCTGCGCCTGTCCGGTTG > NC_000913/4465638‑4465936 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |