Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,224,171 | 0 | A | C | 54.5% | 14.3 / 29.4 | 22 | P162P (CCA→CCC) | pgi | glucose‑6‑phosphate isomerase |
Reads supporting (aligned to +/- strand): ref base A (1/9); new base C (12/0); total (13/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGC > NZ_CP009273/4224083‑4224254 | cGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCAt < 2:435798/90‑1 (MQ=255) gATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATa < 1:344887/90‑1 (MQ=255) cGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGGCCTCCGCCCCCACATgggggg > 2:279130/1‑86 (MQ=255) aGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTTTTACCCCCGCCCCCACCTGGTGGCCGAAg > 1:29353/1‑90 (MQ=255) ttATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTcc < 1:122228/90‑1 (MQ=255) tataCCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTTACCTCCGCCCCCACATGGTGACCGAAGGTCTTCGTCCg > 2:311970/1‑90 (MQ=255) gCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTTTGACCCCCGCCCCCACAAGGTGACCCGAACATTGCGCCCGTGCaaaa > 1:389717/1‑90 (MQ=255) gCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCCGCCCCTACATTGGGAACGGAGCTCTGCGGCCGTACaaaa > 2:12906/1‑90 (MQ=255) gCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGGTCTGACCCCGGCCCCCACATTGGGGCCGAAGGGCTGCGTCCGCGCaaaa > 1:314955/1‑90 (MQ=255) aCGTAGTGAACATCGGGATCGGCGGTTCTTACCCCCGCCCCTACATGGGGGCCGAAGCTCTGCGTCCCTACAAAAAACACCCGAAAAtac > 2:299427/1‑88 (MQ=255) aCGTAGTGAACATCGGGATCGGCGGTTCTGACCCCCGCCCCTACATGGGGGCCGAAGCTCTGCGTCCCTACAAAAAACACCCGAAAATGc > 2:216664/1‑90 (MQ=255) aGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACtt < 1:38524/90‑1 (MQ=255) gTGAACATCGGGATCGGCGGTTCTGACCCCGGCCCCTACATGGGGACCGAAGCTCTGCGTCCGTACAAAAAACACCCGAACATGCACttt > 2:38524/1‑90 (MQ=255) aCATCGGGATCGGCGGTTCTGACCTCGGCCCCTACATGGTGACCGAAGCTCTGCGTCCGTACAAAAAACACCCGAACATGCACTtttttt > 2:33960/1‑90 (MQ=255) aCATCGGGATCGGCGGTTCTGACCCCGGCCCCTACATGGTGAACGAAGCTCTGCGTCCGGACAAAAAACACCCCAACAAGCACTtttttt > 1:163797/1‑90 (MQ=255) ggcggGTCTGACCCCCGCCCCAACATGGGGACCCAAGCTCTTCGTCCGTACAAAAAACACCTCAACATGCCATTTTTTTTTTACGTCCAt > 2:571615/1‑90 (MQ=255) gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg < 1:571615/90‑1 (MQ=255) gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg < 2:34467/90‑1 (MQ=255) gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg < 1:311970/90‑1 (MQ=255) gcggTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATg < 2:575135/90‑1 (MQ=255) cTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATc < 2:265416/90‑1 (MQ=255) cGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATcgc > 1:355360/1‑90 (MQ=255) | CGATTATTTCCGGTGAGTGGAAAGGTTATACCGGCAAAGCAATCACTGACGTAGTGAACATCGGGATCGGCGGTTCTGACCTCGGCCCATACATGGTGACCGAAGCTCTGCGTCCGTACAAAAACCACCTGAACATGCACTTTGTTTCTAACGTCGATGGGACTCACATCGC > NZ_CP009273/4224083‑4224254 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |