| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 16,736 | 0 | T | C | 48.5% | 18.2 / 12.8 | 33 | intergenic (+179/+15) | insL1/hokC | IS186 transposase/toxic membrane protein, small |
| Reads supporting (aligned to +/- strand): ref base T (6/11); new base C (0/16); total (6/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.84e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCA > NC_000913/16590‑16872 | ggAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTcccc < 2:185089/151‑1 (MQ=33) aGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCccgc < 1:333643/151‑1 (MQ=34) cTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCccgccg > 1:636395/1‑151 (MQ=35) tACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCccgcc < 2:391926/148‑1 (MQ=35) ggTGGCAATTCGGGGGACATTACCCCTGACAGCCAACAAAACGCAATTAAAGAACGCGAGGCATCGTCTTAACGAGGCCCCGAGGCGTCGCTTTCTTCAGATGG‑TTCACCCTTT‑AA‑TTTAGCGCTTATGGGATCCCCCCCCGCCGTTGct < 2:453981/150‑1 (MQ=11) cATTCGTGGGACAGTACCCCTGACAGCCTACAAAACGAAATTGAAGAACGCGAGGCATCGTTTTAACGAGGCACCGAGGCGTCGCATTCTCCAGAGGG‑TTCACCCCTT‑AA‑GTTAGCGCTTATGGGATCACCCCCCGCCGTTGCTCTTACTc < 1:522779/151‑1 (MQ=14) tCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCgg > 1:378477/1‑150 (MQ=255) aCAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAAc > 2:297734/1‑150 (MQ=255) ccTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCa > 1:338970/1‑150 (MQ=255) aaaaCGCAATGGAAAAACGCGAGGCATCGTTTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGT‑TTCAACCTTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTc < 2:209347/150‑1 (MQ=255) agcacAATAGAAAAGCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCACCCTTT‑AA‑TTTAGCGCTTATGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGt < 2:128037/147‑1 (MQ=21) aaTTGAAGAACGCGGGGCATCGTCTTAACGGGGCCCGGAGGCGTCGCTTTCTTCAGAGGT‑TTCAACCCTT‑AA‑TTTAGCGTTTAGGGGATCACCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCt < 2:286739/148‑1 (MQ=255) tGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCACCCCTT‑AA‑GTAAGCGCTTAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc < 1:486321/149‑1 (MQ=255) aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTCCAGATGG‑TTCACCCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCgg < 1:197454/151‑1 (MQ=255) gccATCTTCTTACCGAGCCACCGAGGCGTCCCATTCTTCAGATGG‑TTCACCCCTT‑AA‑GTTAGCGCTTATGGAATCACCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGAtgt < 2:58859/146‑1 (MQ=255) tAACGAGGCACCGAGGCGCCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTAGGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACaga < 1:39385/150‑1 (MQ=255) aaCGGGGCCCGGGGGCGTCGCATCCTAATGAGGG‑TTCACCCCTA‑AA‑GTGAGCGCGTGGGTGATCCCCCCCCGCCGTTGCGCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTcaca < 2:349192/147‑1 (MQ=255) actgggCCCCGAGGCGCCGCATCCTGCAGAGGG‑TTAAACCTTT‑AA‑GTCAGCTCTTCGTGACACCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACag < 2:399933/143‑1 (MQ=255) cGAGGCCCCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTAGGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTc < 1:147983/151‑1 (MQ=255) cACCGAGGCGCCGCATTCTTCAGATGG‑TTCACCCCTT‑AA‑GTTAGCGTTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCtttt < 1:474181/150‑1 (MQ=255) cgggtcGTCGCTTTTTTC‑GATGGTTTCACCCCTT‑AA‑GTTAGCCCTTATGGGATCGCTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTtctc < 2:377655/145‑1 (MQ=255) gggCGCCGCTTTCTTCAGATGT‑TACACCCGTT‑AA‑GTTAGCGCTGATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGtt < 2:653311/150‑1 (MQ=255) ggCGTGCCATCCTGCAGTTGA‑TTCACCCCTT‑AG‑TTGAGCTTTTGTGGCACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTa < 2:327101/151‑1 (MQ=255) cccATTCTTCAGATGG‑TTCACCCCTA‑AT‑GTAAGCGCTTATGGGATCACCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAgc < 2:631792/149‑1 (MQ=255) tcAGATGG‑TTCACCCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGcgccgc < 1:547242/148‑1 (MQ=255) cAGATGG‑TTCACCCCTT‑AA‑GTTAGCGCTTAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCgccgcc < 2:466171/148‑1 (MQ=255) ttCCACCCTT‑AATGTTGGCTTTTAGGGCACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTAAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTa < 2:530665/146‑1 (MQ=255) tccaccccTT‑AA‑TTTAGCGTTTAGGGGACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTAcggcg < 1:245672/145‑1 (MQ=255) tCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTAcggcgg > 2:520084/1‑150 (MQ=255) aCCCTTAAA‑GTTAGCTTTTAGGGAACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTAAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGt < 1:674141/149‑1 (MQ=255) cccttaaatttggctcTTAGGGGACCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGAt < 2:120437/136‑1 (MQ=255) gcTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATc > 2:359163/1‑150 (MQ=255) ttAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCa < 1:293789/149‑1 (MQ=255) | GGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGG‑TTCAACCCTT‑AA‑GTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCA > NC_000913/16590‑16872 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |