| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 886,740 | 0 | C | G | 61.9% | 15.2 / 26.0 | 21 | A34G (GCC→GGC) | nfsA | nitroreductase NfsA |
| Reads supporting (aligned to +/- strand): ref base C (7/1); new base G (5/8); total (12/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.69e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAA > NZ_CP009273/886659‑886816 | ttATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGa < 1:56766/90‑1 (MQ=255)ttATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGa < 2:812344/90‑1 (MQ=255) tATTTGTGGCCATCGCTCCATTCGCCATTTCAGTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGAc > 1:24035/1‑90 (MQ=255) tATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGAc > 2:865584/1‑90 (MQ=255) ggCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGtccagt > 1:721282/1‑90 (MQ=255) cATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGtccagttcc > 2:188323/1‑90 (MQ=255) aTTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTtg > 1:364812/1‑90 (MQ=255) gCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTtgcagt < 2:92375/90‑1 (MQ=255) gCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTtgcagt < 2:364812/90‑1 (MQ=255) gCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTtgcagt < 1:4853/90‑1 (MQ=255) ccATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTtgcagtg > 1:59449/1‑90 (MQ=255) cATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTtgcagtgc < 1:68487/90‑1 (MQ=255) aaCCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTc < 2:110139/90‑1 (MQ=255) aCCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCg > 1:82097/1‑90 (MQ=255) tttCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGc > 1:729889/1‑77 (MQ=255) tttCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGc < 2:729889/77‑1 (MQ=255) ccGAAGCGCAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCg > 1:722837/1‑90 (MQ=255) gcgcAGCGTGAGGCGATTATTAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACaaa > 2:574522/1‑90 (MQ=255) aGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTg > 1:859720/1‑90 (MQ=255) attattAACAGCGGCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTgaagaa < 2:481640/90‑1 (MQ=255) attattAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTgaagaa > 1:47053/1‑90 (MQ=255) | TTATTTGTGGCCATCGCTCCATTCGCCATTTCACTGATGAACCCATTTCCGAAGCGCAGCGTGAGGCGATTATTAACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAA > NZ_CP009273/886659‑886816 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |