| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,089,567 | 0 | T | G | 55.6% | 22.5 / 29.4 | 27 | G81G (GGT→GGG) | hisF | imidazole glycerol phosphate synthase subunit HisF |
| Reads supporting (aligned to +/- strand): ref base T (0/12); new base G (15/0); total (15/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.75e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCA > NZ_CP009273/2089478‑2089654 | atCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGt < 2:840542/90‑1 (MQ=255) gCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGggggggggt > 2:825288/1‑87 (MQ=255) cTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTa < 2:1497717/90‑1 (MQ=255) tCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAg < 2:145721/90‑1 (MQ=255) cAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGtc < 2:1474453/90‑1 (MQ=255) aTGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGt > 1:962273/1‑74 (MQ=255) aTGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGt < 2:962273/74‑1 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAA‑ATg > 1:1326999/1‑90 (MQ=255) gtgtGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATTTCTTGAAA‑ATg > 1:1117264/1‑90 (MQ=255) gtgGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGGAG‑ATG‑cc > 1:597501/1‑90 (MQ=255) aaaGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGAAG‑AAG‑CCGC‑Gaaaaatt > 1:1074263/1‑87 (MQ=255) aGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCtt < 2:298359/90‑1 (MQ=255) aGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCtt < 2:1493493/90‑1 (MQ=255) gCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTAATCTCTTGGAG‑ATG‑GCGC‑GGAAAtttttt > 1:482480/1‑90 (MQ=255) tGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTTAATCTTTGGGAG‑ATG‑CCGCGGAAAATTTTTTc > 1:247647/1‑90 (MQ=255) tGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGAAG‑ATG‑CCCC‑GAAAAATtttttc > 1:1178013/1‑88 (MQ=255) gggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGGAG‑ATG‑CCGC‑GGAAAATTTTTC‑Ct > 2:1102819/1‑90 (MQ=255) ggTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGGTTAGGTTCTTGGAG‑ATG‑GCGC‑GGAAATTTTTTT‑Ctt > 2:654248/1‑90 (MQ=255) tctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAAATCTTTC‑CTTTTgc > 1:1213923/1‑90 (MQ=255) ctcGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGgcg < 1:764500/90‑1 (MQ=255) cgTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGAGTAAGTCTCTGGGAG‑AAG‑CCGC‑GGAAAATTTTTT‑CTTTTGCGCggg > 1:1300157/1‑89 (MQ=255) ggAAGTGATCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAATCTCTGGAAG‑AAG‑CCCC‑GAAAAATTTTTT‑CTTTTGGGGGGGTAAAAt > 2:1248624/1‑90 (MQ=255) tCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCa < 1:1446920/90‑1 (MQ=255) tCGACATTCCGTTTTGTGTGGCGGGGGGGGTTAAGATTTTGGGAGAAGG‑CCCG‑GAAAATTTTTTCTTTTTGGCGGGGA‑AAAATTTcctca > 2:747406/1‑87 (MQ=255) ttttGTGTGGCGGGGGGGGTTTAATCTCTTGAAG‑AAGACCCC‑GAAAA‑TCTTTT‑CTTTTGTGCGGGTAAAAATTCCACCAATTCCTCTgcc > 2:1540008/1‑89 (MQ=255) tttGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTgcgc < 2:1405947/90‑1 (MQ=255) gtgtgGCGGGGGGGGTTAAGTCTCTTGGAG‑ATG‑CCGC‑GAAAATTTTTTT‑CTTTTGGGGGGGTAAAATTTTCATCAACTCTCCTGCGCTgg > 2:462405/1‑90 (MQ=255) gggTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAccc < 1:1048117/90‑1 (MQ=255) gggTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAccc < 1:825288/90‑1 (MQ=255) gggTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGAccc < 1:1102819/90‑1 (MQ=255) ggTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCa < 1:1234479/90‑1 (MQ=255) | ATCACCGCTTCCAGCGATGGCCGTGTGGTAGATAAAAGCTGGGTATCTCGCGTGGCGGAAGTGATCGACATTCCGTTTTGTGTGGCGGGTGGGATTAAGTCTCTGGAAG‑ATG‑CCGC‑GAAAATTCTTTC‑CTTTGGCGCGGATAAAATTTCCATCAACTCTCCTGCGCTGGCAGACCCA > NZ_CP009273/2089478‑2089654 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |