Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 172,772 | 0 | A | G | 56.5% | 17.5 / 26.1 | 23 | G393G (GGA→GGG) | clcA | H(+)/Cl(‑) exchange transporter ClcA |
Reads supporting (aligned to +/- strand): ref base A (3/7); new base G (13/0); total (16/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.89e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATA > NZ_CP009273/172686‑172845 | aCTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATg < 2:86116/90‑1 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTGCTATTTCCGGGATgggg > 2:206338/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATtggg > 1:191378/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTGCTTTTTCCGGGATGGGGGGAtt > 2:226427/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTTTATTGCCGGGATGGGGGGAtt > 2:509086/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGGTTGCTATTTCCGGGATGGGGGGAtt > 1:55243/1‑90 (MQ=255) gAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTAc < 1:429504/90‑1 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTTCCGGGATGGGGGGATTACTgg > 1:457531/1‑90 (MQ=255) gCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTGCCGGGATGGGGGGATTACTtgcg > 1:383404/1‑90 (MQ=255) ttCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCa < 1:226427/90‑1 (MQ=255) gAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTGCCGGGATGGGGGGATTACTGGCGGCGTCTa > 2:89244/1‑90 (MQ=255) ttGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTcgcg > 1:146933/1‑90 (MQ=255) tGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTCTATTGCCGGGATGGGGGGATTACTGGCGGCCTCTATTTTCGCGCCgct > 2:92346/1‑88 (MQ=255) cTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACg < 2:146933/90‑1 (MQ=255) gTTTCCGCAATATCACCTTGAGGCGGGGACGTTTTCTTTTGCCGGGATGGGGGGGGTACTGGGGGCCTCTATTCGCGCCGCGTTAACggg > 1:344688/1‑90 (MQ=255) gTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACggg > 1:152521/1‑90 (MQ=255) ttGAGGCGGGGACGTTTGCTATTGCCGGGATGGGGGCATTACTGGCGGGATCTATTCGCGCGCCGTTAAAGGGGGGCAttctggttctgg > 2:300614/1‑90 (MQ=255) gAGGCGGGGGCGTTTGCTATTTCCGGGATGGGGGGATTACTGGCGGGGTCTATTTCCGCGCCGGGAACGGGGGGCAttctggttctgggg > 2:453501/1‑88 (MQ=255) gAGGCGGGGACGGTTGCTATTTCCGGGATGGGGGGATTACTGGCGGCCTCTATTCTCGCGCCGTTAACGGGGGTCAttttgggtttgggg > 2:175411/1‑88 (MQ=255) cGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAcc < 2:224354/86‑1 (MQ=255) cGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGAcc > 1:224354/1‑86 (MQ=255) cGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATa < 2:383404/90‑1 (MQ=255) cGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATa < 2:55243/90‑1 (MQ=255) | ACTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATA > NZ_CP009273/172686‑172845 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |