Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,852,466 | 0 | T | G | 57.1% | 5.1 / 20.5 | 21 | V672G (GTG→GGG) | mutS | DNA mismatch repair protein MutS |
Reads supporting (aligned to +/- strand): ref base T (7/2); new base G (0/12); total (7/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.10e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT > NZ_CP009273/2852378‑2852555 | caAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggtg < 2:245890/90‑1 (MQ=255) gTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATg < 2:297833/90‑1 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa > 1:274170/1‑90 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa > 1:335221/1‑90 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa > 2:169604/1‑90 (MQ=255) gCATTTTTCCCCGGGGGGGCGCGGCAAATGCCCTGGCGCCCGGGCGCTCACCTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTAc < 2:145501/90‑1 (MQ=255) aTCTTTCCCCGCGTGGGCGCGGCAGATGCCCGGGCGCCCGGGCCCTCACCCTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACAt < 2:244724/90‑1 (MQ=255) tccccGCGTAGGCGCGGCAGATGCCCTGGCGTCCGGGCGCCCACCCTTTATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGc < 1:701642/88‑1 (MQ=255) gcgTCGGCCCGGCAGATGCCCGGGCGGCCGGGGGCCCACCCTTTGGGGGGGAGATGACTCAAACCGCCAATATTTTACATAACGCCACCg < 1:414762/90‑1 (MQ=255) ggCGCGGCAGATGCCTGGGCGGCGGGGCGCCCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 2:620662/90‑1 (MQ=255) ggCGCGGCAGATGCCCGGGCGCCGGGGCCCCAACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATAc < 2:335221/90‑1 (MQ=255) gggaaGATGCCCTGGCGCCGGGGCCCCCAACCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 2:274170/86‑1 (MQ=255) cgGCAGATGCCCGGGCGCCGGGGCCCTAACCCTTTAGGGGGGAGATGATTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTc < 2:645426/90‑1 (MQ=255) gATGCCCTGGCGTCCGGGCGCTCACCTTTTATGGGGGAGATGATTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGgtg < 1:385775/90‑1 (MQ=255) cccgggCGCCGGGGCCCTCACCTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa < 1:619416/86‑1 (MQ=255) cccTGGCGTCCGGGCGCTCAACTTTTATGGGGGAGGTGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa < 1:702518/89‑1 (MQ=255) tGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATgg > 2:391192/1‑90 (MQ=255) cgcccGGGCGCTCAACCTTTATGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATg < 1:343993/87‑1 (MQ=255) ccGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATgaga > 1:154434/1‑90 (MQ=255) gcTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGc > 1:722406/1‑90 (MQ=255) gcTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGc > 1:247387/1‑90 (MQ=255) ttAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt < 2:167437/90‑1 (MQ=255) ggtggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc > 2:434820/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:221524/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:239040/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 2:366451/1‑90 (MQ=255) tggAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt > 1:121457/1‑90 (MQ=255) | CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCT > NZ_CP009273/2852378‑2852555 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |