Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 172,755 | 0 | A | G | 53.6% | 16.3 / 31.2 | 29 | T388A (ACG→GCG) | clcA | H(+)/Cl(‑) exchange transporter ClcA |
Reads supporting (aligned to +/- strand): ref base A (2/11); new base G (15/0); total (18/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.63e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGCTGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAC > NZ_CP009273/172670‑172840 | gATGCTGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGtt < 1:447245/90‑1 (MQ=255) cTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGTGGCTTGTTATTTcc > 2:741569/1‑90 (MQ=255) cTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGcc < 1:225943/90‑1 (MQ=255) gTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGaa < 2:304999/90‑1 (MQ=255) aCTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATg < 2:424881/90‑1 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTGTTTTGCCGGGATggggggc > 1:269460/1‑88 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGTTATTGCCGGGATgggggga > 2:120965/1‑88 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGTTTTGCATTGGCGGGATGGGGGGAtt > 2:620369/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTTCTATTTTCGGGATGGGGGGAtt > 2:501240/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGGAtt > 1:781666/1‑90 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTGCTTTTTACGGAATGGGGGCATTACttg > 2:461707/1‑88 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTGCCGGGATGGGGGGATTACttg > 1:45566/1‑88 (MQ=255) gCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGTTTGTTTTTTTCGGGAAGGGGGGGTTAATggcg > 2:57331/1‑90 (MQ=255) gCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTggcg < 1:676469/90‑1 (MQ=255) aaTGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTCTATTGCCGGGATGGGGGCATTAATGGCGGGATCTAt > 1:16507/1‑90 (MQ=255) ggTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTAtttg > 1:334531/1‑88 (MQ=255) tGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTGCTTTTGCCGGGATGGGGGGATTACTTGCGGCATCTATTcgcgc > 2:49068/1‑90 (MQ=255) gTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGTTTTTTATTGTCGGGATGGGGGGGTTACTGGAGGGGTCTATTTGCGGGCCg > 2:435093/1‑90 (MQ=255) tGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATgggg < 2:352025/58‑1 (MQ=255) tGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATgggg > 1:352025/1‑58 (MQ=255) gTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGGATTACTTGCGGGGTCTATTTGCGCGCCGCTAACggg > 2:9559/1‑90 (MQ=255) gTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGCATTACTGGCGGGATCTATTCGCGCGCCCTTAACggg > 1:447629/1‑90 (MQ=255) ccGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATc < 1:405158/90‑1 (MQ=255) ccGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATc < 2:781666/90‑1 (MQ=255) aaTATCCCCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttc < 1:120965/90‑1 (MQ=255) atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct < 1:597489/90‑1 (MQ=255) tatCACCTTGAGGCGGGGGCGTTTGCTATTTCCGGGATGGGGGGATTACTGGCGGGCTCTATTTTCGCGCCCGTTACGGGGGGCAttttt > 2:330173/1‑87 (MQ=255) tatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttctg < 1:435093/90‑1 (MQ=255) ggggGCGTTTTCTATTGCCGGAATGGGGGCGTTACTGGCGGGCTCTATTTGCGCGCCGTTTAGGGGGGTGATTCTTGTTTTGGTGATggc > 2:342629/1‑88 (MQ=255) | GATGCTGGCGCTGGGTACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAC > NZ_CP009273/172670‑172840 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |