Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,092,2250GT27.6% 96.8 / 31.1 59intergenic (+46/+872)PP_1865/aroHtransposase/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.

AAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT  >  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)
                                   tGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTTAATTTTGTCGTGGGCAAAGTAACCAAAACCCTCGGCTCCCAACATACCGCCCCAAcc                                                                                                >  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)
                                                     acGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTcc                                                                                <  1:391628/135‑1 (MQ=255)
                                                      cGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTTACTTTGGTCGTGGCCAAAGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCCGGGTCCCCt                                                                             <  1:149159/136‑1 (MQ=14)
                                                         gTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCt                                                                             <  1:156521/136‑1 (MQ=255)
                                                                       caaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACCAAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt                                                              <  1:527189/136‑1 (MQ=17)
                                                                       caaaaccaaaaGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCACTCCGGGCTTGCt                                                              <  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)
                                                                                                                                cccaaaGTAACCAAAACGGTCGGCTCCCATCATCCGGCCCCAACGCTGCGCTCCGGGGGCCCCTCGCGCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:515335/133‑1 (MQ=14)
                                                                                                                                 ccaaaTGAACCAAAACCGGTGGGGTCCCTTCTTCGGGCCCCTACGCTAGGGTCGGGGGTCCCGTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:65906/134‑1 (MQ=14)
                                                                                                                                 ccaaaGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:547357/134‑1 (MQ=255)
                                                                                                                                 ccaaaGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:326797/134‑1 (MQ=255)
                                                                                                                                 ccaaaGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTACGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:309167/134‑1 (MQ=255)
                                                                                                                                 caaaaGTAACCAAATCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:187475/135‑1 (MQ=255)
                                                                                                                                 caaaaGTAACCAAAACCGGCGGCTCCCATCATCCGGCCCCTACGTTGCGCTCGGGGGTCCCCTCGCTCCGGCCTTGCTCGCGGTAGCAGGCGCTGCTGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  2:291462/135‑1 (MQ=17)
                                                                                                                                 aaaatGTAACCAAAAACGGTCGGCTCCCATCATCCGGCCCCTACGCTGGGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  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)
                                                                                                                                 aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:146760/136‑1 (MQ=255)
                                                                                                                                 aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATcc    <  1:237050/136‑1 (MQ=255)
                                                                                                                                 aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                              >  2:99605/1‑79 (MQ=255)
                                                                                                                                 aaaaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGTCCCCTCGCTCCGGCCTTGCt                                                              <  1:99605/79‑1 (MQ=255)
                                                                                                                                  caaaGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:304562/135‑1 (MQ=255)
                                                                                                                                  caaaGTAACCAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCGCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:285742/135‑1 (MQ=255)
                                                                                                                                  caaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:142824/135‑1 (MQ=255)
                                                                                                                                  aaaaGTAACCAAAACGGTCGGCTCCCATAATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:236821/136‑1 (MQ=17)
                                                                                                                                  aaaaGTAACAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:238238/136‑1 (MQ=255)
                                                                                                                                  aaaaGGAACCAAAACCCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:367260/136‑1 (MQ=255)
                                                                                                                                  aaaaGAAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCa   <  1:37382/136‑1 (MQ=255)
                                                                                                                                   aaaGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCAACGCTGCGCACCGGGGGGCCCACCCGCCCGCCGTGCTCCCCGGGGCACGCCCAGCCCGCCCCCACCAAGGGGCCTCAGCGCGTTAACGGGctccca   >  2:142824/1‑131 (MQ=12)
                                                                                                                                   aaaCGAAACCAAAACCGCCGCCCCCCATCATCCGGCCCCTACGCTGCGCCCCGGGGCCCCCTCGCGCCGGCCTTGCCCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAt  <  1:528360/136‑1 (MQ=17)
                                                                                                                                       |                                                                                                                                    
AAAGCGAAGGCCTGGGTGCAAGGGACGGCTGTACATGCTCATGTTTCGACGGATCGGGTCATCAAAGCCGCCAAAACCAAAAGACCTTGAAAGGGATGGCCCTTACGGCGGGTCCCTTTTTGTCGTGGCAAAAAGGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCAT  >  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: