| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 497,028 | 0 | A | C | 55.6% | 29.0 / 37.5 | 27 | V556G (GTG→GGG) | ybaL | Kef family K(+) transporter |
| Reads supporting (aligned to +/- strand): ref base A (12/0); new base C (0/15); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.75e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACC > NZ_CP009273/496939‑497113 | gCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTca < 2:574349/90‑1 (MQ=255) ccTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAcc > 2:96402/1‑90 (MQ=255) cTGATGCGCCGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATATTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCCCAccc < 1:146428/90‑1 (MQ=255) gCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCg > 2:639685/1‑90 (MQ=255) ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt > 2:170548/1‑90 (MQ=255) ggCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGttt > 1:566140/1‑90 (MQ=255) gTCTTATCATCCCTACATTTTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcagc < 2:621707/90‑1 (MQ=255) ttttATCACCCCTACATTTTTTTCTTTTTTTACACCCGCAACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcagca < 1:96402/88‑1 (MQ=255) ttATCATCCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGt < 2:566140/90‑1 (MQ=255) aTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc > 1:141744/1‑90 (MQ=255) tCATGCCTAAATTTTTTTCATTTTTTCCCCCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTcc < 1:757608/90‑1 (MQ=255) ccgtGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTcc < 1:510694/87‑1 (MQ=255) cctcccTACATTTTTTTAATTTTTTACCCCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCa < 1:221157/86‑1 (MQ=255) gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa > 1:533081/1‑90 (MQ=255) gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa > 2:777110/1‑90 (MQ=255) ccTACTTTTTTTTCATTTTTTCCCCCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt < 1:151381/90‑1 (MQ=255) ccTACTTTTTTTTCATTTTTTACACCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt < 1:311138/90‑1 (MQ=255) ccTACATTTTTTTCATATTTTACCCCCGCCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt < 1:333885/90‑1 (MQ=255) aCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt > 1:338126/1‑90 (MQ=255) atttttttAATTTTTTACACCCGGCAACCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAc < 1:208219/87‑1 (MQ=255) tttttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTAcctg > 2:510694/1‑87 (MQ=255) ttttCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc < 2:466357/90‑1 (MQ=255) tttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg > 2:144837/1‑90 (MQ=255) ttttACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc < 2:763578/90‑1 (MQ=255) aCATCCGGCAACCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc < 2:226798/90‑1 (MQ=255) aTCCGGCACCCCCCTTTCCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc < 1:132445/90‑1 (MQ=255) ttACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt > 1:230124/1‑90 (MQ=255) cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc > 2:130929/1‑90 (MQ=255) | GCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACC > NZ_CP009273/496939‑497113 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |