| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,571,661 | 0 | A | C | 52.4% | 16.4 / 27.6 | 21 | G142G (GGT→GGG) | gntR | gluconate operon transcriptional repressor GntR |
| Reads supporting (aligned to +/- strand): ref base A (10/0); new base C (0/11); total (10/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.84e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGA > NZ_CP009273/3571576‑3571740 | tCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCg > 1:84653/1‑90 (MQ=255) cGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACt > 2:519445/1‑90 (MQ=255) ctgctTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCaa > 2:459035/1‑90 (MQ=255) gctTAAAATTTACAAAACCCGACGGCGATAACAGGGCAGGGCGTCTTG‑CTGTCCATCATTCCCCCCAGGGAAATCCCCGCCACTTCAatc < 1:250026/90‑1 (MQ=255) tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg > 1:273697/1‑90 (MQ=255) tATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAgg > 1:232812/1‑90 (MQ=255) aacccGGCGGCAATACAAAGGCAGGGCGTTTTG‑CTGTCCATCAGTCCCCCCGGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTgcgc < 2:232812/87‑1 (MQ=255) gACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACAACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCgcggctg > 2:517904/1‑88 (MQ=255) ccGGCAATATCAAGGCAGGGCGTCTTG‑CTGCCCATCAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCgtg < 2:70075/89‑1 (MQ=255) atCAGGGCATGGCGCCTTC‑CTGCCCAACAGTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACg < 2:330120/90‑1 (MQ=255) gaGGCATGGTGACTTGCCGGTCCTTCACTTCCCCCGCGGCATTCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGtt < 2:660706/89‑1 (MQ=255) aTGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGa > 1:590476/1‑90 (MQ=255) tGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAg > 1:309065/1‑90 (MQ=255) ccTTG‑CTGCCCATCATCCCCCCCCGGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa < 1:519445/89‑1 (MQ=255) aCTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCa > 1:35158/1‑90 (MQ=255) gaccAGTCCCCCCACGGAAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat < 2:319681/87‑1 (MQ=255) catcatTCCCCCCAGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat < 2:18261/90‑1 (MQ=255) catcatTCCCACCCGGGGAACCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat < 2:482225/90‑1 (MQ=255) cATCAGTCCCCCCCGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGat < 2:332577/90‑1 (MQ=255) ccccagGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGga < 2:273697/88‑1 (MQ=255) ccacGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGgagc > 1:319681/1‑89 (MQ=255) | TCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCATGGCGACTTG‑CTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGA > NZ_CP009273/3571576‑3571740 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |