Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,745,397 | 0 | A | C | 40.5% | 42.1 / 15.5 | 42 | L198F (TTA→TTC) | yiaN | 2,3‑diketo‑L‑gulonate:Na(+) symporter ‑ membrane subunit |
Reads supporting (aligned to +/- strand): ref base A (17/8); new base C (0/17); total (17/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.41e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCG > NC_000913/3745260‑3745530 | tCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTTGCATGATGATTGGGGCAACGCTGATGCTTTCCTGGTGGGGGCGGGCCAGGCGCCTTAATTTTcc > 2:392895/1‑140 (MQ=255) tttATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGGGGCAGGCCAGTCGCCTTAATTTACCCCCCCAGCa > 2:658810/1‑140 (MQ=255) tttATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCa < 1:979227/140‑1 (MQ=255) cTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTTGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCaa < 1:648850/140‑1 (MQ=255) cTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCaa < 2:268068/140‑1 (MQ=255) cTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCaa < 2:678879/140‑1 (MQ=255) ttCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTTCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGGCAGGCCCGTCGCCTTAATTTACCCCCCCCGCAAAAAAGAAAGAAGCAgg > 2:91946/1‑140 (MQ=255) ttCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGGGCAACGCTGATGCTTACCCGGGGGGGGGAGGCCAGTCGCCTTTATTTTCCCCCCCAGCAAAAAACAAAGATGCAgg > 2:1000245/1‑140 (MQ=255) tCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGGGCAACGCTGATGCTTACCTTGTGGTGGCAGGCCCGTCGCCTTAATTTACCCCCCCAGCAAAAAACAACGATGCAGGa > 2:648850/1‑140 (MQ=255) ggATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTg < 1:115191/140‑1 (MQ=255) gATTATCCAGCAGTAAGCTGTTTAGGGGGGGCATTCCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACTTGGTGGTGGCAGCCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTgg < 2:757536/140‑1 (MQ=255) aTTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGc < 1:392895/140‑1 (MQ=255) aTTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGc < 1:952294/140‑1 (MQ=255) tATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTTCCTGGGGGTGGCAGGCCCGTCGCCTTAATTTACCCCGCCCGCCAAAAGCAACGATGCAGGAAATCTGGCAc > 1:644800/1‑140 (MQ=255) tCAGTAAGCTGTTTATGGCGGGCGTTGCCCCTGGCATGATGATGGGCGCAACGCTGATGTTTACCTGGTGGGGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt < 2:824350/140‑1 (MQ=255) cAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTg < 1:658810/140‑1 (MQ=255) aGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGTAGAAGGGCGGAACGCTGATGCTTTCCTGGGGGTGGCAGGGCCGGCGCCTTTATTTTCCCCCCCAGCAAAAAACAACGAAGCAGGAAAACTGGGACTCtttttt > 2:952294/1‑138 (MQ=255) tAAGCTGTTTATGGCGGGCATTGCCCCTGGCGTGATGATTGGCGCCACGCTGATGCTTTCCTGGGGGGGGCCGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCACGGAATCTGGGACTCTTTTTtct > 2:4934/1‑140 (MQ=255) gacaggTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtctc < 1:834696/135‑1 (MQ=255) aGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt < 2:903515/134‑1 (MQ=255) aGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGGGGCAGGCCAGGCGCCTTTATTTACCCCCCCCGCAAAAAACAAAGATGCAGGAAAACTGGCACCCTTTTTTCTcc > 2:749721/1‑140 (MQ=255) aGCTGTTTATGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtttt > 1:903515/1‑134 (MQ=255) tGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGt < 2:878458/133‑1 (MQ=255) tGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGt > 1:878458/1‑133 (MQ=255) ggCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCt < 1:595469/140‑1 (MQ=255) gCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCCGCAAAAAACAACGATGCAGGAAAACTGGCACTCTTTTTTCTCCCGTATCTGGgctgtc > 2:834696/1‑136 (MQ=255) cTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCa < 2:535787/140‑1 (MQ=255) gCATGTTGTGGGGCGCAACGTTGGTGCTTCCGTGGTGGTGGCAGGCGAGTCGCCTTAATTTTCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGtc < 2:661580/140‑1 (MQ=255) atgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTcatcat > 1:54990/1‑140 (MQ=255) tgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCattatt > 1:632258/1‑140 (MQ=255) ggCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTgg > 1:495281/1‑140 (MQ=255) gcgcAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGc < 1:4934/140‑1 (MQ=255) gcgcAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGc > 1:583364/1‑140 (MQ=255) cGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGttt > 2:138105/1‑140 (MQ=255) cGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGttt > 2:652637/1‑140 (MQ=255) aTGCTTACCTGGTGGTGGCAGGCCAGTCGCTTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCGCTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTgg < 2:924636/128‑1 (MQ=255) aTGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCGCTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTgg > 1:924636/1‑128 (MQ=255) cTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg < 2:54990/140‑1 (MQ=255) ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGc < 1:138105/140‑1 (MQ=255) ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGc < 1:652637/140‑1 (MQ=255) tACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGcc > 1:660241/1‑140 (MQ=255) ccTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTccctcctt > 2:39662/1‑134 (MQ=255) ccTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCTCCCGccac > 2:382321/1‑138 (MQ=255) tggtggtggCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt < 2:660241/140‑1 (MQ=255) ggCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg < 1:805504/127‑1 (MQ=255) ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCgg > 2:805504/1‑127 (MQ=255) ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt < 1:515329/133‑1 (MQ=255) ggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGtt > 2:515329/1‑133 (MQ=255) gCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGGCCCGTTTTTCCCtc > 2:295771/1‑139 (MQ=255) cAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTTCTCCGGGCTGTTTTCCCCCAc > 1:285781/1‑140 (MQ=255) cTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGc < 1:295771/140‑1 (MQ=255) tAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCg < 2:495281/140‑1 (MQ=255) | TCCATTCCGTTTATTATCTTCGGCGTTTCCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCG > NC_000913/3745260‑3745530 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |