Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,659,0580TC52.4% 90.0 / 126.6 86D38G (GAT→GGT) gadXacid resistance transcriptional activator GadX
Reads supporting (aligned to +/- strand):  ref base T (30/10);  new base C (0/44);  total (30/56)
Fisher's exact test for biased strand distribution p-value = 1.57e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-16
Rejected as consensus: Frequency below/above cutoff threshold.

CTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATT  >  NZ_CP009273/3658970‑3659144
                                                                                          |                                                                                       
ctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATc                                                                                        <  1:1604704/90‑1 (MQ=255)
     gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa                                                                                   <  2:304614/90‑1 (MQ=255)
     gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa                                                                                   <  1:2066262/90‑1 (MQ=255)
       aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                                 >  2:1473598/1‑90 (MQ=255)
       aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                                 >  2:1021748/1‑90 (MQ=255)
        tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                >  2:1020169/1‑90 (MQ=255)
        tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                >  2:1004229/1‑90 (MQ=255)
        tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                >  1:1529748/1‑90 (MQ=255)
        tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                >  2:1872097/1‑90 (MQ=255)
        tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                                >  1:1620867/1‑90 (MQ=255)
         cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                               >  1:1769771/1‑90 (MQ=255)
         cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                               >  1:1731003/1‑90 (MQ=255)
         cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                               >  1:53869/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:1953694/1‑90 (MQ=255)
                            gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                           >  1:2253825/1‑90 (MQ=255)
                            gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                           >  1:510250/1‑90 (MQ=255)
                            gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                           >  1:791367/1‑90 (MQ=255)
                             tttCCCTTAACACGAATAAAAATTTTCCAACCCTTTTATCTACTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAAt                                                          <  1:2206401/88‑1 (MQ=255)
                             tgccccTTGACACAAATAAAAATTTTTCAACCCACTTATCTCCTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAAt                                                          <  1:210362/87‑1 (MQ=255)
                              gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa                                                         >  2:945032/1‑90 (MQ=255)
                              gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa                                                         >  2:61847/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                                             >  2:600110/1‑90 (MQ=255)
                                          gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                             >  1:168871/1‑90 (MQ=255)
                                          gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                             >  1:450038/1‑90 (MQ=255)
                                           aaTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATc                                                                                        >  1:1547456/1‑47 (MQ=255)
                                           aaTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATc                                                                                        <  2:1547456/47‑1 (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:1941960/90‑1 (MQ=255)
                                                                                tGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc        <  1:867935/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)
                                                                                          |                                                                                       
CTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATT  >  NZ_CP009273/3658970‑3659144

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: