Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,8990CT57.1% 22.9 / 52.0 42A26V (GCC→GTC) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base T (12/12);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGA  >  NZ_CP009273/49763‑50031
                                                                                                                                        |                                                                                                                                    
gttgacgacGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCtgcctacc                                                                                                                                     >  2:487579/4‑138 (MQ=255)
      gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCg                                                                                                                              >  1:12108/1‑139 (MQ=255)
          gTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGccgg                                                                                                                          >  2:448220/1‑137 (MQ=255)
             tACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCACGCCGTGGAACCTGCCTACCGATCTCGCCGGGt                                                                                                                        <  1:487579/138‑1 (MQ=255)
                    tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACg                                                                                                                 >  2:277436/1‑138 (MQ=255)
                    tACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACg                                                                                                                 >  2:283643/1‑138 (MQ=255)
                     aCGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGGTTAAACGc                                                                                                                >  2:429330/1‑138 (MQ=255)
                      cGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  1:193380/138‑1 (MQ=255)
                            gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGccg                                                                                                                           <  1:338208/120‑2 (MQ=255)
                            gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGccg                                                                                                                           >  2:338208/1‑119 (MQ=255)
                            gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCa                                                                                                               >  1:95366/1‑132 (MQ=255)
                            gTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCa                                                                                                               <  2:95366/132‑1 (MQ=255)
                                cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaa                                                                                                     >  2:480338/1‑138 (MQ=255)
                                cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaa                                                                                                     >  1:500690/1‑138 (MQ=255)
                                 gACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTtaaata                                                                                                   <  1:448220/139‑1 (MQ=255)
                                     aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                                                                                               <  1:230513/139‑1 (MQ=255)
                                     aTTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                                                                                               <  1:92008/139‑1 (MQ=255)
                                        ttttttATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                                                                                             <  2:12108/138‑1 (MQ=255)
                                                 gggAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                                                   >  1:217980/1‑139 (MQ=255)
                                                                  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGt                                                                                                                        >  1:451177/1‑85 (MQ=255)
                                                                  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGt                                                                                                                        <  2:451177/85‑1 (MQ=255)
                                                                  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaataaa                                                                                                 >  1:510153/1‑108 (MQ=255)
                                                                  aGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTtaaataaa                                                                                                 <  2:510153/108‑1 (MQ=255)
                                                                            cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                                                               >  2:294293/1‑132 (MQ=255)
                                                                            cggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaa                                                               <  1:294293/132‑1 (MQ=255)
                                                                              gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc                                                      >  2:97591/1‑139 (MQ=255)
                                                                              gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                                                                                             >  1:397553/1‑100 (MQ=255)
                                                                              gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCg                                                                                             <  2:397553/100‑1 (MQ=255)
                                                                                           aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCaa                                                                                                              <  1:474520/70‑1 (MQ=255)
                                                                                           aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCaa                                                                                                              >  2:474520/1‑70 (MQ=255)
                                                                                              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      <  1:277285/139‑1 (MQ=255)
                                                                                               cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      >  2:546583/1‑138 (MQ=255)
                                                                                               cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      >  1:262089/1‑138 (MQ=255)
                                                                                                  tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      >  2:32824/1‑135 (MQ=255)
                                                                                                  tATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      <  1:32824/135‑1 (MQ=255)
                                                                                                    tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAata                                 <  1:97591/138‑1 (MQ=255)
                                                                                                    tCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAata                                 <  1:283643/138‑1 (MQ=255)
                                                                                                                aCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCt                          <  2:390611/133‑1 (MQ=255)
                                                                                                                aCGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCt                          >  1:390611/1‑133 (MQ=255)
                                                                                                                     aTGCCGTGGAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAAc               <  1:47731/139‑1 (MQ=255)
                                                                                                                     aTGCCGTGGAACCTGCCTGCCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAAc               <  1:167097/139‑1 (MQ=255)
                                                                                                                            ggAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                                                                                               <  1:326126/52‑1 (MQ=255)
                                                                                                                            ggAACCTGCCTGTCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAAcc                                                                                               >  2:326126/1‑52 (MQ=255)
                                                                                                                                 ctgcctacCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACg   <  2:65729/139‑1 (MQ=255)
                                                                                                                                  tgcctgtcGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCgg            <  1:282953/129‑1 (MQ=255)
                                                                                                                                  tgcctgtcGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCgg            >  2:282953/1‑129 (MQ=255)
                                                                                                                                  tgcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGa  >  2:145246/1‑139 (MQ=255)
                                                                                                                                   gcctgccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      >  2:7706/1‑102 (MQ=255)
                                                                                                                                   gcctgccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                      <  1:7706/102‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
GTCGACGACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGA  >  NZ_CP009273/49763‑50031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑