| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,431,979 | 0 | A | T | 42.9% | 18.4 / 10.5 | 21 | D5E (GAT→GAA) | msrA | peptide‑methionine (S)‑S‑oxide reductase MsrA |
| Reads supporting (aligned to +/- strand): ref base A (2/10); new base T (9/0); total (11/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCG > NZ_CP009273/4431896‑4432068 | ccGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAAt < 1:392949/90‑1 (MQ=255) cGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATa < 1:182038/90‑1 (MQ=255) gTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTCaaaaaa > 1:146572/1‑90 (MQ=255) gTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTTaaaaaa > 1:432979/1‑90 (MQ=255) tGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGATGCTTTTTTTCAaaaaaaa > 1:271119/1‑89 (MQ=255) tGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAAc > 1:140777/1‑90 (MQ=255) gCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGCTTTTTTTTAAATAAAATCATg > 1:296680/1‑90 (MQ=255) gCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGGAACCAGATGGTTTTTTTTAAATAAAATCATGGTGTCGc > 2:302901/1‑90 (MQ=255) gTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGc‑tc < 1:405589/90‑1 (MQ=255) aTCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCCTTTTTTCAAATAAAATCATGGGGTCGCGCCTCCCAaaaacg > 2:350649/1‑88 (MQ=255) tgtTACGTCCAGGCAGGGCATCGGCGGGGGGAAACAGAAGCTTTTTTTTAAAAAAAATCATGGGGTTGC‑TCTCCCCAAAAACGGGCAttt > 2:243913/1‑89 (MQ=255) tgtTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTg < 1:379690/90‑1 (MQ=255) tCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTaa < 1:293273/90‑1 (MQ=255) gcaggGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTg < 2:427774/90‑1 (MQ=255) gggCATCGGCGGGGGGAACCAGATGCTTTTTTTCAAATAAACTCATGGGGTCGC‑TCTCCCGAAAATCGGTCATTGGGGGTAAGGGTGTaa > 1:306790/1‑90 (MQ=255) gggCATCGGCGGGGGAAACCAGATGCTTTTTTTCAAATAAAATCATGGGGTCGC‑TCTCCCGAAAAACGGTCATTGGGGGTAAGGTTGTaa > 2:159422/1‑90 (MQ=255) gggCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGGGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGGTGTaa > 1:136089/1‑90 (MQ=255) aCCAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGca < 2:433459/90‑1 (MQ=255) cAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGcaca < 1:52422/90‑1 (MQ=255) tttttATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAac < 2:146572/90‑1 (MQ=255) tttttATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAac < 2:296680/90‑1 (MQ=255) aTCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCg < 2:271119/90‑1 (MQ=255) | CCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCATGGTGTCGC‑TCTCCCGAAAATCGGTCATTGGGGTTAAGGTTGTAACAAGAGGCGTATTTGCACACAATAACCACCG > NZ_CP009273/4431896‑4432068 |
| 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: ‑ |