| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,738,734 | 0 | T | G | 55.6% | 10.0 / 15.6 | 18 | W191G (TGG→GGG) | yiaN | 2,3‑diketo‑L‑gulonate transporter large permease YiaN |
| Reads supporting (aligned to +/- strand): ref base T (5/3); new base G (7/3); total (12/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTG > NZ_CP009273/3738648‑3738823 | ccAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGc < 1:55254/90‑1 (MQ=255) tATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTc < 1:149717/90‑1 (MQ=255) gTAAGCTGTTTATGGCGGGCCTTTCCCCCGGCATGATGATGGGCGCAACCCTGATGCTTACCTTGTGGGGGCAGGCCCGGCGCCCTAAtt > 1:376447/1‑90 (MQ=255) aaGCTGTTTATGGCGGGCATTTCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTTCCTGGGGGGGGCAGGCCCGTCGCCCTAAtttt > 2:201580/1‑89 (MQ=255) gTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGGCGCCTTTATTTAccccc > 1:431848/1‑89 (MQ=255) gTTTATGGCGGGCATTGCCCCCCGCATGATGATGGGGGCCACCCTGTTGCTTACCCGGGGGGGGGCGGGCGCGCGCCGTAATTTAccccc > 2:442446/1‑89 (MQ=255) ttATGGCGGGCATTGCCCCTTGCATGATTATGGGGGCAAACCCGGTGCTTACCTGGGGGGGGCAGGCCCGTCGCCCTTATTTACcccccc > 2:293726/1‑90 (MQ=255) ttATGGCGGGCATTGCCCCCGGCATGGTGATGGGGGCAACCCTGATGGTTTTCTGGGGGGGGCAGGCCAGTCGCCTTAATTTTCCCCGcc > 1:18945/1‑90 (MQ=255) cGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGGGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCCGCaaa > 1:27529/1‑90 (MQ=255) ggCATTTCCCCCGGCATGGTGATGGGGGCAACGCTGATGCTTTCCTGGGGGGGGGGGGGCCGGCGCCTTAATTTTTCCCCCCCGCaaaaa > 2:122711/1‑90 (MQ=255) ggCATTGCCCCCGGCATGATGATGGGGGCAACGCTGATGCTTACCTGGGGGGGGCAGGCCCGTCGCCTTAATTTACCCCCCCCGCaaaaa > 2:137762/1‑90 (MQ=255) tGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCa > 1:468933/1‑90 (MQ=255) cATGATGATGGACCCAACGCTGGTCTTTACCGGGTGGGGGCGGGCGCGTCATTTTAATTACCCCCGCCAGCAAAAAGCAACGATGCAGGa < 2:342006/90‑1 (MQ=255) atgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTgg > 1:369301/1‑90 (MQ=255) tGATGCTTACCTGGGGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtct < 2:4740/90‑1 (MQ=255) aTGTTTCCCTGGGGGGGGCGGGCCATTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTcc < 1:261379/90‑1 (MQ=255) cTTCCCTGGTGGGGGCAGGCCAGTCCCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGt < 2:376447/90‑1 (MQ=255) ttACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTa > 2:235090/1‑90 (MQ=255) tggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCAGTCTTTTGTCTCCGGTATCTGGGCGCTg > 1:358990/1‑90 (MQ=255) | CCAGCGGATTATCCATCAGTAAGCTGTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTG > NZ_CP009273/3738648‑3738823 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |