| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,852,466 | 0 | T | G | 48.5% | 11.2 / 11.7 | 33 | V672G (GTG→GGG) | mutS | DNA mismatch repair protein MutS |
| Reads supporting (aligned to +/- strand): ref base T (14/3); new base G (1/15); total (15/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.26e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT > NZ_CP009273/2852384‑2852555 | gTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATg > 1:330862/1‑90 (MQ=255) ttGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGGTGACCTGGCGTCCGGGCGATCAACCTTTATGGTGGAGATGACTGaaa < 1:1500631/90‑1 (MQ=255) tGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAAc > 2:40918/1‑90 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGGGGAGATGACTGAAACCGCCaa > 1:1339072/1‑90 (MQ=255) tcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAata > 1:224937/1‑90 (MQ=255) gatcgCTTCTTTCCCCGGGGAGGCGCGGCAGATGCCCGGGCGTCGGGGCGCTAACCCTTTAGGGGGGAGATGACTGAAACCGCCAATAtt < 2:1083128/90‑1 (MQ=255) tcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATAtttt > 2:1191206/1‑90 (MQ=255) tttCCCCGCGGGGGCGGGGCAGATGCCCGGGCGTCGGGGCCCTCACCTTTTAGGGGGGAGATGACTGAACCCGCCAATATTTTACATAAc < 1:460225/90‑1 (MQ=255) ttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACg > 2:1497105/1‑90 (MQ=255) tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGc > 2:1120106/1‑90 (MQ=255) tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGc > 2:1164387/1‑90 (MQ=255) cccGCGTGGGCGCGGCAGATGCCCGGGCGCCGGGGCCCTCACCTTTTAGGGTGGAGATGACAAAAACCGCCAATATTTTACATAACGCCa < 1:35055/90‑1 (MQ=255) cgcgTAGGCGCGGAAGATGCCCTGGCGCCCGGGCGCTCACCTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAcc < 2:1432947/90‑1 (MQ=255) aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATa > 2:449944/1‑90 (MQ=255) aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATa > 1:831644/1‑90 (MQ=255) ggggcgGCAGATGCCCGGGCGCCGGGGCCCTCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 1:684336/87‑1 (MQ=255) ggCGCGGCAGATGCCCGGCGGCCCGGGCCCCCACCTTTTAGGGGGGAAATGCCTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 2:1276510/90‑1 (MQ=255) gggCAGATGCCGGGGCGCCGGGGCCCTCACCCTTTAGGGGGGAGATGATTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 2:231428/89‑1 (MQ=255) cgGCAGATGCCCGGGCGCCCGGGCCCCCACCTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 2:368120/90‑1 (MQ=255) ggCAGATGCCTGGGCGTCCGGGCGCCCACCCTTTAGGGGGGAGATGACAAAACCCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 2:1339072/90‑1 (MQ=255) ggCAGATGCCCTGGCGTCGGGGCGCTCACCTTTTAGGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 2:896/90‑1 (MQ=255) gCAGATGCCCGGGCGCCCGGGCCCCCACCCTTTAGGGGGGAGATGACGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTg < 1:532040/90‑1 (MQ=255) cccgggCGTCGGGGCCCCCAACTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa < 2:478712/86‑1 (MQ=255) cccTGGCGTCGGGGCGCTCACCCTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa < 2:1195053/89‑1 (MQ=255) ggCGCCCGGGCGCCCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGa < 2:224937/90‑1 (MQ=255) cGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAATCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATg > 2:1304782/1‑90 (MQ=255) cGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATg > 1:970515/1‑90 (MQ=255) gggCGCTCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc < 2:848960/90‑1 (MQ=255) gggCGCTCACCCTTTAGGGGGGGGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc < 2:506379/90‑1 (MQ=255) cgcTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCggg > 2:433686/1‑90 (MQ=255) ccttttAGGGGGGAGATGGCTAAACCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGa < 2:212839/87‑1 (MQ=255) cTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc > 1:1251241/1‑90 (MQ=255) aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc > 2:1216300/1‑90 (MQ=255) tggtggAGATGACTGAAATCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCa < 1:1304782/90‑1 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:362111/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:60035/1‑90 (MQ=255) | GTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT > NZ_CP009273/2852384‑2852555 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |