Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,431,963 | 0 | A | G | 52.4% | 11.3 / 28.0 | 21 | S11P (TCC→CCC) | msrA | peptide‑methionine (S)‑S‑oxide reductase MsrA |
Reads supporting (aligned to +/- strand): ref base A (2/8); new base G (11/0); total (13/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGA > NZ_CP009273/4431879‑4432042 | cATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCa > 2:52074/1‑90 (MQ=255) ttGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGa > 2:335431/1‑90 (MQ=255) tGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGAt < 2:219252/90‑1 (MQ=255) ggTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGGTg > 2:488036/1‑90 (MQ=255) cATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAAAAGCtt > 1:299679/1‑90 (MQ=255) tGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACCGATGCgtttt > 2:189354/1‑90 (MQ=255) tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGGAACCGATGGtttttttta > 1:564329/1‑86 (MQ=255) gACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACCGATGCTTTTTTTCaa > 1:381788/1‑90 (MQ=255) cGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTCAAAAAAACTCAt > 2:635282/1‑90 (MQ=255) gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAAAAAAATCATg > 2:686519/1‑90 (MQ=255) gCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGc < 2:418475/90‑1 (MQ=255) gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGctc < 1:680577/90‑1 (MQ=255) gggCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGgtgt < 1:240377/77‑1 (MQ=255) gggCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAAACAGATGCTTTTTATCAAATAAACTCATGgtgt > 2:240377/1‑77 (MQ=255) aTCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAAATGCTTTTTTTTAAATAAAATCATGGTGTCGCTCTCCCGGAAAACg > 2:519251/1‑90 (MQ=255) tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg < 1:335431/90‑1 (MQ=255) gTCCAGGCAGGGCATCGGCGGGGGGAAACAAATGGTTTTTTTTAAAAAAAATCATGGGGTCGGTCTCCCCAAAAACGGGCGTTGGGGGTa > 1:688445/1‑90 (MQ=255) tCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTaa < 2:456519/90‑1 (MQ=255) ggcaggGCATCGGCGGGGGGAAACCGATGCCTTTTTTCTAAAAAAAACATGGGGTCGGGCTCCCCAAAAACGGTCCTTGGGGGTAAGGtt > 1:456519/1‑90 (MQ=255) gCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACa < 1:632343/90‑1 (MQ=255) tCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAaga < 2:688445/90‑1 (MQ=255) | CATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGA > NZ_CP009273/4431879‑4432042 |
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: ‑ |