Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,092,225 | 0 | G | T | 27.6% | 96.8 / 31.1 | 59 | intergenic (+46/+872) | PP_1865/aroH | transposase/3‑deoxy‑7‑phosphoheptulonate synthase |
Reads supporting (aligned to +/- strand): ref base G (13/29); new base T (1/15); total (14/45) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.36e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAT > NC_002947/2092090‑2092356 | aaaGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCATTTTGTCGTGGCAAAAAgg > 1:251089/1‑136 (MQ=255) ccTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCATTCGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCATAAAc > 2:390925/1‑136 (MQ=255) cTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAAcc > 2:92167/1‑136 (MQ=255) ggACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCGTAAGGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCca > 2:419497/1‑136 (MQ=255) ggACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCca > 2:156521/1‑136 (MQ=255) gACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCcat > 1:69705/1‑136 (MQ=255) gCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCatcatc > 2:121345/1‑136 (MQ=255) gCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCatcatc > 2:118559/1‑136 (MQ=255) aCATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGcc < 2:321428/136‑1 (MQ=255) tGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGT‑TAATTTTGTCGTGG‑GCAAAGTAACC‑AAAACCCTCGGCTCCCAACATACCGCCCCAAcc > 2:149159/1‑135 (MQ=14) cTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTAc < 1:102794/136‑1 (MQ=255) cGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGTATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCgg < 2:220018/136‑1 (MQ=255) cGGATCGGGTCATAAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTGCGCTACGCTCCgg < 2:69705/136‑1 (MQ=255) tCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTc < 1:562271/136‑1 (MQ=255) acGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTcc < 1:391628/135‑1 (MQ=255) cGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTTAC‑TTTGGTCGTGGC‑CAAAGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCCGGGTCCCCt < 1:149159/136‑1 (MQ=14) gTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCt < 1:156521/136‑1 (MQ=255) caaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:527189/136‑1 (MQ=17) caaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:529860/136‑1 (MQ=17) aaatccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCTTCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:150882/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:190311/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:258388/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:286062/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:57111/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:173439/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:92167/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGGCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:474626/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGGCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:448812/136‑1 (MQ=255) agccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTcc < 1:545482/134‑1 (MQ=255) aaaGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTg < 2:311348/136‑1 (MQ=255) aGGGATGGCCGTAAGGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCa < 1:419497/136‑1 (MQ=255) aTGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGcc < 2:21927/136‑1 (MQ=255) ggcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCGCGCTCCGGCCGTGCTCCCGGGGGCACGCGCTGCAGGCCCCATCCATggg > 2:367260/1‑136 (MQ=255) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:343803/84‑1 (MQ=255) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt > 2:343803/1‑84 (MQ=255) cccaaaGTAACC‑AAAACGGTCGGCTCCCATCATCCGGCCCCAACGCTGCGCTCCGGGGGCCCCTCGCGCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATcc < 1:515335/133‑1 (MQ=14) ccaaaTGAACCAAAACCGGTGGGGTCCCTTCTTCGGGCCCCTACGCTAGGGTCGGGGGTCCCGTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATcc < 1:65906/134‑1 (MQ=14) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:547357/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:326797/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:309167/134‑1 (MQ=255) caaaaGTAACC‑AAATCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:187475/135‑1 (MQ=255) caaaaGTAACC‑AAAACCGGCGGCTCCCATCATCCGGCCCCTACGTTGCGCTCGGGGGTCCCCTCGCTCCGGCCTTGCTCGCGGTAGCAGGCGCTGCTGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 2:291462/135‑1 (MQ=17) aaaatGTAACCAAAAACGGTCGGCTCCCATCATCCGGCCCCTACGCTGGGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATcc < 2:70503/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:145429/79‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt > 2:396212/1‑79 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt > 2:145429/1‑79 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:396212/79‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATcc < 1:146760/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATcc < 1:237050/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGTCCCCTCGCTCCGGCCTTGCt > 2:99605/1‑79 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:99605/79‑1 (MQ=255) caaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAt < 1:304562/135‑1 (MQ=255) caaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCGCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAt < 1:285742/135‑1 (MQ=255) caaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:142824/135‑1 (MQ=255) aaaaGTAACC‑AAAACGGTCGGCTCCCATAATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAt < 1:236821/136‑1 (MQ=17) aaaaGTAA‑CAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAt < 1:238238/136‑1 (MQ=255) aaaaGGAACCAAAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:367260/136‑1 (MQ=255) aaaaGAAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCa < 1:37382/136‑1 (MQ=255) aaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCAACGCTGCGCACCGGGGGGCCCACCCGCCCGCCGTGCTCCCCGGGGCACGCCCAGCCCGCCCCCACCAAGGGGCCTCAGCGCGTTAACGGGctccca > 2:142824/1‑131 (MQ=12) aaaCGAAACCAAAACCGCCGCCCCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGCCCCCTCGCGCCGGCCTTGCCCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAt < 1:528360/136‑1 (MQ=17) | AAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGC‑TTGACGGGCATCCAT > NC_002947/2092090‑2092356 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |