Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,436,8530TC63.3% 39.9 / 49.7 54D83G (GAC→GGC) ldhAD‑lactate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (8/10);  new base C (16/15);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATACGCGCGGTGAATACGGCGGTTCAGCGTCATCATCATACCGATGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCGCTTCGCAGCCATTG  >  NZ_CP009273/1436729‑1436975
                                                                                                                            |                                                                                                                          
aTACGCGCGGTGAATACGGCGGTTTCGCGGCATCATCATACCCATGGCGGGGTCAGGAACGGGCTCTGGGTCATAGGCTGGGACACGGGCTACTTTTAGCCCCCGTTCTTTTGCCGCGTCAAGGGCGGCGTTATTGgaa                                                                                                              >  2:849620/1‑139 (MQ=255)
aTACGCGCGGTGAATACGGCGGTTCAGCGTCATCATCATACCGATGGCGGGTTCAGCAACCGCCTCTGGATCATAGGCCGGAACACGGACTACTTTCAACCCCAGTTCTTTTTCCGCGTCCAGGTCGACGTTATTGaaa                                                                                                              >  1:180992/1‑139 (MQ=255)
aTACGCGCGGTGAATACGGCGGTTCAGCGTCATCATCATAACGATGGCGTGTTCAGCAACGGGCTCTGGATCATAAGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTTCCGCGTCAAGGCCGGCGTTATTGaa                                                                                                               >  1:420134/1‑138 (MQ=255)
     cgcgGTGAATACGGCGGGTTAGCGTCATCATCCTACCGATGGCGGGTTCCGCAAAGGCCTCTGGATCATAAGCTGGCAAACGGACAACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGAGGGTTTTGGAAACGGc                                                                                                         >  1:637792/1‑139 (MQ=255)
      gcgGTGAATACGGCGGTTCAGCGTCATCATCATACCGATGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGACGACGTTATTGAAACCGGca                                                                                                        <  1:508129/139‑1 (MQ=255)
       cgGTGAATACGGCGGTTCAGCGTCATCATCATACCGCTGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCCGTTCTTTTTCCGCGTCAAGGTCGACGTTATTGAAACCGGccc                                                                                                       >  1:563095/1‑137 (MQ=255)
                  gcggTTCAGCGTCATCATCATACCGATGGGGGGTTCAGCAACGGCCTCTGGATCATAGGCTGGAAAACGGACCACTTTCAGCCCCAGTTCTTTTTCCGCGTCAAGGGCGGCGGTATTGAAAACGGCCCAGAGCAggggg                                                                                            >  2:49697/1‑137 (MQ=255)
                   cggTTCAGCGTCATCATCATACCGATGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGa                                                                                           <  2:671545/139‑1 (MQ=255)
                      ttCAGCGTCATCATCATACCGATGGCGTGTTCAGCAAAGGCCTCTGGATCATAGGGGGGAACACGGGCTACCTTCAGCCCCCGCTCTTTTGCCGCGTCCAGGGCGACGTTATTGAAAACGGCACAGCGCCGGGCGgta                                                                                         >  2:827844/1‑138 (MQ=255)
                                   tcatACCGATGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTg                                                                           <  1:547361/139‑1 (MQ=255)
                                           aTGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCttttt                                                                     >  2:600238/1‑137 (MQ=255)
                                           aTGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCg                                                                                            >  1:422276/1‑114 (MQ=255)
                                           aTGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCg                                                                                            <  2:422276/114‑1 (MQ=255)
                                              gCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 <  1:770556/138‑1 (MQ=255)
                                              gCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGat                                                                                          <  1:2981/113‑1 (MQ=255)
                                              gCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGat                                                                                          >  2:2981/1‑113 (MQ=255)
                                              gCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 <  2:384575/138‑1 (MQ=255)
                                              gCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 <  1:849620/138‑1 (MQ=255)
                                              gCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 <  1:554096/138‑1 (MQ=255)
                                                   ttCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTAtt                                                                                                                  <  2:4907/84‑1 (MQ=255)
                                                   ttCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTAtt                                                                                                                  >  1:4907/1‑84 (MQ=255)
                                                   ttCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTc                                                            >  2:73680/1‑138 (MQ=255)
                                                   ttCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTc                                                            >  1:81131/1‑138 (MQ=255)
                                                        cAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAggg                                                                                              >  2:497658/1‑99 (MQ=255)
                                                        cAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAggg                                                                                              <  1:497658/99‑1 (MQ=255)
                                                             gCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCa                                                                                                 >  1:461374/1‑91 (MQ=255)
                                                             gCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCa                                                                                                 <  2:461374/91‑1 (MQ=255)
                                                                  tGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGggcggc                                              <  2:617459/137‑1 (MQ=255)
                                                                  tGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 <  2:270111/118‑1 (MQ=255)
                                                                  tGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 >  1:270111/1‑118 (MQ=255)
                                                                       cATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGt                                        >  2:480490/1‑138 (MQ=255)
                                                                        aTAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCa                                      >  1:667076/1‑139 (MQ=255)
                                                                         tAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 <  1:811909/111‑1 (MQ=255)
                                                                         tAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGc                                                                 >  2:811909/1‑111 (MQ=255)
                                                                          aGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATc                                    >  2:186083/1‑139 (MQ=255)
                                                                          aGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAgg                                                                                               <  1:673392/80‑1 (MQ=255)
                                                                          aGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAgg                                                                                               >  2:673392/1‑80 (MQ=255)
                                                                                gAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTa                               <  2:397702/138‑1 (MQ=255)
                                                                                gAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTa                               <  2:81131/138‑1 (MQ=255)
                                                                                    acGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGaaa                          <  1:186083/139‑1 (MQ=255)
                                                                                      ggACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGat                                                                                          >  1:54184/1‑73 (MQ=255)
                                                                                      ggACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGat                                                                                          <  2:54184/73‑1 (MQ=255)
                                                                                      ggACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGaaaa                         >  2:399215/1‑138 (MQ=255)
                                                                                        actactTTCAGCCCCAGTTCTTTTGCCTCGTCAAGGTCGACGTTATTGAAACCGGca                                                                                                        <  1:854679/57‑1 (MQ=255)
                                                                                        actactTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGGCGACGTTATTGAAACCGGca                                                                                                        >  2:854679/1‑57 (MQ=255)
                                                                                         ctactTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAataca                     <  1:480490/139‑1 (MQ=255)
                                                                                           actTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAatacata                   >  1:432317/1‑139 (MQ=255)
                                                                                           actTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAatacata                   >  1:217500/1‑139 (MQ=255)
                                                                                            ctTTCAGCCCCAGTGCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAatacatac                  <  2:819716/139‑1 (MQ=255)
                                                                                               tCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCg                <  1:827844/138‑1 (MQ=255)
                                                                                               tCAGCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCg                <  1:73680/138‑1 (MQ=255)
                                                                                                  gCCCCAGTTCTTTTGCCGCGTCAAGGCCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCGCTTc            >  2:837484/1‑139 (MQ=255)
                                                                                                        gTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCGCTTCGCAGc       >  1:125772/1‑138 (MQ=255)
                                                                                                              ttGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCGCTTCGCAGCCATTg  <  1:49697/137‑1 (MQ=255)
                                                                                                                            |                                                                                                                          
ATACGCGCGGTGAATACGGCGGTTCAGCGTCATCATCATACCGATGGCGTGTTCAGCAACGGCCTCTGGATCATAGGCTGGAACACGGACTACTTTCAGCCCCAGTTCTTTTGCCGCGTCAAGGTCGACGTTATTGAAACCGGCACAGCGCAGGGCGATATATTTAACGCCGTGCTTTTTCAGCTCTTCCAGCACCGGGCGGCTGCCGTCATCGTTTACGAAAATACATACCGCTTCGCAGCCATTG  >  NZ_CP009273/1436729‑1436975

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