| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,673,544 | 0 | G | A | 53.8% | 18.6 / 32.7 | 26 | K287K (AAG→AAA) ‡ | ydgH | DUF1471 family protein YdgH |
| Reads supporting (aligned to +/- strand): ref base G (6/6); new base A (0/14); total (6/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.01e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAT > NZ_CP009273/1673457‑1673632 | gatgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCggagaa < 2:888079/90‑1 (MQ=255) atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa < 2:534130/90‑1 (MQ=255) ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt < 1:795734/90‑1 (MQ=255) gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa < 1:254991/90‑1 (MQ=255) gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa < 1:269367/90‑1 (MQ=255) gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa < 2:767712/90‑1 (MQ=255) gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa < 2:332697/90‑1 (MQ=255) aTCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACcac > 2:513374/1‑90 (MQ=255) ctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGcc > 1:189004/1‑90 (MQ=255) aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTTCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGt > 1:518329/1‑90 (MQ=255) aaGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAAcc > 2:918744/1‑90 (MQ=255) tATCGGGTTGCAAAACGTGCCGGGGAAAAGGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:331592/90‑1 (MQ=255) tATCCGGTTGCAAACCGTGCGGCGAAAAAAGGTTCAAAATATTCCCACAACCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 2:248778/90‑1 (MQ=255) tATCAGTTTGCAAAGCGTCCGGGGGAAAAAGGTGCTAATTATTCCCACCCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 2:385855/90‑1 (MQ=255) tATCAGGTTGCGAAACGTGCCGGGGAAAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:718059/90‑1 (MQ=255) tATCAGGTTGCGAAACGTGCCGGGGAAAAAGGTGCAAAGTTTTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:543938/90‑1 (MQ=255) tATCAGGTTGCAAACCGTGCCGGGGAAAAAGGTGCTAATTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTc < 1:83897/90‑1 (MQ=255) ggTTGCGAAACGTGCCGGGGAAAAAGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:488006/90‑1 (MQ=255) ggTTGCAAAACGTGCCGCGGAAAAAGGTGCAAAGTATTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc < 1:798222/90‑1 (MQ=255) tGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGa > 2:7446/1‑90 (MQ=255) gcaaaaCGGGCCGCGAAAAAGGGTGTAAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGAt < 2:215264/87‑1 (MQ=255) cGAACCGTGCGGGGGAAAAGGGTGCAAAGTATTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc < 2:189004/90‑1 (MQ=255) cGAAACGTGCCGGGGAAAGAGGTGATAATTTTTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc < 2:513508/90‑1 (MQ=255) gTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGGGGGAATAACCTGACCGTCAGCGCAGATCTGtata > 1:120445/1‑90 (MQ=255) tGCCGGGGAAAAAGGTGCAAAGTATTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataa < 1:729533/90‑1 (MQ=255) tGCCGCGGAAAAAGGTGCTAATTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataa < 1:347230/90‑1 (MQ=255) cgcgGAAAAGGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata < 1:513374/90‑1 (MQ=255) agaaAGTTGCTAATTTTTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAt < 2:128142/90‑1 (MQ=255) agaaAGGTGCTAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAt < 1:728947/90‑1 (MQ=255) | GATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAT > NZ_CP009273/1673457‑1673632 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |