Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 172,772 | 0 | A | G | 63.3% | ‑2.8 / 17.4 | 30 | G393G (GGA→GGG) | clcA | H(+)/Cl(‑) exchange transporter ClcA |
Reads supporting (aligned to +/- strand): ref base A (6/5); new base G (19/0); total (25/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.24e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATA > NZ_CP009273/172687‑172845 | cTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGCCGGGATgg > 1:456126/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGGGCTGTTTCCGCAATATCACCTTGAGGCGGGGGGGGTTGCTATTTTCGGGAGgggg > 2:1221552/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATgggg > 2:423689/1‑90 (MQ=255) gtgCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTTCCGGGATgggg > 1:1036297/1‑90 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGGTATTGCCGGGGTgggggga > 2:1010677/1‑88 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATgggggga > 1:536341/1‑88 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTCTTTTGCCGGGATgggggga > 1:1128294/1‑88 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCa < 1:229637/90‑1 (MQ=255) cTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTCTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGAATgggggga > 2:486822/1‑88 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTGCCGGGATgggggggtt > 1:630/1‑86 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTCTATTGCCGGGATgggggggtt > 1:810553/1‑86 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGGAtt > 2:734609/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTCTTTTGCCGGGATgggggggtt > 1:599377/1‑86 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTTCTATTTCCGGGATGGGGGGAtt > 1:959083/1‑90 (MQ=255) gggAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGGATGGGGGGAtt > 2:419427/1‑90 (MQ=255) aaCCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGGATGGGGGGATTACt > 2:28927/1‑90 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTTTTTTTTGCGGGATGGGGGGATTTCttg > 2:904057/1‑88 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGTCGGGATGGGGGGATTACttg > 1:664566/1‑88 (MQ=255) ccGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGGTTGCTATTGCCGGAATGGGGGGATTACTgg > 1:190356/1‑90 (MQ=255) gCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTGCCGGGATGGGGGCATTACTggcg > 2:81907/1‑90 (MQ=255) cTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTggcgg < 1:751366/90‑1 (MQ=255) ttGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTTCTTTTTCCGGAATGGGGGGGTTTCTGGCGGGGTCTATTctcg > 1:764278/1‑90 (MQ=255) tGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTcgcgc < 2:777422/90‑1 (MQ=255) cGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGGCGTTTGCTATTTCCGGGATGGGGGCATTACTGGCGGCATCTATTCGCGCGcc > 2:551213/1‑90 (MQ=255) ttCCGCAATATCACCTTGAGGCGGGGACGTTTTCTATTGCCGGGATGGGGGCATTACTGGCGGCATCTTTTTGCGCGCCGGTAAAggggg > 2:591914/1‑89 (MQ=255) atatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttct < 2:871296/90‑1 (MQ=255) tatCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCAttctg > 2:351399/1‑90 (MQ=255) ggggGCGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTTGCGCGCCGTTAACGGGGGTCATTCTGGGTCTTGGGATGAc > 1:311302/1‑90 (MQ=255) ggggACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGGTAACGGGGGTCATTCTTGTTTTGGGGATGAc > 1:512180/1‑90 (MQ=255) cGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATa < 2:952743/90‑1 (MQ=255) | CTGTGCTGGGAACCGCTTTCGGAATGGTTGCCGTTGAGCTGTTTCCGCAATATCACCTTGAGGCGGGGACGTTTGCTATTGCCGGAATGGGGGCATTACTGGCGGCATCTATTCGCGCGCCGTTAACGGGGATCATTCTGGTTCTGGAGATGACCGATA > NZ_CP009273/172687‑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: ‑ |