Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,129,305 | 0 | G | A | 47.1% | 37.9 / 10.5 | 34 | L16L (CTG→CTA) | flgF | flagellar basal‑body rod protein FlgF |
Reads supporting (aligned to +/- strand): ref base G (16/2); new base A (0/16); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.99e-08 | |||||||||||
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) tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAAc < 2:877998/50‑1 (MQ=255) tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAAc > 1:877998/1‑50 (MQ=255) tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGcc > 1:3020/1‑90 (MQ=255) tCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGcc > 2:630636/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) | 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: ‑ |