Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,737,017 | 0 | A | C | 52.2% | 17.2 / 36.4 | 23 | T309P (ACC→CCC) | yiaK | 3‑dehydro‑L‑gulonate 2‑dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (11/0); new base C (0/12); total (11/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.40e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCG > NZ_CP009273/3736888‑3737123 | aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTgg > 1:192395/1‑139 (MQ=255) aTTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTgg > 2:15995/1‑139 (MQ=255) gAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTCCCCTGCTGGcc < 1:560815/138‑1 (MQ=255) gACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAAcc > 2:287442/1‑139 (MQ=255) gCTTACCGACGGTCCCCCCCGCGATGCCAAGCTGAAACGCATCATGAATTACTTTACTGGGGCCGGGCGGGCTGACAAAAACCGGGCCATCCGCTTCCCCGGCAATAAATTTCCTCCCCTGCTGGCAAAAAACCGCCGt < 1:131170/139‑1 (MQ=255) cTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGt > 2:196930/1‑138 (MQ=255) cTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGt > 2:131170/1‑138 (MQ=255) taccGACGGTCCCCCCCGCGATGCCAAGCTGCAACGCATCATGGTTTCCGTTACTAGTGCCGGGCGTGCTGACAAAAATCAGGCCATTCGCTTCCCCGGCCATGATTTTACTCCCCTGCTGGCCGAAAACCGCCGTaa < 1:452901/135‑1 (MQ=255) aCGGTCCCCCCCGCGACGCCAAGCTGCAACGCATCATGGATTGCGTTACTGGCGCGGGGCGGGATGACAAAAATCGGGCCATTCTCTTACCCGGCAATAATTTTCCTCCCCGGTTGGCAGAAAGCCGCCGTAACGGCAt < 1:577360/139‑1 (MQ=255) ggTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCa > 2:135815/1‑139 (MQ=255) tCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACt > 2:60589/1‑139 (MQ=255) gATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTtgatga > 1:112841/1‑137 (MQ=255) catcatGGATTACTTTACTGTTGCCGGGCGTGCTGACAAAAACCAGCCCATCCGCTTCCCCGGCCATAATTTTACTCCCCTGCGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATc < 2:192395/139‑1 (MQ=255) atGGATTACTTTACTAGTGCCGAGCGTGCTGACAAAAACCAGCCCATTCGCTTCCCCGGCCATAATTTTACCCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAg < 1:135815/138‑1 (MQ=255) actttACTGGTGCCGGCCGGGCTGACAAAAAACAGCCCATCCGCTTCCCCGGCCTTGATTTTCCTCCCCGGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTAt < 2:112421/135‑1 (MQ=255) cTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATgag > 1:302057/1‑135 (MQ=255) gggCGTGCTGACAAAAACCAGCCCATCCCCTTCCCCGCCCATAATTTTCCTCCCCTGCGGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCAt < 1:287442/137‑1 (MQ=255) gggCGTGCTGAAAAAAACCGGGCATTCCGCTCCCCCGGCCATAATTTTACTCCCCGGTGGGCAAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCAt < 2:302057/137‑1 (MQ=255) aGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCAt > 2:369598/1‑138 (MQ=255) gCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGAtttt > 2:136367/1‑139 (MQ=255) aaaCCAGCCCTTCCGCTTCCCCGGCCATAATTTTACTCCCCTGCTGCCGAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACata < 1:196930/139‑1 (MQ=255) cAGCCCATCCCCTTACCCGGCAATAATTTTCCTCCCCGGCTGGCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATcgc < 2:457335/139‑1 (MQ=255) agcccATCCGCTTCCCCGGCCATAATTTTACTCCCCTGCTGCCCAAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATcgcg < 1:60589/136‑1 (MQ=255) | ATTGAAGTGGACAAGCTTATCGACGGTCCCACCCGCGATGCCAAGCTGCAACGCATCATGGATTACGTTACTAGTGCCGAGCGTGCTGACGAAAATCAGGCCATTCGCTTACCCGGCCATGAATTTACTACCCTGCTGGCCGAAAACCGCCGTAACGGCATCACTGTTGATGACAGCGTGTGGGCCAAAATCCAGGCGTTATGAGGAGATAAGTCATGATTTTTGGACATATCGCG > NZ_CP009273/3736888‑3737123 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |