| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,964,708 | 0 | G | C | 62.1% | 20.7 / 35.5 | 29 | A27P (GCG→CCG) | ygdR | lipoprotein YgdR |
| Reads supporting (aligned to +/- strand): ref base G (6/5); new base C (7/11); total (13/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.66e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGT > NZ_CP009273/2964620‑2964797 | ggATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGc > 2:250461/1‑90 (MQ=255) gATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGa < 2:635459/90‑1 (MQ=255) gATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGa > 2:126442/1‑90 (MQ=255) aTGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGAt > 1:82758/1‑90 (MQ=255) aaTGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGGCCAAAGATGGCCGTa < 2:243786/90‑1 (MQ=255) cTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGaa > 2:47338/1‑90 (MQ=255) cTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGaa > 1:500952/1‑90 (MQ=255) cTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGaa > 1:500910/1‑90 (MQ=255) cTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGaa > 2:33645/1‑90 (MQ=255) tGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATtg < 1:464747/90‑1 (MQ=255) tGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATtg < 2:268017/90‑1 (MQ=255) ggCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATtga < 1:310534/90‑1 (MQ=255) tttGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATtgatgat < 2:601680/90‑1 (MQ=255) tttGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATtgatgat > 2:357022/1‑90 (MQ=255) ttCCGGGTGTTCCAGTGATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCgg < 1:1959/90‑1 (MQ=255) tgtTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTg < 1:250461/90‑1 (MQ=255) gtTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGa < 1:472343/90‑1 (MQ=255) gATTACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCAc > 1:618962/1‑90 (MQ=255) ttACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACg < 1:356932/89‑1 (MQ=255) ttACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACg > 2:356932/1‑89 (MQ=255) ttACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGa < 2:302200/90‑1 (MQ=255) ttACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGa < 1:478555/90‑1 (MQ=255) ttACGTCATGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGa > 1:458645/1‑90 (MQ=255) tACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGAt > 1:278110/1‑90 (MQ=255) aCGTGATGCCGACCAAAGATGGCCGTATGATTTTGACCGAGGGAAAGCCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATc < 1:442745/90‑1 (MQ=255) tCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGc < 2:278110/90‑1 (MQ=255) aTGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGa < 2:514821/83‑1 (MQ=255) aTGCCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCaa < 1:268563/90‑1 (MQ=255) tGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAg > 2:195635/1‑90 (MQ=255) gCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGt < 1:47338/90‑1 (MQ=255) | GGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACTGGCGTTTGCTGTTTCCGGGTGTTCCAGTGATTACGTCATGGCGACCAAAGATGGCCGTATGATTTTGACCGATGGAAAACCTGAAATTGATGATGATACCGGGCTGGTGAGTTATCACGATCAGCAAGGT > NZ_CP009273/2964620‑2964797 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |