| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,129,313 | 0 | A | G | 46.9% | 29.8 / 11.0 | 32 | Q19R (CAG→CGG) | flgF | flagellar basal‑body rod protein FlgF |
| Reads supporting (aligned to +/- strand): ref base A (15/2); new base G (0/15); total (15/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.42e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCC > NZ_CP009273/1129233‑1129391 | cGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAAc > 1:1371998/1‑90 (MQ=255)cGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAAc > 2:1408828/1‑90 (MQ=255) gCTAATCGCGGACGGGATACCACAATGGACCACGCAATTTATCCCGGGGGGGGAGCAGCCAGCCAAACACTAAAACAACGGGGGGTAAcc < 2:77463/90‑1 (MQ=255) cTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGTCGGTAACCGCCAGTaa > 2:740245/1‑90 (MQ=255) aTAGCTAAATGGATCACCAATTTTATCCCGGGAGGGGAGCACCCAGCCAGACACTAAACCAACGGGCGGAAACCGCCAGTAATCTGGCCa < 1:846627/90‑1 (MQ=255) cAATGGATCACGCAATTTACACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCt > 1:535269/1‑90 (MQ=255) aTGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCa > 1:131283/1‑90 (MQ=255) tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGcc > 2:630636/1‑90 (MQ=255) tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGcc > 1:3020/1‑90 (MQ=255) aTTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGtttt > 2:1295822/1‑90 (MQ=255) aTTTATACCGCGATGGCAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGtttt > 2:1051673/1‑90 (MQ=255) tcccGGGAGGGGAGCACCCAGCAAGACACTAAACCAACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgc < 1:1489116/88‑1 (MQ=255) taCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgc > 2:598641/1‑90 (MQ=255) taCCGCGAGGGGAGCAGCCAGCCAGACACTAAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgc < 1:489544/90‑1 (MQ=255) ccGGGAGGGGAGCAGCCAGCCAGACCCTAAATCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc < 1:302554/90‑1 (MQ=255) ccGGGAGGGGAGCAGCCAGCCAGACACTAAATAAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc < 1:1519607/90‑1 (MQ=255) ccGCGATGGGAGCACCCAGCCAGACACTAAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc < 1:630636/90‑1 (MQ=255) ccGCGAGGGGAGCAGCCAGCCAAACACTAAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc < 1:339043/90‑1 (MQ=255) gcgATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAg > 2:371646/1‑90 (MQ=255) cgATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGt > 1:702898/1‑90 (MQ=255) ggAGCACCCAGCAAAACACTAAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAAt < 1:763124/90‑1 (MQ=255) agcagcCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGc > 1:766418/1‑90 (MQ=255) agccagccAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt > 2:110619/1‑90 (MQ=255) gccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg > 1:1248435/1‑90 (MQ=255) ccagACCCAAAATAAACGGGGGGAACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgc < 2:3020/90‑1 (MQ=255) ccagACACAAAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgc < 1:735384/90‑1 (MQ=255) ccagACACAAAACAAACGGGGGGTACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgc < 2:22895/90‑1 (MQ=255) caagACACAAAACCAACGGGGGGAACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgc < 2:535269/88‑1 (MQ=255) cagACACAGAATAAACAGGCGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgcg < 1:364578/90‑1 (MQ=255) agACACTAAATCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCgg < 1:170732/90‑1 (MQ=255) aaccacTAAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCgg < 1:1138350/87‑1 (MQ=255) acTAAACCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGcc < 1:544458/90‑1 (MQ=255) aGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc > 2:179439/1‑54 (MQ=255) aGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc < 1:179439/54‑1 (MQ=255) | CGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCC > NZ_CP009273/1129233‑1129391 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |