Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA40.9% 47.8 / 22.4 44L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (19/7);  new base A (0/18);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 7.53e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAG  >  NZ_CP009273/4058870‑4059039
                                                                                         |                                                                                 
aaTTTGCCTTTAGGTCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                  <  2:1053388/90‑1 (MQ=255)
   ttGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAg                                                                               >  2:518866/1‑90 (MQ=255)
   ttGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAg                                                                               >  2:1186986/1‑90 (MQ=255)
     gCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGAt                                                                             >  1:689377/1‑90 (MQ=255)
         ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                         >  2:1098647/1‑90 (MQ=255)
          tAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATcc                                                                        >  2:1269126/1‑90 (MQ=255)
             cccaaaTATTTTCCCGCGAGTCAAAAAAAAATCCGGTGAGCATATCGGAAAACGCGAAAAAAACCTCCCCCACAAGAAAAATAACCCCGc                                                                     <  1:1428199/86‑1 (MQ=255)
               ccAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCa                                                                   >  1:41527/1‑90 (MQ=255)
                cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                  <  1:774791/74‑1 (MQ=255)
                cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                  >  2:774791/1‑74 (MQ=255)
                      ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg                                                            >  1:1319978/1‑90 (MQ=255)
                      ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg                                                            >  2:1428199/1‑90 (MQ=255)
                       tttCCGCGAGTCGAAAAAAAATCCGGTGAGCATATCGGAAACCGCGGAAAAAAACTTCGCCACCAGAAAGATAACCCCGCCATAGTAAGc                                                           <  2:1318141/90‑1 (MQ=255)
                          ccGCGAGTCGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGAAAAAACCTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAgg                                                        <  2:1006727/90‑1 (MQ=255)
                           cgcgAGTCAAGAAAAAACCCGGGGACAATACCGGAAACCGCGGAAAAAACCTCCCCCCCAAAAAAAAAAATCCCGCCATAGTAAGCAGGc                                                       <  2:973003/90‑1 (MQ=255)
                            gggagtcgaaaAAAAACCCGGAGAAAATACCGAAAAACGCGAAAAAAACCCCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa                                                      <  1:1186986/88‑1 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                 <  2:1334663/63‑1 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                 >  1:1334663/1‑63 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                      >  1:1082219/1‑90 (MQ=255)
                            gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa                                                      >  1:912333/1‑90 (MQ=255)
                            gcgAGTCGAGAAAAAACCCGGGGAGCATACCGGAAACCGGGAAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCa                                                      <  1:76220/90‑1 (MQ=255)
                            gcgAGTCGAAAAAAAACCCGGTGAGCATACCGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa                                                      <  1:518866/90‑1 (MQ=255)
                            gcgAGTAGAGAAAAAACCGGGGGAGAATACCGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCa                                                      <  1:1194779/90‑1 (MQ=255)
                               agtcgagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATcc                                                   >  1:796740/1‑90 (MQ=255)
                                   gagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCga                                                                     >  1:772262/1‑67 (MQ=255)
                                   gagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCga                                                                     <  2:772262/68‑2 (MQ=255)
                                   gagtAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTaa                                               >  1:1111243/1‑90 (MQ=255)
                                        aaaaCCGGGTGAGCATATGGAAAAACGGGAAAAAAATCTCCCCCACAAAAAAGATATTCCCGCCATAGTAAGCAGGCATCCCTAATTCAt                                          <  2:1082219/90‑1 (MQ=255)
                                                     cATATCGGTAAACGCGAAAAAAAACTTCGCCACAAAAAAGAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                             <  2:1278851/90‑1 (MQ=255)
                                                         tCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                         >  2:1193181/1‑90 (MQ=255)
                                                          cGGAAAACGCGGAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                        <  1:174640/90‑1 (MQ=255)
                                                          cGGAAAACGCGAAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                        <  1:169821/90‑1 (MQ=255)
                                                           ggTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGACATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgaggag                       >  1:1044661/1‑90 (MQ=255)
                                                            gaaaaCGCGAAAAAAAACTTCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                      <  1:1349084/88‑1 (MQ=255)
                                                             aaaaCGCGAAAAAAACCTCCCCCCCCAGAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGAt                     <  2:355528/89‑1 (MQ=255)
                                                                  gcgTTAAAAAATTTCCCCACAAAAAAGATACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga                <  2:1319978/90‑1 (MQ=255)
                                                                  gcgGAAAAAAACTCCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga                <  2:1075358/90‑1 (MQ=255)
                                                                    gaaaaaaacCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt              <  2:2234/87‑1 (MQ=255)
                                                                    gaaaaaaaaCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt              <  2:615374/87‑1 (MQ=255)
                                                                    aaaaaaaacCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt              <  1:432210/87‑1 (MQ=255)
                                                                     gaaaaaaaCCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt              <  1:1331312/88‑1 (MQ=255)
                                                                     gTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTa             >  1:948871/1‑90 (MQ=255)
                                                                      tAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAc            >  1:69269/1‑90 (MQ=255)
                                                                                ccccaccaAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCagag  <  2:1212597/88‑1 (MQ=255)
                                                                                         |                                                                                 
AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAG  >  NZ_CP009273/4058870‑4059039

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