Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,058,9580GA48.6% 29.3 / 19.9 35L57L (CTG→TTG) yihOMFS transporter
Reads supporting (aligned to +/- strand):  ref base G (14/4);  new base A (0/17);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058869‑4059038
                                                                                          |                                                                                
aaaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAg                                                                                  <  1:788497/90‑1 (MQ=255)
 aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                 <  1:831476/90‑1 (MQ=255)
 aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa                                                                                 <  1:337161/90‑1 (MQ=255)
          ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                >  2:257472/1‑82 (MQ=255)
          ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa                                                                                <  1:257472/82‑1 (MQ=255)
          ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc                                                                        >  1:1087120/1‑90 (MQ=255)
            aGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATccc                                                                      >  2:763449/1‑90 (MQ=255)
                        tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc                                                          >  1:718687/1‑90 (MQ=255)
                        ttcccGCGAGTCGAGAAAAAATCCGGTGAGAATATCGGAAAACGCGGAAAAAAACTTCCCCACAAAAAAGATAACCCCGCCATAGTAAGc                                                          <  1:641110/87‑1 (MQ=255)
                            cgcgAGTCGAGAAAAAACCCGGTGAGAATATGGGAAACCGTGGAAAAAAACTCCCCCCCCAAAAAAATAACCCCGCCATAGTAAGCAGGc                                                      <  1:1078809/90‑1 (MQ=255)
                             gcgAGTCGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTCCGCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCa                                                     <  1:763449/90‑1 (MQ=255)
                             gcgAGTCGAGAAAAAACCGGGTGAGAATATCGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa                                                     <  2:186949/90‑1 (MQ=255)
                             gagagtcgagaAAAAACCGGGTGAGCATACCGGAAAACGCGAAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCa                                                     <  2:921814/88‑1 (MQ=255)
                                          aaaTCCGGTGAGCATATGGGTAAACGGGGTAAAAAACTTCCCCCCCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                        <  1:173896/90‑1 (MQ=255)
                                          aaaTCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc                                        <  1:1149494/90‑1 (MQ=255)
                                                ggTGAGCATATGGGAAAACGCGGAAAAAAACTCCCCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                  <  1:820592/90‑1 (MQ=255)
                                                ggTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata                                  >  1:67885/1‑90 (MQ=255)
                                                     gCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAtt                             >  2:23847/1‑90 (MQ=255)
                                                      cATATCGAAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                            <  2:718687/90‑1 (MQ=255)
                                                      aatatCGGAAAAGGGGGAAAAAACCTCCCCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                            <  2:1058533/89‑1 (MQ=255)
                                                      aatatCGGAAAACGGGAAAAAAAACTCCGCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt                            <  1:764676/89‑1 (MQ=255)
                                                       atatCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                           >  2:691192/1‑90 (MQ=255)
                                                          tCGGAAACCGCGAAAAAAAACTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg                        <  1:1183216/90‑1 (MQ=255)
                                                           cGGAAACCGCGAAAAAAACCTCCGCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga                       <  1:791014/90‑1 (MQ=255)
                                                             gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     >  2:576333/1‑90 (MQ=255)
                                                             gTAAACGCGGTAAAAAACTTCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga                     <  1:23847/90‑1 (MQ=255)
                                                                aaCGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                  >  1:377816/1‑90 (MQ=255)
                                                                 accgggTAAAAAATCTCCCCCACGAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc                  <  1:1086587/86‑1 (MQ=255)
                                                                   gcgGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga               >  2:195457/1‑90 (MQ=255)
                                                                     gaaaaaaacCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt             <  2:67885/87‑1 (MQ=255)
                                                                       aaaaaaaCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAc           <  2:1138510/89‑1 (MQ=255)
                                                                                 tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  2:397464/1‑90 (MQ=255)
                                                                                 tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  2:612738/1‑90 (MQ=255)
                                                                                 tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  1:1174374/1‑90 (MQ=255)
                                                                                 tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  >  1:731075/1‑90 (MQ=255)
                                                                                 cccccacaaAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga  <  1:1003850/87‑1 (MQ=255)
                                                                                          |                                                                                
AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA  >  NZ_CP009273/4058869‑4059038

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: