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

ATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATT  >  NZ_CP009273/3658977‑3659144
                                                                                      |                                                                                  
aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                                >  2:1236003/1‑90 (MQ=255)
aTCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaa                                                                                >  1:1280071/1‑90 (MQ=255)
 tCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCaaaaa                                                                               >  2:649027/1‑90 (MQ=255)
       gaAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGt                                                                         >  2:682764/1‑90 (MQ=255)
            aatAAAAGCGTGCCCCTTGACACAAATAAAAAATTTTCAACACACTTATCTACTCGATTTTGGCTTCCACCCGCAAAACCCAGGTCACCg                                                                    <  1:1248765/90‑1 (MQ=255)
                 aaGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCAtt                                                               >  2:85271/1‑90 (MQ=255)
                 aaGCGTGTCCCTTGACACAAATAAAAATTTTTCAACACTTTTATCTACTCGAATTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCAtt                                                               <  1:738676/90‑1 (MQ=255)
                    cGTGCCCCTTGACACAAATAAAAAATTTTCAACCCACTTATCTCCTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaa                                                            <  1:1425/90‑1 (MQ=255)
                     gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                           >  2:1248765/1‑90 (MQ=255)
                     gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                           >  2:1321473/1‑90 (MQ=255)
                     gtgtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTaaaa                                                           >  2:183163/1‑90 (MQ=255)
                       gtCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATa                                                         >  2:608178/1‑90 (MQ=255)
                         cccTTGACACAAATAAAAATTTTTCAACACACTTATCTCCTCGAATTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGc                                                       <  1:1204198/90‑1 (MQ=255)
                         cccTTGACACAAATAAAAAATTTCCAACACACTTATCTACTCAATTTTGGCTTCCACCCAAAAAAACCAGGCCACCGCCATTAAAATAGc                                                       <  2:1280071/90‑1 (MQ=255)
                           cTTGACACAAATAAAAAATTTTCAACCCACTTATCTCCTCGATTTTGGCTTGCACCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGa                                                     <  1:57017/90‑1 (MQ=255)
                           cTTGACACAAATAAAAAATTTTAAACACACTTATCTCCTCGATTTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGa                                                     <  1:1236003/90‑1 (MQ=255)
                           cTTGACAAAAAAAAAAATTTTTCAACACCCTTATCTACTCAAATTTGGTTTGCACCCCAAAAACCCAGGCCACCGCCATTAAAATAGCGa                                                     <  1:1321473/90‑1 (MQ=255)
                           cTTCACACAAAAAAAAAATTTCCAACCCATTTATCTACTCGATTTTGGTTTGCACCCAAAAAAACCAGGTCCCCGCCATTAAAATAGCGa                                                     <  2:664199/90‑1 (MQ=255)
                                cacaAATAAAAATTTTTCAACACTCTTATCTACTCAATTTTGGCTTGCACCCGAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTc                                                <  2:481272/90‑1 (MQ=255)
                                cacaAATAAAAAATTTTCAACACCCTTATCTACTCGATTTTGGCTTGCATCCGAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTc                                                <  1:201849/90‑1 (MQ=255)
                                cacaAATAAAAAAATTTAAACCCACTTATCTACTCGAATTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTc                                                <  1:525569/90‑1 (MQ=255)
                                  caaaTAAAAAATTTTCAACACCCTTATCTACTCGATTTTGGCTTGCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAc                                              <  1:85271/88‑1 (MQ=255)
                                   gAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCAcc                                             >  2:933790/1‑90 (MQ=255)
                                      aaaaaaaTTTCCAACACATTTACCTACACGAATTTGGCTTCCCCCCAAAAAACCCGGCCCCCCGCCATTAAAATAGCGATATTCACCAtt                                          <  2:1276995/87‑1 (MQ=255)
                                       aaaaaaTTTTCCAACCCTTTTTTCTTCTCAATTTTGGCTTCCACCCAAAAAAACCGGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                          <  2:1305219/88‑1 (MQ=255)
                                       aCAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTa                                         >  2:179293/1‑90 (MQ=255)
                                                tCAACACATTTATCTCCTCGATTTTGGCTTGCACCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                <  1:1333076/90‑1 (MQ=255)
                                                 taaccactcTTATCTACTCAATTTTGGTTTGCATCCAAAAAACCCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                <  2:1188120/86‑1 (MQ=255)
                                                   aaCCCATTTATCTACTCGATTTTGGTTTCCACCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgag                              <  1:1319133/90‑1 (MQ=255)
                                                   aaCACTCTTATCTCCTCGAATTTGGCTTGCATCCCAAAAACCCAGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgag                              <  1:1048407/90‑1 (MQ=255)
                                                   aaCAAATTTATCCCCTCGATTTTGGCTTGCAACCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgag                              <  1:183163/90‑1 (MQ=255)
                                                              cTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTAt                   >  2:1262819/1‑90 (MQ=255)
                                                                 cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc                >  2:468319/1‑90 (MQ=255)
                                                                 cTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTc                >  1:590459/1‑90 (MQ=255)
                                                                     aaTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                >  2:242565/1‑70 (MQ=255)
                                                                     aaTTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGc            >  1:208548/1‑90 (MQ=255)
                                                                     aaTTTGGCTTGCATCCGAAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                <  1:242565/70‑1 (MQ=255)
                                                                         tGGCTTCCACCCAAAAAACCCGGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATAc        <  1:768484/90‑1 (MQ=255)
                                                                               gCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAAtt  >  2:563351/1‑90 (MQ=255)
                                                                                      |                                                                                  
ATCGGGAGAAATAATGAAAGCGTGTCCCTTGACACGAATACAAAATTTTCAACACACTTATCTACTCGAATTTGGCTTGCATCCGCAAAAACCAGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATT  >  NZ_CP009273/3658977‑3659144

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