Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 792,056 | 0 | C | T | 53.8% | 13.2 / 28.9 | 27 | T16I (ACT→ATT) | modC | molybdenum ABC transporter ATP‑binding protein ModC |
Reads supporting (aligned to +/- strand): ref base C (7/5); new base T (0/14); total (7/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.20e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGT > NZ_CP009273/791974‑792144 | ggcCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAAt > 2:2274709/1‑90 (MQ=255) gAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAAc < 1:1112824/90‑1 (MQ=255) gAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAAc < 1:1165055/90‑1 (MQ=255) cAGCCGTGAACGGGGGGGGCGCTAATCATGCGGGAACTGATTTTTCCCCAGACGTTGGGCAACCATTGCCTGATTATTAATGAAACGCTg < 1:2302751/90‑1 (MQ=255) gagAAGGGGGGGGGCGCTAATCATGTTGGAACTAAATTTTCCCCAGGCGTTGGGCACCCATTGCCTGATTATTAATGAAACGCTGCCCGc < 1:794606/88‑1 (MQ=255) gAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCa > 2:1060175/1‑90 (MQ=255) gggggcgCTAATCATGCTGGAACTAATTTTTCCCCAGAGGTTGGGCACCCTTTGCCTGATTATTAATAAACCGCTGCCCGCCAATGGCAt < 1:1923759/89‑1 (MQ=255) gggggccCTAATCATGCTGAAATTTATTTTTCCCCAGCCGTTGGCCACCCTTTCCCTGATTATTAATAAACCGCTGCCCGCCAATGGCAt < 1:83234/89‑1 (MQ=255) ggggccgCTATTCATGCTGGAACAAATTTTTCCCCAGAGGTGGGACAACCATTGCCTGATTATTAATAAACCGCTGCCCGCCAATGGCAt < 1:2017065/89‑1 (MQ=255) ggggggCTAATCGGGCGGGAACTGATTTTTCCCGGGACGTTGGCCAACCATTGCCTGAGTATTAATAAACCGCTGCCCGCCAATGGCATc < 2:129723/90‑1 (MQ=255) ggggcgCTAACATTGCTGGAACTAATTTTTCCCCAGACGTTGGGCAACATTTGCCTCCCTATAAATAAACCGCTCCCCGCCAATGGAATc < 1:544239/90‑1 (MQ=255) ggcgTTAATCATGCTGAAACTAATTTTTTCCCAGACGTGGGGCAACCTTTGCCTGATTATTAATAAAACGCTGCCCGCCAATGGCATCAc < 2:1219210/90‑1 (MQ=255) gAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCgg > 1:527073/1‑90 (MQ=255) gACGTGGGGCAACCATTGCCTGATTTTTAATAAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAAc < 2:983010/90‑1 (MQ=255) tGGACACCCATTGCTTGTTTATTAACGAAACCCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGc < 1:2002159/90‑1 (MQ=255) aaCCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGAtt > 2:453190/1‑90 (MQ=255) ccctttGCTTGATTTTAAAAGAACCGCCCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTa < 1:1484922/86‑1 (MQ=255) aCCTTTGCCTGATTATAAAGGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTa < 1:770183/90‑1 (MQ=255) aCCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTa < 1:1246099/90‑1 (MQ=255) aTTGCCTGATTATAAATAAAACGCGGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACg < 1:1622071/90‑1 (MQ=255) tGCCTGTTTTTTAATAAACCCCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGcc < 2:1326130/90‑1 (MQ=255) tGCCTGATTATTAATGAACCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGcc < 2:2458380/90‑1 (MQ=255) tGCCTGATTATAAAAGAACCGCTCCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGcc < 1:1691059/90‑1 (MQ=255) tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTa > 1:1081982/1‑85 (MQ=255) tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTa < 2:1081982/85‑1 (MQ=255) gACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAg > 2:1631956/1‑90 (MQ=255) aCTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGt > 1:512072/1‑90 (MQ=255) | GGCCAGAATCAGCCGTGAACGGGCGGGGCGCTAATCATGCTGGAACTGAATTTTTCCCAGACGTTGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGT > NZ_CP009273/791974‑792144 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |