Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 497,020 | 0 | A | C | 41.4% | 19.9 / 10.1 | 29 | *559E (TAA→GAA) | ybaL | Kef family K(+) transporter |
Reads supporting (aligned to +/- strand): ref base A (14/3); new base C (0/12); total (14/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.05e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGA > NZ_CP009273/496931‑497105 | ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTa < 1:662013/90‑1 (MQ=255) cATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcc < 1:1483444/90‑1 (MQ=255) ctctCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcac < 1:412363/90‑1 (MQ=255) tctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc < 2:887882/90‑1 (MQ=255) gATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGAAACCACCGTTTACCCCGTCACCACCTCACccgc > 1:366308/1‑90 (MQ=255) tGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgccg > 1:1475518/1‑90 (MQ=255) aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTgg > 1:160917/1‑90 (MQ=255) aCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTgg > 1:800011/1‑90 (MQ=255) ttCTTATCATCCCTACATTTTTTTCATATTTTACCCCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCcagc < 2:27969/89‑1 (MQ=255) gTCTTATCATGCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc < 2:381950/90‑1 (MQ=255) ttATCATCCCTACATTTTTTTAATTTTTTACACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGt < 2:1475518/90‑1 (MQ=255) aTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc > 1:1421276/1‑90 (MQ=255) aTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc > 1:1161161/1‑90 (MQ=255) gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa > 2:505850/1‑90 (MQ=255) tACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATgg > 1:394613/1‑90 (MQ=255) atttttttAATTTTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc < 2:1112125/87‑1 (MQ=255) tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg > 2:1251032/1‑90 (MQ=255) ttttCATTTTTTCCACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc < 2:1213782/90‑1 (MQ=255) ttttAATTTTTTACACCCGGCAACCACCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc < 2:440742/90‑1 (MQ=255) cATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt > 1:69406/1‑90 (MQ=255) cATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt > 1:333918/1‑90 (MQ=255) atTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCa > 1:403300/1‑90 (MQ=255) cccccGCCACCCACCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg < 1:135478/87‑1 (MQ=255) cgggCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc < 2:160917/88‑1 (MQ=255) ccGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc < 2:1252457/90‑1 (MQ=255) ccGGCACCCCCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc < 2:800011/90‑1 (MQ=255) cGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCa > 1:267883/1‑90 (MQ=255) gCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc < 2:1421276/90‑1 (MQ=255) acccaccTTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc < 2:394613/88‑1 (MQ=255) gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGa > 2:1143751/1‑90 (MQ=255) | GGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGA > NZ_CP009273/496931‑497105 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |