Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,431,979 | 0 | A | T | 23.3% | 105.1 / 10.0 | 43 | D5E (GAT→GAA) | msrA | peptide‑methionine (S)‑S‑oxide reductase MsrA |
Reads supporting (aligned to +/- strand): ref base A (6/27); new base T (10/0); total (16/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.18e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCG > NZ_CP009273/4431893‑4432068 | tGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCa > 1:572657/1‑90 (MQ=255) gACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTCaa > 1:477948/1‑90 (MQ=255) ccGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAAt < 1:313575/90‑1 (MQ=255) gTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGGGGGGGAAACCAGATGCTTTTTATCAAATaa < 1:237824/90‑1 (MQ=255) gTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTaaaaaa > 2:198934/1‑90 (MQ=255) tGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGCTTTTTTTTAaaaaaaa > 2:64692/1‑89 (MQ=255) cGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCAt < 2:168840/90‑1 (MQ=255) gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGTTTTTTTTAAAAAAAATCATg > 2:87678/1‑90 (MQ=255) gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAAt < 2:227477/81‑1 (MQ=255) gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAAt > 1:227477/1‑81 (MQ=255) gCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGc < 1:370165/90‑1 (MQ=255) gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGctc < 2:486215/90‑1 (MQ=255) gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGctc < 2:387622/90‑1 (MQ=255) ggCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAAAAAAATCATGGGGTTGCtctc > 2:190435/1‑90 (MQ=255) gggCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGaaa > 1:164500/1‑90 (MQ=255) ggggTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTc < 2:163108/90‑1 (MQ=255) gtgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCAtt < 2:477948/90‑1 (MQ=255) tgtTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGGTTTTTTTTAAAAAAAATCATGGGGTTGCTCTCCCCGAAAACGGGCAttt > 2:74057/1‑89 (MQ=255) tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg < 1:64692/90‑1 (MQ=255) tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTg < 2:423273/90‑1 (MQ=255) ttACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTggg < 2:572657/90‑1 (MQ=255) tCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAAAAAACACATGGTGGCGCTCCCCCCAAAAACGGGCATTGGGGGTaa > 2:318260/1‑90 (MQ=255) caggcaggGCATCGGCGGGGGGAAACAGGTGGTTTTTTTTAAATAAAATCAAGGGGGGGCTCTCCCGAAAAACGGTCATTGGGGGTTAgg > 2:259125/1‑90 (MQ=255) caggcaggGCATCGGCGGGGGGAAACAGATGCTTTTTTTTAAAAAAACTCATGGGGGCGCTCTCCCCAAAAAAGGTCATTTGGGGTAAgg > 2:131994/1‑90 (MQ=255) caggcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGGGTCGCTCTCCCGAAAATCGGGCATTGGGGGTAAgg > 1:195796/1‑90 (MQ=255) gggCATCGGCGGGGGGAAACAGAAGGTTTTTTTTAAAAAAAATCATTGGGTCGCTCTCCCCAAAAACGGGCATTGGGGGTAAGGTTGTaa > 2:349365/1‑90 (MQ=255) gcggGGGAAACCAGATGCTTTTTATCAAATAACCTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGc < 2:77849/90‑1 (MQ=255) gcggGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGc < 1:56171/90‑1 (MQ=255) gcggGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGc < 2:164500/90‑1 (MQ=255) aaCCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGc < 1:484420/90‑1 (MQ=255) aCCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGca < 2:533030/90‑1 (MQ=255) aCCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGca < 2:579904/90‑1 (MQ=255) cAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGcaca < 1:198934/90‑1 (MQ=255) cAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCAGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGcaca < 1:87678/90‑1 (MQ=255) gATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGcacaca > 1:150057/1‑90 (MQ=255) tttttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGcac > 1:222699/1‑82 (MQ=255) tttttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGcac < 2:222699/82‑1 (MQ=255) tttttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAac < 1:259125/90‑1 (MQ=255) tttttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAac < 1:318260/90‑1 (MQ=255) tttttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAac < 2:106772/90‑1 (MQ=255) tttttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAac < 1:90651/90‑1 (MQ=255) ttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAaccac < 2:340280/90‑1 (MQ=255) ttATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAaccac < 1:253767/90‑1 (MQ=255) aTCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCg < 2:232289/90‑1 (MQ=255) aTCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCg < 1:190435/90‑1 (MQ=255) | TGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGCTCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCG > NZ_CP009273/4431893‑4432068 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |