| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,089,571 | 0 | A | G | 55.6% | 11.3 / 29.0 | 27 | I83V (ATT→GTT) | hisF | imidazole glycerol phosphate synthase subunit HisF |
| Reads supporting (aligned to +/- strand): ref base A (2/10); new base G (15/0); total (17/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.82e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCA > NZ_CP009273/2089495‑2089648 | tGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGa > 2:343820/1‑90 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGaag > 1:840/1‑88 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑Gaa < 1:742543/90‑1 (MQ=255) aGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTTTCTGGAAGATGCCG‑C‑Gaa > 2:592356/1‑90 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATGTCTGGGAGAGGGCG‑C‑GGAAAttt > 1:546146/1‑89 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTGAAGAGTCTGGGAAATGGCGCC‑GAAAAtt > 1:508558/1‑90 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAATTTTTGGAAGAAGGCG‑CGGAAAAtt > 1:974309/1‑90 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAg > 1:357894/1‑59 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGGTATAGCTGTGGGAGAAGCCG‑C‑GGAAATTTTTTTCt > 1:477415/1‑90 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAg < 2:357894/59‑1 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGTTCTGGAAGAAGACG‑C‑GAAAATTTTTTTCt > 1:512682/1‑90 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGGCTCTTGGAGAAGCCC‑C‑GCAAAATTTTTCCt > 1:366874/1‑90 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTAAAACTCTGGGAGAAGCCC‑C‑GAAAATTCTTTTCt > 2:608182/1‑90 (MQ=255) tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGgc < 2:273301/90‑1 (MQ=255) ctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCC‑TTGgcgc < 2:242402/90‑1 (MQ=255) tcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCG‑C‑GAAAAATCTTTTCTTTGgcgc > 2:393367/1‑90 (MQ=255) tcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGgcgc < 2:732225/90‑1 (MQ=255) aaGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAAttt > 2:359355/1‑90 (MQ=255) aTCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATc < 1:608182/90‑1 (MQ=255) gccATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAAc < 2:622865/88‑1 (MQ=255) gACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAAc < 1:592356/90‑1 (MQ=255) ttCCGTTTTGTGTGGCGGGTGGGGTTAAGTCTCTGGAAGATGCCC‑C‑GAAAAATCTTTTCTTTTGCGGGGGGAAAATTTTCATCAActccc > 1:152934/1‑88 (MQ=255) ttCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAGATGCCG‑C‑GAAAAATTTTTTCTTTGGCGCGGGTAAAATTTTCATCAAAtctc > 1:373560/1‑90 (MQ=255) gTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTTCTTTGGCGCGGGTAAAAATTTCATCAACTCTCCTgc > 2:17304/1‑90 (MQ=255) ttGTGTGGCGGGGGGGGTTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTTTTTTCTTTTGCGCGGGTAAAAATTTCATCAACTCTCCTGCGCt > 1:315183/1‑90 (MQ=255) gtgGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCa < 1:811277/90‑1 (MQ=255) gtgGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCa < 1:230486/90‑1 (MQ=255) | TGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAGATGCCG‑C‑GAAAATTCTTTCCTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCA > NZ_CP009273/2089495‑2089648 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |