Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,348,810 | 0 | A | C | 51.7% | 33.4 / 42.0 | 29 | T243P (ACC→CCC) | gltB | glutamate synthase large subunit |
Reads supporting (aligned to +/- strand): ref base A (13/1); new base C (0/15); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.06e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCAC > NZ_CP009273/3348721‑3348884 | tGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAaca > 1:584964/1‑90 (MQ=255) caccaGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATc < 1:936565/90‑1 (MQ=255) cgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGt > 2:191266/1‑90 (MQ=255) cTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAAt > 1:821674/1‑90 (MQ=255) aaaCCCCGTCCCGCGCTGGCCGCGGGCGCACCCGTCCCGCTATCGGGCGCAAAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATg < 2:233199/89‑1 (MQ=255) aaCACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCTTCACCGGTAACCGCCAATgg > 2:167804/1‑90 (MQ=255) ccGTCCCGCGCGGGCCGCGGGCGCACCCTTTCCGCTTTCTGGCGCATAACGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGgcgc < 2:731405/90‑1 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGGGCATAACGGAGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:552181/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 1:598228/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 1:233199/1‑90 (MQ=255) cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg > 2:211565/1‑90 (MQ=255) gTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgt > 1:623768/1‑90 (MQ=255) cgcgcgGGCCGCGGGCGCACCCTTCCCGCTATCTGGCGCATAAGGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGgcgcgtgcgc < 1:72348/90‑1 (MQ=255) cgcgCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc > 1:86097/1‑90 (MQ=255) cgcgCTGGCCGCTGGCGCAACCGTTCCGCTATCCGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc > 2:455289/1‑90 (MQ=255) gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg > 1:428212/1‑90 (MQ=255) cgctggccgctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcgt > 2:547984/1‑90 (MQ=255) ctggccgcgggcCCAACCGTTCCGCTATCTGCCCAATAGCGGTGAAACCAACCCCATCACCGTTAACCGCCAATGGGCGCGTGCGCGTAc < 2:584964/90‑1 (MQ=255) cgggcCCACCCGTCCCGCTTTCGGGCGCATAGGGGAAAAACCAACCCCACCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATaaa < 1:609674/88‑1 (MQ=255) cccgttccgttATCGGGCGCATAAGGGTAAAACCAACCCCTCCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac < 1:167804/89‑1 (MQ=255) cccGTCCCCCTATCGGGCGCAAAAGGGTAAAACCACCCCCCCCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac < 1:547984/89‑1 (MQ=255) cttcccGCTATCGGGCGCAAAACGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac < 2:557511/86‑1 (MQ=255) ccttATCTGGCGCATAACGGAAAAACCAACCCCACCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:542463/87‑1 (MQ=255) cGCTTTCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:185528/90‑1 (MQ=255) cGCTATCGGGCGCATAACGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt < 2:578992/90‑1 (MQ=255) tCTGGCGCATAAGGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATccc < 2:178858/90‑1 (MQ=255) tGGCGCATAACGGTAAAATCAACCCCCTCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTg < 2:758580/90‑1 (MQ=255) gcgcATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAc > 2:929782/1‑90 (MQ=255) aagggTAAAACCACCCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 1:211565/87‑1 (MQ=255) aagggTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac < 2:348864/87‑1 (MQ=255) | TGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCAC > NZ_CP009273/3348721‑3348884 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |