| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,976,010 | 0 | C | A | 42.9% | 27.9 / 18.5 | 28 | intergenic (‑142/+25) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
| Reads supporting (aligned to +/- strand): ref base C (3/13); new base A (12/0); total (15/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATG > NZ_CP009273/1975883‑1976146 | tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAATGAAGAACAGCCAACCGTTATGGC‑GGTTTGCGCT‑‑GG‑GGGGGAAAAAAACTaaaaattt > 1:287416/1‑134 (MQ=255) gttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGaaaa < 1:346643/138‑1 (MQ=255) 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) tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAAttgttgta > 1:767239/1‑96 (MQ=255) tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCa < 2:767239/98‑1 (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) | TCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGC‑GGTTGGCTAT‑‑GG‑TGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATG > NZ_CP009273/1975883‑1976146 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |