Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,852,466 | 0 | T | G | 60.0% | ‑1.6 / 13.5 | 20 | V672G (GTG→GGG) | mutS | DNA mismatch repair protein MutS |
Reads supporting (aligned to +/- strand): ref base T (7/1); new base G (0/12); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTC > NZ_CP009273/2852378‑2852550 | caAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAtggtg < 1:211205/90‑1 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa > 2:436917/1‑90 (MQ=255) tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCaa > 2:114214/1‑90 (MQ=255) tATCGACCGCATCTTTCCCCGCGTGGGCGCGGCAGATGCCCTGGCGCCCGGGCGCTCACCCTTTATGGGGGAGATGACTGAAACCGCCaa < 1:81024/90‑1 (MQ=255) atcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAat > 2:483085/1‑90 (MQ=255) tttCCCCGCGTGGGCGCGGCAGATGCCCGGGCGCCCGGGCGCTCAACTTTTATGGGGGAGATGACTAAAACCGCCAATATTTTACATAAc < 2:328279/90‑1 (MQ=255) tccccGCGCGGGCGCGCCAGGGCACCGGGCGCCCGGGCCCCCCACTTTTAGGGGGGAATTGACAAACACCGCCAATATTTTACATAACGc < 2:332516/88‑1 (MQ=255) cccccGCGTAGGCGCGGCAGATGCCCTGGCGTCCGGGCGCTCACCCTTTATGGGGGAGTTGACTAAAACCGCCAATATTTTACATAACGc < 1:543936/88‑1 (MQ=255) ccGCGTGGGCGCGGCAGATGACCTGGCGCCCGGGCCCTCAACCTTTAGGGTGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc < 1:293988/90‑1 (MQ=255) ccGCGGGGGCGCGGCAGATGCCTGGGCGCCCGGGCGCTCAATTTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc < 2:322143/90‑1 (MQ=255) ccGCGAGGGCGCGGCAGATGCCCGGGCGCCCGGGCCCACACCCTTTAGGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCAc < 2:280468/90‑1 (MQ=255) cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCAcc > 2:368784/1‑90 (MQ=255) gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCg > 2:122214/1‑90 (MQ=255) gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCg > 2:543936/1‑90 (MQ=255) ggCAGATGCCCGGGCGCCGGGGCGCTCAACCTTTATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt < 1:35/90‑1 (MQ=255) ggCAGATGACCTGGCGTCCGGGCGCTCATCCTTTAtggt > 2:289235/1‑39 (MQ=255) ggCAGATGACCTGGCGTCCGGGCGCTCATCCTTTAtggt < 1:289235/39‑1 (MQ=255) ggCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCt > 1:29822/1‑90 (MQ=255) gaaGATGCCCGGGCGCCGGGGCGCTCAACTTTTAGGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTg < 1:122214/88‑1 (MQ=255) gggCCCTCACCTTTTATGGGGGAGATGACGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATc < 2:325967/90‑1 (MQ=255) cttttAGGGGGGAGATGATGAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa < 2:352247/88‑1 (MQ=255) ttATGGGGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGt < 1:181425/90‑1 (MQ=255) tATGGGGGAGATGACTAAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc < 1:436917/90‑1 (MQ=255) | CAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTC > NZ_CP009273/2852378‑2852550 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |