| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 19 | 0 | C | A | 25.0% | 66.7 / 14.2 | 28 | intergenic (–/‑171) | –/thrL | –/thr operon leader peptide |
| Reads supporting (aligned to +/- strand): ref base C (12/9); new base A (4/3); total (16/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAC > NC_000913/1‑137 | acgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGt > 1:516718‑M1/107‑139 (MQ=255) cgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgctattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAAc > 2:409375‑M1/106‑127 (MQ=255) cgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAAc < 1:409375‑M1/22‑1 (MQ=255) cgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCa > 2:2405‑M1/106‑132 (MQ=255) cgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCa < 1:2405‑M1/27‑1 (MQ=255) cgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGtc < 1:335540‑M1/34‑1 (MQ=255) cgcctggggcttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAACTTATCATTCTGACTGCAACGGGCAATATGtc < 2:394687‑M1/34‑1 (MQ=255) cttttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTgg > 1:165532‑M1/97‑139 (MQ=255) tttagagcaacgagacacggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaagtaaaaaacgccttagtaagtatttctcAGCCTTGCATTCTGACTGAAACGGGAAATATGTCTCTGTGTGGAt > 2:327657‑M1/95‑139 (MQ=255) cggcaatgttgcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTg < 2:681892‑M1/62‑1 (MQ=255) gcaccgtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCtt > 1:681892‑M1/68‑139 (MQ=255) gtttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGaa > 1:394687‑M1/63‑139 (MQ=255) tttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGAAACGGGCAATATGTCTCTGTATGGATTAAAAAAAGAGTGTCTGATagc < 1:200027‑M1/67‑1 (MQ=255) tttgctgcatgatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGAAACGGGCAATATGTCTCTGTATGGATTAAAAAAAGAGTGTCTGATagc > 2:200027‑M1/62‑128 (MQ=255) gatattgaaaaaaatatcaccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTg < 2:357301‑M1/88‑1 (MQ=255) caccaaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTaaaa < 2:287564‑M1/105‑1 (MQ=255) aaataaaaaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAAtttt > 1:14633‑M1/31‑139 (MQ=255) aaacgccttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACt > 1:357301‑M1/24‑139 (MQ=255) ttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGaa < 1:234678‑M1/77‑1 (MQ=255) ttagtaagtatttttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGaa > 2:234678‑M1/17‑93 (MQ=255) ttcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTaa < 1:72268‑M1/136‑1 (MQ=255) tcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAc > 2:533368‑M1/3‑139 (MQ=255) tcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAc > 1:287564‑M1/3‑139 (MQ=255) tcAGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAc > 2:72268‑M1/3‑139 (MQ=255) tttCATTCTGACTGAAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAc < 1:492032/133‑1 (MQ=255) tttCATTCTGACTGAAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAc > 2:492032/1‑133 (MQ=255) ttCATTCTGACTGAAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGcc > 2:133818/1‑85 (MQ=255) ttCATTCTGACTGAAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGcc < 1:133818/85‑1 (MQ=255) | ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCTTTTCATTCTGACTGCAACGGGCAATATGTCTCTGTGTGGATTAAAAAAAGAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAAC > NC_000913/1‑137 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |