Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,659,0620CA55.6% 4.3 / 46.9 36A37S (GCG→TCG) gadXacid resistance transcriptional activator GadX
Reads supporting (aligned to +/- strand):  ref base C (10/6);  new base A (0/20);  total (10/26)
Fisher's exact test for biased strand distribution p-value = 3.15e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAT  >  NZ_CP009273/3658974‑3659143
                                                                                         |                                                                                 
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)
                                                                            tGGCTTGCACCCGAAAAACCCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc       <  2:235516/90‑1 (MQ=255)
                                                                            tGGCTTCCACCCAAAAAACCCGGTTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc       <  2:345461/90‑1 (MQ=255)
                                                                                 tGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAt  <  2:47256/90‑1 (MQ=255)
                                                                                         |                                                                                 
AGCATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAT  >  NZ_CP009273/3658974‑3659143

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