Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 228,191 | 0 | A | G | 60.0% | 7.5 / 20.3 | 20 | S195G (AGC→GGC) | yafD | endonuclease/exonuclease/phosphatase family protein |
Reads supporting (aligned to +/- strand): ref base A (2/6); new base G (12/0); total (14/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.22e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCG > NZ_CP009273/228104‑228276 | gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGc < 2:586707/90‑1 (MQ=255) aGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGa < 2:1806/90‑1 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGGTTTTAATGCCCGGGGGCGGAGAAgg > 2:416959/1‑90 (MQ=255) aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGGGTTTAATTCCCGGGGGCGGAGAAgg > 2:568155/1‑90 (MQ=255) gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGAATTTAATTACTGGGGGCGGAGAAAgaagaa > 2:702279/1‑90 (MQ=255) ccTTCCTATTGGCGATCAGAAAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc < 2:617430/89‑1 (MQ=255) gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGGGATTTTAATTCCTGGGGGCGGAGAAGGATGAACGCGGTTTATTGCttt > 2:27905/1‑90 (MQ=255) aTAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCgg < 1:702279/90‑1 (MQ=255) caccacAGCGGCCCGGTCATTATGGC‑GGGGGGTTTTTATTTCTGGGGGCCGAGAAGGGGGAAAGCGGTTTATTGCTTTTCGCGGGgaaag > 1:262527/1‑88 (MQ=255) caccacAGCGGCCCGGTCATTATGGC‑GGGAGATTTTAATGCCTGGGGCCCTAGAAGGAAGAAAGCGTTATATCGCTTTTTGCGGGgaaag > 2:766778/1‑88 (MQ=255) gCGGCCCGGTCATTATGGCGGGGGGGTTTAAAACCCTGGGGCC‑CAGAGGGAGGAAAGCCGTATATTGCTTTTTGGGGGGAATGTCGCTgg > 1:627571/1‑89 (MQ=255) gCGGCCCGGTCATTATGGC‑GGGAGGTTTTAATGCCTGGGGGCGGAGAAGGAGGAAAGCGTTATATCGCTTTTCGCGGGGAAAGGtgctgc > 2:450280/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGGGGTTTTAATTCCCGGGGGCGGAGAAGGGGGAACGCGGTTTATTTCTTTTCGCGGGGAAAGGCGCTGCGcc > 1:351665/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGGGGGTTTAATGCCCGGGGGCGGAGAAGGGGGAAAGCGGTTTATCGCTTTTTGCGGGGAAAGGCGCGGCGcc > 1:742450/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGGGATTTTAATTCCCTGGGCCGCCGAAGGGTGAAACCGGTATATCCCTTTTCGCGGGGGAAGTCGCCGCGcc > 1:228866/1‑90 (MQ=255) gCCCGGTCATTATGGC‑GGGAGGTTTTAATTCCTGGGGCCGGAGAAGGGGGAAAGCGGTATATCGCTTTTTGCGGGGGAGGTCGCTGCGcc > 1:602450/1‑90 (MQ=255) aTTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGt < 1:755572/90‑1 (MQ=255) aTGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt < 2:322974/90‑1 (MQ=255) aTGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGcc > 1:397067/1‑90 (MQ=255) cTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCg > 1:397816/1‑90 (MQ=255) | GTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGC‑GGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCG > NZ_CP009273/228104‑228276 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |