Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 792,067 | 0 | A | C | 62.5% | 14.6 / 16.8 | 16 | T20P (ACG→CCG) | modC | molybdenum ABC transporter ATP‑binding protein ModC |
Reads supporting (aligned to +/- strand): ref base A (4/2); new base C (0/10); total (4/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.24e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCG > NZ_CP009273/791979‑792155 | gAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAAc < 2:341714/90‑1 (MQ=255) gAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAAc < 2:234439/90‑1 (MQ=255) ggggctCTAATCATGCTGAAACTAATTTTTCCCCAGAGGTGGGCCACCCTTTCCCTGACTATTAATAAACCGCTGCCCGCCAATGGCATc < 1:853019/90‑1 (MQ=255) cATGCTGGAACTGATTTTTCCCCAGACGTTGGGCAACCATTGCCGGACTATTAATGAACCGCTGCCCGCCAATGGCATCACTGCTATCtt < 2:135207/90‑1 (MQ=255) tgaaaCTAATTTTTCCCCAGCCGTGGGGCACCCTTTCCCTGATTATTAATGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCg < 1:776875/87‑1 (MQ=255) gAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGtc > 1:497915/1‑90 (MQ=255) gAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCgg > 2:610989/1‑90 (MQ=255) cAGACTTTGGGCAACCATTGCCTGACTATTAATGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGaaaa < 2:779309/90‑1 (MQ=255) tGGGCACCCTTTGCCTGACTATTAATAAACCGCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGc < 1:166259/90‑1 (MQ=255) tGGGCAACCATTGCCTGACTATTAAAGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGc < 2:203981/90‑1 (MQ=255) tGGGCAACCATTGCCTGACTATTAAAGAACCGCCCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGc < 1:415667/90‑1 (MQ=255) aaCCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGAtt > 1:840407/1‑90 (MQ=255) ctaaCTATAAATAAACCCCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATc < 2:840407/87‑1 (MQ=255) gattataAATGAACCGCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAg < 2:762143/87‑1 (MQ=255) aattATTAATGAACCCCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAg < 1:610989/87‑1 (MQ=255) aTTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGa > 1:240733/1‑90 (MQ=255) aaaCGCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGAcgc < 2:887955/90‑1 (MQ=255) aaCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGAcgcg > 1:882937/1‑90 (MQ=255) aaCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGAcgcg > 1:392149/1‑90 (MQ=255) aaCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGAcgcg > 1:24405/1‑90 (MQ=255) aaCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGAcgcg > 1:166164/1‑90 (MQ=255) | GAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCG > NZ_CP009273/791979‑792155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |