Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,738,757 | 0 | A | C | 55.6% | 11.1 / 19.2 | 21 | L198F (TTA→TTC) | yiaN | 2,3‑diketo‑L‑gulonate transporter large permease YiaN |
Reads supporting (aligned to +/- strand): ref base A (6/2); new base C (0/10); total (9/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.51e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCA > NZ_CP009273/3738673‑3738836 | gTTTATGGCGGGCATTGCCCATGGCATGATGATGGGGGCAACGCTGATGCTTTCCTGGTGGGGGCAGGCCAGGCGCCTTAATTTTCCCCg > 1:187280/1‑90 (MQ=255) ttATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGcc < 1:86881/90‑1 (MQ=255) tGGCGGGCATTGCCCCCGGCATGATGATGGGGGCAACGCTGATGCTTACCTGGTGGGGGGAGGCCAGGCGCCTTAATTTTCCCCCCCCGc > 1:364505/1‑90 (MQ=255) tGGCGGGCATTGCCCCCGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGGGGCAGGCCAGTCGCCTTAATTTTCccccccccc > 1:314701/1‑87 (MQ=255) gatgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAAc > 2:294253/1‑78 (MQ=255) gatgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAAc < 1:294253/78‑1 (MQ=255) gatgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAAt > 2:279642/1‑90 (MQ=255) gatgatgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAAt > 2:264273/1‑90 (MQ=255) tgGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGc > 1:189750/1‑90 (MQ=255) ggCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCAc < 1:270352/90‑1 (MQ=255) ccaCGCTGATGCTACCCGGGGGGGGGGGGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCtt < 2:359624/88‑1 (MQ=255) gCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGt < 1:68905/90‑1 (MQ=255) cTGATGCTTCCCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtc < 1:101072/90‑1 (MQ=255) cTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtc < 1:312264/90‑1 (MQ=255) gATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGtctc > 2:118711/1‑90 (MQ=255) aTGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTcc < 1:279673/90‑1 (MQ=255) gCTTCCCTGGTGGTGGCAGGCCAGTCCCCTTATTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCgg < 2:187280/90‑1 (MQ=255) tACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACCCTTTTGTCTCCGGTAt > 2:259064/1‑90 (MQ=255) tggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTg < 1:264273/90‑1 (MQ=255) tggCAGCCCAGTCGCTTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTg < 2:100643/90‑1 (MQ=255) gCCTTAATTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGtca < 1:279642/90‑1 (MQ=255) | GTTTATGGCGGGCATTGCCCCTGGCATGATGATGGGCGCAACGCTGATGCTTACCTGGTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCA > NZ_CP009273/3738673‑3738836 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |