| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,976,032 | 0 | G | A | 44.0% | 20.0 / 14.3 | 25 | intergenic (‑164/+3) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
| Reads supporting (aligned to +/- strand): ref base G (2/12); new base A (11/0); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.02e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCC > NZ_CP009273/1975898‑1976163 | tgtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCAt < 1:773644/138‑1 (MQ=255) aCAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCttt < 2:806111/138‑1 (MQ=255) tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGGCT‑‑GG‑TGGGGAAAAAAACTAAAATTTTGCAGAAAAAAAAAACAAAAAGAGCGGTTTGCTTTTTTTTTGt > 1:842751/1‑139 (MQ=255) tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGTTTT‑‑TG‑TGGGGAAAAAAGCCAAAATTTTTCAGAAAAAAAGAACAGACAGGGGGGTTCGCTTTTTTCttttta > 2:111653/1‑135 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTATGGGG‑GGGGAAAAAACCCAAAATT‑TT‑CAGAAAAAAAAACCAAAAAGAGGGGTTTGCTTTTTTTTTTTAAAGGGCGAAGaac > 2:668270/1‑136 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTAT‑‑TG‑TGGGGAAAAAAGCCAAATTTTTGTAGAAAAAAACATCAAACAATGCGGTTTGCTTTTTTCTTGTAAAGGGCGGAGAACa > 1:527141/1‑139 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTAT‑‑GG‑TGGGGAAAAAAGCTTAAATTTTTCAGAAAAAAAAATCAGACAATGTGGTTTGCTTTTTTCTTTTAAAGGTCGAAGAtta > 1:847371/1‑137 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTAT‑‑GG‑TGGGGAAAAAAACTAAAATGTTGGAGAAAAAAGAAACAGACAGTGGGGTTTGCTTTTTTCTTGTTACGGGCGAAGAACa > 2:154199/1‑139 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTAT‑‑GG‑GGGGGAAAAAAGCCAAATTTTTTCAGAAAAAAACATCAGAAAGTGCGGTTCGCTTTTTTCTTTTAACGGGCGAAGaac > 1:343401/1‑136 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGGAAAAAAGCTAAATTGTTTCAGAAAAAAACAACAGAAAGTGGGGTTCGCTTTTTTTTTTTAAAGGTCGAAGaac > 1:832499/1‑136 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGGTGTGCAT‑‑GGCTGGGAAAAACAGCTAAA‑TATTGCTCAGAAAAACATCACACAGAGCGGTTTGCCGTTTTCGTGTTACGGGCTAAGAtag > 2:62769/1‑137 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGATGGGTAT‑‑GG‑GGGGGAAAAACAACAAAATATTTTAGAAAAAAGGACCAGAAAGGGGGGTTTGCTTTTTGTTTGTTAAGGGCGagacaca > 2:636301/1‑133 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGGAT‑‑GG‑GGGGGAAAAAAACAAAAATTTTGGAGAAAAAAAAAACAAACAGGGCGGTTTGCGTTTTTTTTGTTAAGGGCGAAGaaaac > 2:345969/1‑135 (MQ=255) gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGCTTG‑‑GG‑GGGGAAAAAAAACAAAATT‑TTGTAGAAAAAAAAAATAAAAAGAGGGGGTTGCTTTTTTTTTGTAAAGGTCGAAAAACACCGcccgcca > 2:837747/1‑133 (MQ=255) ttGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGGTAT‑‑GG‑GGGGGAAAAAAACTAAAATGTTTTAGAAAAAAACATCACACAGAGGGTTTTGCTTTTTTTTTTTAACGGTCGAAGAACAACGCTGCCAGCAggggt > 2:733483/1‑135 (MQ=255) cGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 1:154199/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 2:842751/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 2:572410/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 2:305382/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 1:868672/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 1:345969/139‑1 (MQ=255) gC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAg < 1:163293/138‑1 (MQ=255) ggTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCATAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCAAAGATGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGgag < 1:668453/139‑1 (MQ=255) aaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATg < 2:343401/139‑1 (MQ=255) gAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGcc < 1:656945/139‑1 (MQ=255) | TGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGGCGTTTTCCACGGTGCC > NZ_CP009273/1975898‑1976163 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |