| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,092,230 | 0 | A | . | 33.3% | 41.4 / 23.1 | 22 | intergenic (+51/+867) | PP_1865/aroH | transposase/3‑deoxy‑7‑phosphoheptulonate synthase |
| Reads supporting (aligned to +/- strand): ref base A (4/10); new base . (0/7); total (4/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.55e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT > NC_002947/2092107‑2092356 | gCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGc > 1:4021/1‑136 (MQ=255) gCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGCGTCCCATTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCcaacaac > 2:302965/1‑134 (MQ=255) cTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATcc > 2:300924/1‑136 (MQ=255) gTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCgg < 1:140186/136‑1 (MQ=255) gACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCCACGCAGCGCGcc > 2:326525/1‑136 (MQ=255) ggTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCGGTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTcccc < 1:302965/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCTAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTTCGCTCCGGGGTCCCCTCACTCCGGGCTTGCt < 1:219122/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGGCTTGCt < 2:94526/136‑1 (MQ=255) aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt < 1:300924/136‑1 (MQ=255) aaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTcc < 1:342827/136‑1 (MQ=255) aGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGTGAGc < 2:277176/136‑1 (MQ=255) ccTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGGA‑cc < 1:243996/136‑2 (MQ=255) ttGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑cgcg < 1:133064/136‑1 (MQ=255) ggTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGGACCGCGCTG‑AGGCCCCATGGATGGGGCCt < 1:204920/136‑1 (MQ=11) gccaaaGTAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:122276/133‑1 (MQ=14) gccaaaGTAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:380860/133‑1 (MQ=14) gccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTGCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:201398/133‑1 (MQ=14) gccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:159611/133‑1 (MQ=255) taaaaGTAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc < 1:163651/135‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:388459/134‑1 (MQ=255) ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa < 1:130072/134‑1 (MQ=255) aaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGCCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt < 1:159606/136‑1 (MQ=255) | GCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCA‑CGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT > NC_002947/2092107‑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: ‑ |