Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,743,657 | 0 | A | C | 42.5% | 65.2 / 12.0 | 40 | T309P (ACC→CCC) | yiaK | 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent |
Reads supporting (aligned to +/- strand): ref base A (19/4); new base C (0/17); total (19/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.94e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCT > NC_000913/3743508‑3743807 | tttCACAAATTTTTATTGCCATTGAATTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTtactac > 1:403299/1‑151 (MQ=255) aTTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTgg < 1:29006/151‑1 (MQ=255) aTTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTgg > 1:750569/1‑151 (MQ=255) ttttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGcc > 1:326789/1‑151 (MQ=255) ttATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCg < 2:227941/150‑1 (MQ=255) tATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGa < 1:879152/150‑1 (MQ=255) aTTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGaaa < 2:222613/151‑1 (MQ=255) ccATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccg > 2:342770/1‑151 (MQ=255) ttGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAccgccg > 2:458450/1‑151 (MQ=255) gAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCTGAAAACCGCCGt > 1:86669/1‑150 (MQ=255) cGACGGCCCCCCCCGGGATGCAAAGCTACAACCCATCATGGTTTATGTTCTTGGTGCGGGGGGTGTTGACAAAAATCAGCCTTTTCGCTTCCCCGGCCATAAATTTACTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTtg < 2:661971/149‑1 (MQ=255) ggcccccccGGCGTTGCCAGCCTCCACCCCATCAGGAATTACTTTCTTGGTGCCGAGCGTCCTGACAAAAACCAGCCCTTCCGTTTCCCCGACCTAGATTTTACTCCCCTGCGGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACa < 1:452725/148‑1 (MQ=255) ccccccCGCGATGCCAAGCTGCAACGCATCATGGTTTCCGTTACTAGTGCCGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGGCCATGAATTTACTCCCCTGCGGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCg < 2:302252/151‑1 (MQ=255) ccACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgt > 2:801275/1‑151 (MQ=255) cACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtg > 2:464711/1‑151 (MQ=255) aCCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCgtg > 1:877351/1‑150 (MQ=255) cAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATc > 1:626251/1‑151 (MQ=255) aaGCTACAACGCATCAGGGTTTCCTTTACTGGTGCGGGGCGTGCTGACAAAAACCGGGCCTTCCGCTTCCCCGGCCTTGATTTTACTCCCGTGCGGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATc < 1:423739/150‑1 (MQ=255) aGCTGCAACGCATCAGGGATTACTTTCCTGGGGCCGGGCGTGCTGACAAAAACCAGGCCATCCGCTTCCCCGGCCATGATTTTACTCCCCTGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCaaaa < 2:599129/147‑1 (MQ=255) gCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg > 2:170742/1‑151 (MQ=255) atcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg > 1:552294/1‑149 (MQ=255) atcatGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg > 1:755034/1‑149 (MQ=255) tgtaTTACTTTACTGGCCCCGAGCGTGCTAACAAAAATCACCCCTTCCGCTTCCCCGCCCATAATTTTACTCCCCTGCTGGCCAAACACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATg < 2:61875/142‑1 (MQ=255) ggATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgagg > 1:529481/1‑147 (MQ=255) ttcctttcctggTGCCGGGCGTGCTAAAAAAAACCGGGCCATCCGCTCCCCCGGCAATAATTTTCCTCCCCGGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATa < 1:135926/138‑1 (MQ=255) ttACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATg > 2:534683/1‑151 (MQ=255) acggttgcggggCGTGCTGACAAAAAACGGCCCCTTCGCTCCCCCGGCCATAAATTTCCTCCCGTGCGGGCAGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAt < 2:185476/140‑1 (MQ=255) cggggCGTCCTGAAAAAAACCAGCCCTTTCGCTTCCCCGGCAATAATTTTCCTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtttt < 2:414608/143‑1 (MQ=255) aGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACat > 1:668168/1‑151 (MQ=255) aaaaaaaTCGGCCCTTTCCCTTCCCGGGCAATAAATTTCCTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAg < 1:751167/148‑1 (MQ=255) aaaaaTCAGGCCTTCCGCTTCCCGGGCTATAATTTTCCTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAt < 2:316004/129‑1 (MQ=255) aaaaccggccctttCGTTTCCCCGGCAATAAATTTACTCCCCGGCTGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATcg < 1:746811/131‑1 (MQ=255) aactctggtatttCGCTTCCGCGGTCATTATTTTATTCCCTTCCTGGCTAAAACCCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCa < 1:753943/135‑1 (MQ=255) aTTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCg > 2:489578/1‑151 (MQ=255) ttCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGt > 1:474913/1‑151 (MQ=255) tcccccGCCCATGTTTTTCCTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGttt < 1:348066/145‑1 (MQ=255) ccccGGCCTTAATTTTACTCCCTAGCCGCCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGt < 1:612361/138‑1 (MQ=255) ccccGGCCATGAATTTACTCCCCTGGGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccg < 2:750662/150‑1 (MQ=255) aCCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGccc > 2:748026/1‑150 (MQ=255) cccGGCAATAATTTTACTCCCCTGCGGGCAGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCccgc < 2:877469/151‑1 (MQ=255) attagtttgcgacCCTGGTGGGCGAAAAGCGCCGAAACCGGATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCAtt > 1:276156/12‑151 (MQ=255) ttACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGaaa > 1:244529/1‑149 (MQ=255) acCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGAGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCt > 2:93128/1‑151 (MQ=255) | TTTCACAAATTTTTATTGCCATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCGCAGCCCAATCCGTGCCGTTTGCCCGCCGCCATTGAAAAGGCGCT > NC_000913/3743508‑3743807 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |