Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,716,5870AG63.9% 44.4 / 70.7 61G70G (GGA→GGG) ydhKFUSC family protein
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (21/18);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TACCATTGCCATGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGCCTGTGCGCACTTCAC  >  NZ_CP009273/1716458‑1716721
                                                                                                                                 |                                                                                                                                      
tACCATTGCCATGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCt                                                                                                                               >  2:657949/1‑139 (MQ=255)
     ttGCCATGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAg                                                                                                                          >  1:32380/1‑139 (MQ=255)
     ttGCCATGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAg                                                                                                                          >  1:1281986/1‑139 (MQ=255)
        ccATGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGttt                                                                                                                       >  2:1084100/1‑139 (MQ=255)
          aTGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGc                                                                                                                     >  1:211329/1‑139 (MQ=255)
          aTGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGc                                                                                                                     >  1:278001/1‑139 (MQ=255)
           tgtgTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCt                                                                                                                    <  2:878548/139‑1 (MQ=255)
            gtgtCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTc                                                                                                                   >  1:390716/1‑139 (MQ=255)
                  ggCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGc                                                                                                              >  1:761747/1‑138 (MQ=255)
                  ggCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGc                                                                                                              >  2:1075980/1‑138 (MQ=255)
                  ggCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGc                                                                                                              >  1:1087939/1‑138 (MQ=255)
                      cTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTg                                                                                                         <  1:22728/139‑1 (MQ=255)
                         aCGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCgg                                                                                                      <  1:1084100/139‑1 (MQ=255)
                           ggTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCa                                                                                                    <  2:98977/139‑1 (MQ=255)
                                       tttAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGc                                                                                         <  2:63054/138‑1 (MQ=255)
                                              cTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATa                                                                                 >  2:352620/1‑139 (MQ=255)
                                                 gATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGc                                                                              >  2:166807/1‑139 (MQ=255)
                                                 gATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGcc                                                                                        <  2:668957/129‑1 (MQ=255)
                                                 gATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGcc                                                                                        >  1:668957/1‑129 (MQ=255)
                                                       cccTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAAt                                                                         >  2:200842/1‑138 (MQ=255)
                                                               ggCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTgg                                                                >  2:363705/1‑139 (MQ=255)
                                                                 cGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGtt                                                              >  2:795765/1‑139 (MQ=255)
                                                                  gATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTg                                                                 >  2:907433/1‑135 (MQ=255)
                                                                  gATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTg                                                                 <  1:907433/135‑1 (MQ=255)
                                                                   aTGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGtttt                                                            <  2:609223/139‑1 (MQ=255)
                                                                    tGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGc                                                                                         <  2:1049458/109‑1 (MQ=255)
                                                                    tGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGc                                                                                         >  1:1049458/1‑109 (MQ=255)
                                                                    tGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAAt                                                                         >  2:495278/1‑125 (MQ=255)
                                                                    tGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAAt                                                                         <  1:495278/125‑1 (MQ=255)
                                                                    tGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGAGTTATCAGCAAAAGACTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGAGGAACTGCTTCTCGCCGGGCATACGCTCGATGAGCCGTGGttttt                                                           >  2:263100/1‑139 (MQ=255)
                                                                      aCCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGtttttt                                                          <  2:118653/138‑1 (MQ=255)
                                                                      aCCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTc                                                         >  2:687398/1‑139 (MQ=255)
                                                                                cAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCAt                                               >  1:905405/1‑139 (MQ=255)
                                                                                 aGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGc                                                                                         <  1:512287/96‑1 (MQ=255)
                                                                                 aGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGc                                                                                         >  2:512287/1‑96 (MQ=255)
                                                                                     gTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCgg                                          >  1:702547/1‑139 (MQ=255)
                                                                                      ttAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGc                                         <  1:1254402/139‑1 (MQ=255)
                                                                                      ttAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGc                                         <  1:363705/139‑1 (MQ=255)
                                                                                       tAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCg                                        >  2:2685/1‑139 (MQ=255)
                                                                                             tCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGtggctt                                  >  1:725520/1‑139 (MQ=255)
                                                                                             tCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGtggctt                                  >  1:699943/1‑139 (MQ=255)
                                                                                               ccACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCa                                                                                   <  1:1074551/88‑1 (MQ=255)
                                                                                               ccACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCa                                                                                   >  2:1074551/1‑88 (MQ=255)
                                                                                               ccACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGtggcttgg                                <  2:699943/139‑1 (MQ=255)
                                                                                                  ccGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTggcttggctt                             <  1:657949/139‑1 (MQ=255)
                                                                                                    gTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCtttt                           <  2:1162984/139‑1 (MQ=255)
                                                                                                    gTTGGCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCtttt                           <  2:390716/139‑1 (MQ=255)
                                                                                                       ggCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTa                        >  2:378034/1‑139 (MQ=255)
                                                                                                       ggCGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTa                        >  2:725457/1‑139 (MQ=255)
                                                                                                         cGGTGTTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTAcc                      <  2:32380/139‑1 (MQ=255)
                                                                                                             gtTATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTggg                  >  1:1056831/1‑139 (MQ=255)
                                                                                                              ttATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTggg                  <  1:930924/138‑1 (MQ=255)
                                                                                                              ttATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTggg                  <  1:1049050/138‑1 (MQ=255)
                                                                                                              ttATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGc                 <  1:200842/139‑1 (MQ=255)
                                                                                                              ttATCAGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGc                 <  1:2685/139‑1 (MQ=255)
                                                                                                                   aGCAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGCCtgtg            <  2:1190152/139‑1 (MQ=255)
                                                                                                                     cAAAAGCCTCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGGTTCTTGGCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGCCTGTgcg          <  2:1182489/139‑1 (MQ=255)
                                                                                                                     cAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGCCTGTgcg          <  1:795765/139‑1 (MQ=255)
                                                                                                                        aaGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGt                                                               <  2:834501/83‑1 (MQ=255)
                                                                                                                        aaGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGt                                                               >  1:834501/1‑83 (MQ=255)
                                                                                                                             tCGGGCGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGCCTGTGCGCACTTCAc  >  1:1163366/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TACCATTGCCATGTGTCTGGCGCTGACGGTTGCCTATTATTTAAATCTGGATGAACCCTATTGGGCGATGACCTCGGCTGCAGTGGTTAGCTTTCCCACCGTTGGCGGTGTTATCAGCAAAAGCCTCGGACGCATCGCTGGCAGTTTGCTCGGAGCCATTGCGGCACTGCTTCTTGCCGGGCATACGCTCAATGAGCCGTGGTTTTTTCTATTGAGCATGTCGGCGTGGCTTGGCTTTTGTACCTGGGCCTGTGCGCACTTCAC  >  NZ_CP009273/1716458‑1716721

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑