Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,975,994 | 0 | C | G | 53.8% | 8.4 / 13.7 | 13 | intergenic (‑126/+41) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
Reads supporting (aligned to +/- strand): ref base C (1/5); new base G (7/0); total (8/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGA > NZ_CP009273/1975912‑1976072 | cgcTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGGAtgggg > 2:362721/1‑88 (MQ=255) gTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCGATGGTGGGGGAAAAAGccga > 1:206708/1‑87 (MQ=255) gACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAGCTAAAATTTTGCa > 1:131618/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTTGGTTTGGGGGGGAAAAAAGCCAAAATGTTTCAGaaa > 1:623537/1‑90 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGGTTGGGGGGGGAAAAAAACAAATTGTTGAAGaaa > 2:6429/1‑90 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTTTGGGGGGGAAAAAAACTAAAATTTTTTAGaaaa > 1:182311/1‑90 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGGGGGTTTGGGGGGGGAAAAAGATAAAATTTTTCAGaaaa > 2:493411/1‑90 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 1:378836/90‑1 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGaaaaa < 2:933990/90‑1 (MQ=255) ggTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAACTAAATTGTTGGAGAAAAAAACAACAg > 1:399297/1‑90 (MQ=255) aaCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTaa < 1:913301/90‑1 (MQ=255) aCCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAAc < 1:902433/90‑1 (MQ=255) aTGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGa < 1:493411/90‑1 (MQ=255) | CGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGA > NZ_CP009273/1975912‑1976072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |