Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,976,010 | 0 | C | A | 40.0% | 37.4 / 11.5 | 26 | intergenic (‑142/+25) | otsB/araH | trehalose‑phosphatase/arabinose ABC transporter permease AraH |
Reads supporting (aligned to +/- strand): ref base C (3/12); new base A (10/0); total (14/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTT > NZ_CP009273/1975883‑1976142 | tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTGT‑GGGGGGGAAAAAAGACAAAatgtt > 2:25302/1‑139 (MQ=255) aTTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGaaaa < 1:319474/136‑1 (MQ=255) gtgtTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGGGGGGAAAAAAGACAAAATTTTTG‑AGaaaaaaaaa > 2:319474/1‑135 (MQ=255) gcgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGa < 1:480980/138‑1 (MQ=255) cgcAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTGT‑GGGGGGGAAAAAAACAAAAAATATGCAGAAAAAAGAATCAGAc > 2:108147/1‑138 (MQ=255) aCAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑ttt < 1:25302/138‑1 (MQ=255) tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAT‑GGGGGGGAAAAAAGCTAAAATGTTGTAGAAAAAAAAAACAGACAAAGCGTT‑‑TTGT‑TTTTTTTTTGTa > 2:258604/1‑138 (MQ=255) aGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTTT‑TGTGGGGGAAAAAGCAAAAATGTTGTAGTAAAAAACAAAAGACAATGCGGT‑‑TTGCTTTTTTGTTTGTaaggtgg > 2:411912/1‑133 (MQ=255) aGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGGTT‑GGTGGGGAAAAAAGCTAAATTGTTGCAGAAAAAAAAAAAAAACAGTGCGGT‑‑TTGC‑TTTTTTTTTGTAACGGGCGAAGaa > 2:116968/1‑138 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTAT‑GGGGGGGAAAAAAGACAAAATTTTGGAGAAAAAAGCATCAAAAAGGGCGGT‑‑TCGC‑TTTTTTTTTGTTTAGGTCGAAGaac > 2:121841/1‑136 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGGTT‑GGTGGGGGAAAAAGCTAAAATGTTGCTGAAAAAAAAAACAGACAGTGGGGT‑‑TTGC‑TTTTTTGTTTTAAAGGGCGCAGaata > 2:293074/1‑136 (MQ=255) tgtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑TGGGGGGAAAAAAGCCAAAATGTTGCAGAAAAAAACATCAGACAGTGCGTT‑‑TTGC‑TTTTTTCTTGTAAAGGGCGGAAATCa > 1:329020/1‑139 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGGAT‑GGGGGGGAAAAAAGCAAAAATTTTTTTGAAAAAAAAATAAGACAATGGGGT‑‑TTTC‑TTTTTGTTTTTTAAGGTCGAAGaaaaa > 2:429911/1‑135 (MQ=255) gtgACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGGGGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACAACAGACAGAGCGGT‑‑TTGC‑TTTTTTCTTTTAAAGGGCGAAGATCAc > 1:494255/1‑139 (MQ=255) cGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCt < 1:136390/139‑1 (MQ=255) tttGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAtt < 1:116968/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 1:95088/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCcgc < 1:48537/139‑1 (MQ=255) cGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGACAGCAGGATTAAGCcgc < 1:614312/139‑1 (MQ=255) ccAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTa < 1:121841/139‑1 (MQ=255) gtggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAg < 1:258604/138‑1 (MQ=255) tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCTCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAg < 1:108147/139‑1 (MQ=255) tggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAg > 1:269225/1‑137 (MQ=255) ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa < 2:14401/120‑1 (MQ=255) ggAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGaa > 1:14401/1‑120 (MQ=255) aaaaaaCGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGtt > 2:393636/1‑139 (MQ=255) | TCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTAT‑GGTGGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTT‑‑TCGC‑TTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTT > NZ_CP009273/1975883‑1976142 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |