Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,091,435 | 0 | T | G | 57.9% | 14.3 / 19.7 | 19 | intergenic (+77/‑183) | sodA/kdgT | superoxide dismutase [Mn]/2‑keto‑3‑deoxygluconate transporter |
Reads supporting (aligned to +/- strand): ref base T (1/7); new base G (11/0); total (12/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.59e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAT > NZ_CP009273/4091382‑4091521 | ggCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAAATTTAAAAAAAAAAAAAAGaac > 2:365367/1‑88 (MQ=255) gCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATcc > 2:16493/1‑90 (MQ=255) gTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGa > 2:229414/1‑90 (MQ=255) ataGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGata < 1:212882/90‑1 (MQ=255) ggCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTGTGATCAATTTTAAAATAAAAAAAAAGATCCGataaaa > 1:458224/1‑90 (MQ=255) ggCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAATGATCCGgtaaaa > 2:532534/1‑90 (MQ=255) ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGTTTTTTTGGGATGAATTTTAAAATAAAAAAAATGATCCGataaaaaa > 1:529329/1‑89 (MQ=255) ccGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGcaaaaaaa > 2:170998/1‑89 (MQ=255) cGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGGGGGTTTTTTGTGATCAATTTTAAAAAAAAAAAAATGATCCGataaaaaaa > 1:358418/1‑88 (MQ=255) gCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAATAAAAAAAAAGATCCGataaaaataa > 1:288270/1‑90 (MQ=255) cATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGataaaaaaaaa > 2:32046/1‑90 (MQ=255) tatCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAAc < 1:243533/90‑1 (MQ=255) atCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACa < 2:203281/90‑1 (MQ=255) gcttaaaaaTGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGtt < 2:411354/86‑1 (MQ=255) aaaaaTGAACCATCGCCAACGGCGGTGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAATAAAAAAACGTTTCaa > 2:243533/1‑90 (MQ=255) aCCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGAtgtg < 2:169784/90‑1 (MQ=255) tCGCCAACGGCGGGGGGTTTTTTTTGATCAATTTTAAAAAAAAAAAAAAGATCCGATAAAAAAAAAAAAACGTTTTAATTGATGTGGttt > 2:524597/1‑90 (MQ=255) ccAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGa < 2:307660/90‑1 (MQ=255) aaCGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATc < 1:92895/90‑1 (MQ=255) tggtTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAt < 2:413455/90‑1 (MQ=255) | GGCGTATAGGCCGCATATCAGCTTAAAAAATGAACCATCGCCAACGGCGGTGGTTTTTTTGTGATCAATTTCAAAATAAAAACAATGATCCGATAAAAATAAAACAGCGTTTCAATTGATGTGGTTTTGATCACTTTTAT > NZ_CP009273/4091382‑4091521 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |