Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,092,225 | 0 | G | T | 26.4% | 96.9 / 12.4 | 55 | intergenic (+46/+872) | PP_1865/aroH | transposase/3‑deoxy‑7‑phosphoheptulonate synthase |
Reads supporting (aligned to +/- strand): ref base G (11/28); new base T (2/12); total (13/42) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.73e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT > NC_002947/2092090‑2092356 | aaaGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAgg > 2:294925/1‑136 (MQ=255) gCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGc > 2:115224/1‑136 (MQ=255) gCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGc > 1:400848/1‑136 (MQ=255) ggACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAAACCGTCGGCTCCcat > 2:386329/1‑136 (MQ=255) ggACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCca > 1:141727/1‑136 (MQ=255) gACGGCTGTACATGCTCATGTTTCGACGGATCGTGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAAACCGTCGGCTCCcata > 2:534549/1‑135 (MQ=255) tGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAAc < 2:324054/117‑1 (MQ=255) tGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAAc > 1:324054/1‑117 (MQ=255) gTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCgg < 1:115224/136‑1 (MQ=255) gTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCgg < 1:480964/136‑1 (MQ=255) gtCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTAc < 2:180441/135‑1 (MQ=255) cTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTAc < 1:87119/136‑1 (MQ=255) cTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTAc < 2:400848/136‑1 (MQ=255) tttttGAGGGATCGGGTCTTCATAGCCTCCAAAACCAATATTCTTTGAAAGGGATGGCCCTTACTGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAATTCCGTCGGCTCCCATCATCCGGCCCCTACGCTgc < 2:211321/135‑1 (MQ=18) tCGACGGATCGGGTCATCTAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCt < 2:418746/136‑1 (MQ=255) aCGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCg < 2:338998/136‑1 (MQ=255) tCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTc < 1:175492/136‑1 (MQ=255) gTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCt < 2:290494/136‑1 (MQ=255) caaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:386329/136‑1 (MQ=17) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:329543/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 2:207450/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:172875/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:447422/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:10044/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGC‑CAAAGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCTcc < 1:149733/136‑1 (MQ=14) aaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTcc < 1:400401/136‑1 (MQ=255) aaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTGCGCTACGCTCCGGGGTCCCCTCGCTCCGGGCTTGCTCCCGGGa < 1:437455/136‑1 (MQ=255) aGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGc < 1:125205/136‑1 (MQ=255) gcggGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCGACGCGGCGCGGCGGGGTGCCCGCGCTCCGGCCGTGCTCCCGGGCGCACCCGCAGCCGG‑CCCCCACCAAgggg > 2:332891/1‑136 (MQ=14) ggTCCCTTTTTGTCGTGCAAAAAAGTAACAAAAAACCGTCGGCTCCCATCATCCGGCCCCTCCGCGGCGCTCCGGGGGCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCt < 1:64598/136‑1 (MQ=17) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCATGCTCCCGGGAGCACGCGCTGCAGG‑CCCCAACCATGGGGCCAAAGCGCATGaacgg > 2:328100/1‑132 (MQ=17) gTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCGACGCGGCGCGCCGGGGGGCCCACGCGCCGGCCGTTCTACCGGGCGCACCCCCAGCCGG‑CCCCCACCAAGGGGGCGCCGCGCATGTCggg > 2:498823/1‑136 (MQ=9) cccaaaGTAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCGGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:397023/133‑1 (MQ=14) gcaaaaGAACCCAAACCCGTCGGCTCCCATCATCCGCCCCCAACGCTGCGCTCCGGGGCCCCCTCGCGCCGGCCTTGCCCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:374035/134‑1 (MQ=17) cccaacgAACCCAAACCCGCCGCCACCCATCACCCGCCCCCTACGCGGCGCGCCGGGGCCCCCTCGCCCCGGCCTTGCCCCCGGGCGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:51664/130‑1 (MQ=17) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:105050/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:252357/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:394507/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:46375/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGGCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:170656/134‑1 (MQ=255) aaaatGTAACAAAAACCCGTCGGCTCCCATCATCCGGCCCCTACGGTGCGCGCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:9141/136‑1 (MQ=17) aaaaaGGAACCAAAACCCGGCGGCCCCCATCATCCGGCCCCTACGCTGCGCTCGGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:252355/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:208353/79‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt > 2:208353/1‑79 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:332891/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:48222/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:498823/136‑1 (MQ=255) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCAACGCTGCGCTGCGGGGGGCCCACGCTCCGGCCGTTCTGCCGGGGGCACCCGCGGCAGGCCCCCAACAATGGGG‑CGCAGCGCGTGAACGGCATcc > 2:15487/1‑136 (MQ=9) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCAACGCTGCGCTCCGGGGGGCCCACGCGCCGGCCGTGCTCCCCGGAGCACCCCCTGCCCG‑CCCCCAACAAGGGGCCGGCGCGCTTGACGGGCATcc > 2:299304/1‑136 (MQ=9) aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCAACGCGGCGCGCCGGGGGGCCCCCGCGCCGGCCGTGCTCCCGGGCGCACGCCCGGCCGG‑CCCCCCCCCAGGGGGCGCCGCGCTTGACggggctcc > 2:252355/1‑131 (MQ=12) aaaaaGGAACCAAAAACCGGCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:502535/136‑1 (MQ=255) caaaGAAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt < 1:253401/135‑1 (MQ=255) aaaaGTAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGCCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt < 1:166085/136‑1 (MQ=17) aaaaGTAACAAAAAACGGTCGGCTCCCATCATCCGGCCCCTCCGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:84916/136‑1 (MQ=255) aaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt < 1:299304/136‑1 (MQ=255) aaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt < 1:15487/136‑1 (MQ=255) aaaGCAACAAAAACCCGTCGGCTCCAATCATCCGGCCCCTACGCTGCGCTCCGGGGGCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt < 1:200123/136‑1 (MQ=17) | AAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGG‑CCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT > 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: ‑ |