Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,745,509 | 0 | A | T | 41.7% | 40.6 / 13.2 | 37 | T236S (ACC→TCC) | yiaN | 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein |
Reads supporting (aligned to +/- strand): ref base A (5/16); new base T (15/0); total (21/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCT > NC_000913/3745361‑3745645 | gCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCCGCCCGTTTTcc > 2:24815/1‑151 (MQ=255) ccTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCCGCCCGTTTACCCCTAc > 2:389628/1‑151 (MQ=255) tggtggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCGGCCTGGTTTcccc > 2:538374/1‑146 (MQ=255) ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCgg < 2:24215/151‑1 (MQ=255) ggccgtacGCCTTATTTGCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGcc < 1:501347/144‑1 (MQ=255) gcccAGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGcc < 1:566898/149‑1 (MQ=255) gCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGACATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTGTCGCTTCGGCCTGGTTTTTCCTACCGCAGACGGGGCCGGtgct > 1:115185/1‑149 (MQ=255) aaTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCt < 2:169908/151‑1 (MQ=255) aaTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCCCCGGCCTGTTTATCCCCACCGAAGGCGGGGCCGTTGCTGCTTTCt > 2:207237/1‑151 (MQ=255) aTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCCGGCCTTTTTTTCCCCACCCAATACCCGGCGGCTGCCGCTTTCt > 1:237541/1‑150 (MQ=255) tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGGCT‑CGGGCTGGTTGTCCCTTACGAAGCCCGGGGCGGtggttgtt > 1:454934/1‑140 (MQ=255) ccccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCCGGCCTGTTTTCCCCTACCCAAGCCGGGGCCGTTGCTGCTTTTTATtcg > 1:196152/1‑150 (MQ=255) cGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTCTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGt < 2:542664/151‑1 (MQ=255) gCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGtt < 1:183858/151‑1 (MQ=255) ccAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCCGGCCTGTTTTCCCCCACCCAAACCGGGGCCGTTGCTGCTTTTTATGCGCTGttt > 1:443544/1‑151 (MQ=255) cAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTTCCCCCACCGAAGCCGGGGGCGTTGCTGCTTTTTATGCGCTGttt > 2:465511/1‑147 (MQ=255) cAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCGGCCTGTTTTTCCCCACCGAAGCCGGGGCCGTTGCTGCTTTTTTTGCGCTGGTTTTCGcccc > 1:477995/1‑146 (MQ=255) gATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTAtt < 2:116745/151‑1 (MQ=255) aaaTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCCGCCTGTTTTCCCCCACCGAAACCGGGGGCGTTGGTGCTTTCTTTGCGCTGTTTTTTGCCCCCGTTTTTTAcc > 1:208961/1‑147 (MQ=255) aTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCGGCCCGTTTTCCCCTCCCCAAGCCCGGGCCGGTGCTGCCTTCTTTTCGCTGTTTGTCGCCCCAGTTTTTTTCCGCGaa > 1:458293/1‑150 (MQ=255) tCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTTTTCCGGCCTGTTTAACCCTACCTAAGGCAGGGACGGTGGTGCTTTTTATGCGCTGGTTTTCGGCACAGGTTTTTATCGttgaa > 2:35860/1‑146 (MQ=255) ggCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGaaa < 2:26106/147‑1 (MQ=255) cACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTTCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTTTCGCCACAGTTATTTTCCGTGAAATGAcc > 2:275865/1‑150 (MQ=255) tttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGc < 2:111448/149‑1 (MQ=255) ctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCTCTCCGGCCTGTTTATCCCCACCCAAACCGGGGGCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCCCCCTTTTTTACCGTGAAATGACCTTCGCCccac > 1:75238/1‑149 (MQ=255) gTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTTGCTCCGGCCTGTTTTTCCCCACCGAAACCGGGGGCGGTGGTGCTTTTTATGCGCTGTTTTTTGCCACACTTATTTTCCGTGGAAAGACCTTTGCCCCACttttggt > 1:4687/1‑147 (MQ=255) aTCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCGGCCCGTTTTTCCCCACCCAAACCGGGGCCGGTGCTGCCTTTTATTCGCTGTTTTTTGCCACAGTTATTTACCCTGGAAAGACCTTCGCCACACTTTGGcct > 2:194471/1‑148 (MQ=255) tCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGGTTTGCTCCGGCCTGTTTTCCCCCCCCGAAACCCGGGGCGGTGCTGCTTTCTTTGCGCTGTTTTTCGCCACAGTTATTTTTCGCGAAAAGACCTTCGCCACACTTTTGGATGt > 1:130177/1‑151 (MQ=255) gggCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTAc < 1:465418/150‑1 (MQ=255) gcgcTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCa < 1:210374/151‑1 (MQ=255) gTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTg < 1:70411/149‑1 (MQ=255) tcttcCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGgcgc < 1:117122/150‑1 (MQ=255) tcCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGcgcg < 1:539792/148‑1 (MQ=255) aGTCATCATTATTGGC‑GGGTTTTGCTCCGGCCCGTTTTTCCCTACCCCAGCCGGGGGCGGTTCTGCTTTCTATGCGCTGTTTTTTGCCACAGTTATTTTCCGCGAAATGACCCTCGCCCCCCTTTTGGCTGTACTCATTTGCgcgggg > 1:338847/1‑146 (MQ=255) ggC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTg < 1:240521/150‑1 (MQ=255) ttCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGtt > 2:111999/1‑151 (MQ=255) cTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCt < 2:40191/150‑1 (MQ=255) | GCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGC‑GGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCT > NC_000913/3745361‑3745645 |
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: ‑ |