Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,089,2030A.30.0% 59.7 / 31.4 30intergenic (+51/+867)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (3/18);  new base . (0/9);  total (3/27)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT  >  NC_002947/2089085‑2089329
                                                                                                                      |                                                                                                                              
ggACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCca                                                                                                               >  2:252862/1‑136 (MQ=255)
ggACGGCTGTACATGCTCATGTTTCGACGGAACGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCca                                                                                                               >  2:371251/1‑136 (MQ=255)
             tGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCt                                                                                                  <  1:165876/136‑1 (MQ=255)
                                                 caaaacccaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCGACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt                                                             <  1:424854/136‑1 (MQ=17)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  1:133982/136‑1 (MQ=255)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  2:199068/136‑1 (MQ=255)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  1:2031/136‑1 (MQ=255)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  1:441639/136‑1 (MQ=255)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  1:387809/136‑1 (MQ=255)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  1:371251/136‑1 (MQ=255)
                                                  aaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt                                                             <  1:252862/136‑1 (MQ=255)
                                                              ccTTGAAAGGGATGGCCCTTACGGCTGGTCCCTTTTTGTCGTGGCAAAAAGTAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCGGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGGAcc                                                 <  2:399556/136‑2 (MQ=255)
                                                                               cTTCCGGCGGGTCCC‑TTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCccatcca                               <  1:83901/136‑1 (MQ=255)
                                                                                       gggTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGcc                        <  1:224461/136‑1 (MQ=255)
                                                                                                  tGTCGTGGAAAAAAGGAACAAAAAACGGTCGGCTCCCATCATCGGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCTTGGGGCCTCAGCGCTTGa             <  1:257571/136‑1 (MQ=17)
                                                                                                          gccaaaGTAACC‑AAACCCGGCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:106440/133‑1 (MQ=14)
                                                                                                           ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:38903/134‑1 (MQ=255)
                                                                                                           ccaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:111634/134‑1 (MQ=255)
                                                                                                           ccaaaGAAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:87913/134‑1 (MQ=255)
                                                                                                           ccaaaGAAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:135733/134‑1 (MQ=255)
                                                                                                           caaaaGAAACC‑AAAACCGGCGGCCCCCATCATCCGCCCCCTACGCGGCGCTCCGGGGTCCCCGCGCCCCGGCCTTGCTCCCGGGCGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:258688/135‑1 (MQ=17)
                                                                                                           aaaaaGTAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCGGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:397570/136‑1 (MQ=255)
                                                                                                           aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCGGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             >  2:289042/1‑79 (MQ=255)
                                                                                                           aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCGGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                             <  1:289042/79‑1 (MQ=255)
                                                                                                            caaaGTAACC‑AAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:88349/135‑1 (MQ=255)
                                                                                                            aaaaGTAACC‑AAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:177877/136‑1 (MQ=255)
                                                                                                            aaaaGTAACC‑AAAACCGCCGGCTCCCATCAACCGCCCCCTACGCTGCGCGCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:78130/136‑1 (MQ=14)
                                                                                                            aaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:414711/136‑1 (MQ=255)
                                                                                                            aaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGGCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:238550/136‑1 (MQ=255)
                                                                                                             aaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCGGCGCGCCGGGGGCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:342508/136‑1 (MQ=255)
                                                                                                                      |                                                                                                                              
GGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT  >  NC_002947/2089085‑2089329

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: ‑