Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,659,0620CA56.3% 64.6 / 158.7 87A37S (GCG→TCG) gadXacid resistance transcriptional activator GadX
Reads supporting (aligned to +/- strand):  ref base C (31/7);  new base A (0/49);  total (31/56)
Fisher's exact test for biased strand distribution p-value = 3.50e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-15
Rejected as consensus: Frequency below/above cutoff threshold.

GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAT  >  NZ_CP009273/3658975‑3659151
                                                                                         |                                                                                          
gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa                                                                                          <  1:2066262/90‑1 (MQ=255)
gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa                                                                                          <  2:304614/90‑1 (MQ=255)
  aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                                        >  2:1021748/1‑90 (MQ=255)
  aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                                        >  2:1473598/1‑90 (MQ=255)
   tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                       >  2:1020169/1‑90 (MQ=255)
   tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                       >  2:1004229/1‑90 (MQ=255)
   tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                       >  2:1872097/1‑90 (MQ=255)
   tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                       >  1:1620867/1‑90 (MQ=255)
   tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                       >  1:1529748/1‑90 (MQ=255)
    cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                                      >  1:53869/1‑90 (MQ=255)
    cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                                      >  1:1731003/1‑90 (MQ=255)
    cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                                      >  1:1769771/1‑90 (MQ=255)
     gggAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAcc                                                                                     <  1:1743942/90‑1 (MQ=255)
      ggAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCa                                                                                   <  2:1235274/90‑1 (MQ=255)
              aatGAAAGCGTTTCCCTTGCCAAAAATAAAAAATTTTCAACACCCTTATCTACCCAAATTTGGCTTCCACCCGAAAAAACCGGGTCCCCg                                                                           <  2:1769771/90‑1 (MQ=255)
              aatAAAAGCGTGCCCCTTAACACAAAAAAAAAATTTTCAACACACTTATCTCCTCAATTTTGGCTTGCACCCGAAAAAACCAGGTCACCg                                                                           <  2:749602/90‑1 (MQ=255)
                 gAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCa                                                                        >  2:2241980/1‑90 (MQ=255)
                   aaGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCAtt                                                                      >  2:1104870/1‑90 (MQ=255)
                   aaGCGTGCCCCTTGACACAAATAAAAAATTTTCAACCCACTTATCTACTCAAATTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCAtt                                                                      <  2:53869/90‑1 (MQ=255)
                     gCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaa                                                                    >  2:2206401/1‑90 (MQ=255)
                      cGTGTCCTTTCACACGAATAAAAATTTTCCAACACCCTTATCAACAGGATTTGGGCTTCCCCCCAAAAAAACCAGGTCACCGCCATTaaa                                                                   <  2:2258203/90‑1 (MQ=255)
                       gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                                  >  1:510250/1‑90 (MQ=255)
                       gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                                  >  1:2253825/1‑90 (MQ=255)
                       gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                                  >  1:791367/1‑90 (MQ=255)
                       gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                                  >  1:1953694/1‑90 (MQ=255)
                        tttCCCTTAACACGAATAAAAATTTTCCAACCCTTTTATCTACTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAAt                                                                 <  1:2206401/88‑1 (MQ=255)
                        tgccccTTGACACAAATAAAAATTTTTCAACCCACTTATCTCCTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAAt                                                                 <  1:210362/87‑1 (MQ=255)
                         gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa                                                                >  2:61847/1‑90 (MQ=255)
                         gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa                                                                >  2:945032/1‑90 (MQ=255)
                           cccTTGACACAAAAAAAAAATTTTCAACCCCCTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGc                                                              <  2:183754/90‑1 (MQ=255)
                           cccTTGACACAAAAAAAAAATTTACAACCCCTTTATCTACTCGAATTTGGCTTGCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGc                                                              <  2:510250/90‑1 (MQ=255)
                             cTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCAATTAAAATAGCGa                                                            <  1:1645082/90‑1 (MQ=255)
                             cTTGACACGAATAAAAAATTTTCAACCCACTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGa                                                            <  1:111983/90‑1 (MQ=255)
                             cTTGACAAAAATAAAAAATTTTCAACCCTTTTATCTCCTCAATTTTGGCTTGCACCCCAAAAACCCAGGTCACCGCCATTAAAATAGCGa                                                            <  1:1872097/90‑1 (MQ=255)
                             cTTGACAAAAAAAAAAATTTTTCAACCCCTTTATCTCCTCAATTTTGGCTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGa                                                            <  2:450038/90‑1 (MQ=255)
                                  cacaAATAAAAATTTTTCAACACATTTTTCTCCTCAATTTTGGCTTGCATCCCAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  1:1185767/90‑1 (MQ=255)
                                  cacaAATAAAAATTTTTCAACACATTTATCTACCCGATTTTGGCTTGCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  1:1314401/90‑1 (MQ=255)
                                  cacaAATAAAAATTTTCCAACACTTTTTTCTACTCAATTTTGGCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  2:942134/90‑1 (MQ=255)
                                  cacaAATAAAAAATTTTCAACCCCCTTATCTACTCGATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  1:61847/90‑1 (MQ=255)
                                  cacaAATAAAAAATTTCCAACACACTTATCTACTCAAATTTGGCTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  1:2119308/90‑1 (MQ=255)
                                  cacaAATAAAAAATTTACAACACACTTATCTACTCGAATTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  1:475992/90‑1 (MQ=255)
                                     gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                                    >  1:450038/1‑90 (MQ=255)
                                     gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                                    >  2:600110/1‑90 (MQ=255)
                                     gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                                    >  1:168871/1‑90 (MQ=255)
                                      aaTAAAAAATTTACAACACACTTATCTCCCCGAATTTGGCTTGCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCa                                                   <  1:1072982/90‑1 (MQ=255)
                                       aaaaaaaaTTTTCAAACACCTTTATCTCCCAAATTTGGGCTTCCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCa                                                   <  2:696151/86‑1 (MQ=255)
                                       aTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCAt                                                  >  1:1313265/1‑90 (MQ=255)
                                        aaaaaatTTTTCAACCCCCTTATCTACCCGTTTTTGGCTTGCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                 <  1:2194557/87‑1 (MQ=255)
                                        aaaaaatTTTTCAACACACTTATCTACTCAATTTTGGCTTGCATCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                 <  1:1104870/87‑1 (MQ=255)
                                        aaCAAAATTTCCACCCACCTTATATCCTCAAATTGGCCTTGCCCCCGAAAAAACCAGGGCACCGCCATTAAAATAGCGATATTCACCAtt                                                 <  1:2211357/89‑1 (MQ=255)
                                         aaaaatTTTTCACCCCCCTTATCTCCTCAATTTTGGTTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTa                                                <  1:1021748/88‑1 (MQ=255)
                                            aaaTTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                             >  2:1758430/1‑90 (MQ=255)
                                             atttttCAACACCTTTATCTCCTAAATTTTGGCTTGCACCCGAAAAAACCGGGACCCCGCCATTAAAATAGCGATATTCACCATTAACCa                                            <  2:2147878/88‑1 (MQ=255)
                                             attttcCAACACATTTATCTACTCAATTTTGGCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCa                                            <  1:2091576/88‑1 (MQ=255)
                                              aTTTCCAACACATTTTTCTCCTCAAATTTGGCTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCAt                                           <  1:1689437/90‑1 (MQ=255)
                                               ttttCAACCCTTTTTTCTCCTAAATTTTGGCTTCCACCCAAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATg                                          <  2:712288/90‑1 (MQ=255)
                                                  tCAACACTTTTATCTCCTCAATTTTGGTTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                       <  1:1473598/90‑1 (MQ=255)
                                                  tCAACACACTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                       <  2:1622417/90‑1 (MQ=255)
                                                  ccAACCCTCTTATCTACTCAATTTTGGCTTCCACCCGAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                       <  1:1020169/89‑1 (MQ=255)
                                                  acAACACATTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                       <  2:831115/89‑1 (MQ=255)
                                                   caaaCACCTTTATCTACAAGAATTTGGTTTGCATCCAAAAAACCCAGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                       <  1:2209104/88‑1 (MQ=255)
                                                     aaCACTTTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgag                                     <  2:2183918/90‑1 (MQ=255)
                                                       cccacTTATCTCCTCAATTTGGGTTTCCACCCAAAAAACCAAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                   <  2:174764/88‑1 (MQ=255)
                                                       cacccTTATCTCCTCGATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                   <  1:600110/90‑1 (MQ=255)
                                                       cacatTTATCTCCTCGATTTTGGTTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                   <  1:196756/90‑1 (MQ=255)
                                                       cacatTTATCTACTCTATTTTGTCTTGCACCCACAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                   <  1:2180028/90‑1 (MQ=255)
                                                       cacatTTATCTACTCAATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                   <  1:557091/90‑1 (MQ=255)
                                                        ccctttATCTCCCCAATTTTGGTTTGCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                  <  2:2101397/86‑1 (MQ=255)
                                                        ccctttATCTACTAAAATTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCCTGGTGAGAat                                  <  2:2010101/86‑1 (MQ=255)
                                                        ccccTTACCTACTCGAATTGGGCTCGCCCCCAAAAAACCGAGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                  <  2:1211634/87‑1 (MQ=255)
                                                        acctttATCTCCTCAATTTTGGTTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                  <  1:945032/86‑1 (MQ=255)
                                                        acatTTATCTACCCAATTTTGGCTTGCAACCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                  <  1:1528881/90‑1 (MQ=255)
                                                              aTCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAttt                            >  2:440448/1‑90 (MQ=255)
                                                                cTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTAt                          >  2:358837/1‑90 (MQ=255)
                                                                     cGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt                     >  1:824956/1‑90 (MQ=255)
                                                                           tGGCTTGCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc               <  2:1368489/90‑1 (MQ=255)
                                                                           tGGCTTGCACCCAAAAAAACCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc               <  1:732089/90‑1 (MQ=255)
                                                                           tGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc               <  1:867935/90‑1 (MQ=255)
                                                                           tGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc               <  1:1941960/90‑1 (MQ=255)
                                                                           tGGCTTACACCCAAAAAACCCAGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc               <  1:2287735/90‑1 (MQ=255)
                                                                            tgctttGCACCCAAAAAACCCAGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc               <  1:2006045/86‑1 (MQ=255)
                                                                                  gCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAtt         >  2:128622/1‑90 (MQ=255)
                                                                                  gCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAtt         >  2:1627073/1‑90 (MQ=255)
                                                                                      ccGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                       >  2:2030572/1‑40 (MQ=255)
                                                                                      ccGAAAAACCAAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc     <  1:1997793/90‑1 (MQ=255)
                                                                                      ccGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                       <  1:2030572/40‑1 (MQ=255)
                                                                                        gCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa   >  2:1903568/1‑90 (MQ=255)
                                                                                         cAAACACCAGTGCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAt  <  1:18414/90‑1 (MQ=255)
                                                                                         |                                                                                          
GCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAT  >  NZ_CP009273/3658975‑3659151

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: