Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927349,9030TA58.1% 14.3 / 38.7 31D27E (GAT→GAA) folAtype 3 dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (6/7);  new base A (10/8);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGATCGCG  >  NZ_CP009273/49769‑50036
                                                                                                                                      |                                                                                                                                     
gacgGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGGACTCg                                                                                                                                   >  1:369684/1‑139 (MQ=255)
                    taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAcc                                                                                                               <  2:246729/139‑1 (MQ=255)
                    taGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACAcc                                                                                                               >  1:29355/1‑139 (MQ=255)
                          cGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAAAGCCATGCCGTGGCACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTtaaa                                                                                                          >  2:374737/1‑138 (MQ=255)
                            aCAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTtaaata                                                                                                        >  2:117892/1‑138 (MQ=255)
                                           gggAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                                                        <  2:187192/139‑1 (MQ=255)
                                           gggAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATggg                                                                                        <  1:265378/139‑1 (MQ=255)
                                                       tGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg                                                                            >  2:191445/1‑139 (MQ=255)
                                                         aTCAGTCTGATTGCGGCGTTAGCGGTAGCGCGCGTTATCGGCATGGAAAACGCCATGCGGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCATCTCCTTAAATAAACCCGTTATTATGGGCCGCCATTCCTgg                                                                           <  1:374737/138‑1 (MQ=255)
                                                             gTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaat                                                                       >  2:5619/1‑138 (MQ=255)
                                                                     gcggcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt                                                               <  1:5619/138‑1 (MQ=255)
                                                                        gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt                                                               <  1:246732/135‑1 (MQ=255)
                                                                        gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgt                                                               >  2:246732/1‑135 (MQ=255)
                                                                        gcgTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGGTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGgtcgtc                                                           >  1:20297/1‑139 (MQ=255)
                                                                                       cgcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                                                                   >  2:438170/1‑116 (MQ=255)
                                                                                       cgcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCAGTGATTATGGGCCGCCATACCTGGGaatcaat                                                                   <  1:438170/116‑1 (MQ=255)
                                                                                        gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                           <  2:290050/139‑1 (MQ=255)
                                                                                         cgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaa                                           >  1:60841/1‑138 (MQ=255)
                                                                                          gTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCaaa                                          >  1:176281/1‑138 (MQ=255)
                                                                                                                 gCCGTGGAACCTGCCTACCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCg                  >  2:471025/1‑139 (MQ=255)
                                                                                                                       gAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                                                                   >  1:117622/1‑84 (MQ=255)
                                                                                                                       gAACCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                                                                   <  2:117622/84‑1 (MQ=255)
                                                                                                                         aCCTGCCTGCCGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGa          <  2:176281/139‑1 (MQ=255)
                                                                                                                           ctgcctacCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACg        <  1:191445/139‑1 (MQ=255)
                                                                                                                           ctgcctacCGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACg        <  2:60841/139‑1 (MQ=255)
                                                                                                                            tgcctgccGAACTCGCCGGGTTTAAACGCAACAcc                                                                                                               <  2:111798/35‑1 (MQ=255)
                                                                                                                            tgcctgccGAACTCGCCGGGTTTAAACGCAACAcc                                                                                                               >  1:111798/1‑35 (MQ=255)
                                                                                                                             gcctgccGAACTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGAt      >  2:51848/1‑139 (MQ=255)
                                                                                                                             gcctaccGATCTCGCCGGGTTTAAACGCAACCCCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                                                                   <  2:145471/78‑1 (MQ=255)
                                                                                                                             gcctaccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatcaat                                                                   >  1:145471/1‑78 (MQ=255)
                                                                                                                                  ccGATCTCGCCGGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGATcgcg  <  2:29355/138‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
GACGGTTTACGCTTTACGTATAGTGGCGACAATTTTTTTTATCGGGAAATCTCAATGATCAGTCTGATTGCGGCGTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCAATCGGTCGTCCGTTGCCAGGACGCAAAAATATTATCCTCAGCAGTCAACCGGGTACGGACGATCGCG  >  NZ_CP009273/49769‑50036

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