Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,659,0460AT60.0% 13.4 / 43.7 35I42N (ATT→AAT) gadXacid resistance transcriptional activator GadX
Reads supporting (aligned to +/- strand):  ref base A (11/3);  new base T (0/21);  total (11/24)
Fisher's exact test for biased strand distribution p-value = 8.72e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-07
Rejected as consensus: Frequency below/above cutoff threshold.

CTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTT  >  NZ_CP009273/3658964‑3659132
                                                                                   |                                                                                      
ctcctcctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCt                                                                                 >  1:490629/1‑90 (MQ=255)
     cctTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCAt                                                                            <  2:214507/90‑1 (MQ=255)
          aGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCa                                                                       >  2:565687/1‑90 (MQ=255)
           gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa                                                                      >  2:351940/1‑90 (MQ=255)
           gCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaa                                                                      <  1:565687/90‑1 (MQ=255)
             aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                    >  2:389317/1‑90 (MQ=255)
              tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                   >  2:240057/1‑90 (MQ=255)
               cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                  >  1:700279/1‑90 (MQ=255)
               cGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAAc                                                                  >  1:147006/1‑90 (MQ=255)
                         aatGAAAGCGTGTCCCTTGACACGAATAAAAAATTTTCAACACACTTATCTCCTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCg                                                        <  2:195961/90‑1 (MQ=255)
                                   tctCCCTTGACACAAATAAAAAATTTTCAACCCCCTTATCTACTCGATTTTGGCTTGCATCCGAAAAACCCAGGTCACCGCCATTAAAAt                                              <  1:389317/88‑1 (MQ=255)
                                      cccTTGACCCAAAAAAAAAATTTAAAACACCTTTACCTCCGCAATTTTGGCTTCCTCCCGCAAAACCCAGGGCACCGCCATTAAAATAGc                                           <  2:410952/90‑1 (MQ=255)
                                      cccTTGACACAAATAAAAATTTTTCAACACACTTATCTCCTCGATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGc                                           <  2:483988/90‑1 (MQ=255)
                                      cccTTGACACAAAAAAAAATTTTTCAACCCATTTATCTCCTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGc                                           <  2:653237/90‑1 (MQ=255)
                                        cTTGACAAAAATAAAAATTTTCCAACACTCTTTTCTACCCGATTTTGGCTTCCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGa                                         <  2:599291/90‑1 (MQ=255)
                                            acacAAAAACAAATTTTTCAACCCATTTATCTCCTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATAtt                                     <  2:390476/90‑1 (MQ=255)
                                             cacaAATAAAAATTTTTAAACACACTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                    <  2:147006/90‑1 (MQ=255)
                                               caaaTAAAAATTTTTCAACACCCTTATCTCCTCGATTTTGGCTTGCACCCGCAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCAc                                  <  2:67993/88‑1 (MQ=255)
                                               cGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAc                                  >  2:134017/1‑90 (MQ=255)
                                                 aaTACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCa                                >  2:574409/1‑90 (MQ=255)
                                                   taaaaatTTTTCAACACTTTTTTCTACTCAATTTTGGCTTCCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                              <  2:804402/87‑1 (MQ=255)
                                                    taaaaaTTTTCCACCCCACTTTTCTCCTCAATTTTGGCTTCCCCCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCAtt                              <  1:134017/88‑1 (MQ=255)
                                                    taaaaaTTTTCCAACCCACTTTTCTACTCGATTTTGGCTTCCACCCCAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCAtt                              <  2:723336/88‑1 (MQ=255)
                                                             tCAACCCTCTTTCCTCCTCGATTTTGGCTTGCACCCAAAAAACCCGGCTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                    <  2:490629/90‑1 (MQ=255)
                                                             tCAACCCACTTATCTACTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                    <  2:342034/90‑1 (MQ=255)
                                                             tCAACACCTTTTTCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                    <  1:255854/90‑1 (MQ=255)
                                                             tCAACACATTTATCTACTAAATTTTGGCTTGCACCCGAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                    <  1:41912/90‑1 (MQ=255)
                                                             tCAACACACTTATCTCCTCGATTTTGGCTTGCACCCGAAAAACCCGGGCCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                    <  2:713078/90‑1 (MQ=255)
                                                                   ccatttATCTACTAAATTTTGGTTTCCACCAGAAAAACCGGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat               <  1:491659/86‑1 (MQ=255)
                                                                   acatTTATCTACTCAATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat               <  2:172944/90‑1 (MQ=255)
                                                                   acacTTATCTACTCAATTTTGGCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat               <  1:118433/90‑1 (MQ=255)
                                                                    cctttATCTACTCAATTTTGGCTTGCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAata              <  2:585642/87‑1 (MQ=255)
                                                                          tctcCTCAATTTTGGCTTCCACCCGAAAAACCCAGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTa        <  2:348337/90‑1 (MQ=255)
                                                                              cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc    >  1:642079/1‑90 (MQ=255)
                                                                                cGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCtt  >  2:333640/1‑90 (MQ=255)
                                                                                   |                                                                                      
CTCCTCCTTGAGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTT  >  NZ_CP009273/3658964‑3659132

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