Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273146,0450TC62.5% 34.3 / 54.6 48E15G (GAA→GGA) panB3‑methyl‑2‑oxobutanoate hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (11/7);  new base C (16/14);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTATCCGTAGCAGCACTGCCA  >  NZ_CP009273/145916‑146177
                                                                                                                                  |                                                                                                                                    
aGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTg                                                                                                                              >  1:523553/1‑138 (MQ=255)
       aGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTg                                                                                                                      <  2:701890/139‑1 (MQ=255)
       aGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTg                                                                                                                      <  2:41807/139‑1 (MQ=255)
            tGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGt                                                                                                                  >  1:408038/1‑138 (MQ=255)
            tGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGt                                                                                                                  >  2:428288/1‑138 (MQ=255)
                         tCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                                                    <  2:347894/139‑1 (MQ=255)
                         tCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                                                    >  1:347894/1‑139 (MQ=255)
                         tCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                                                    >  1:475620/1‑139 (MQ=255)
                           aTGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTg                                                                                                                              >  2:789020/1‑111 (MQ=255)
                           aTGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTg                                                                                                                              <  1:789020/111‑1 (MQ=255)
                            tGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCgg                                                                                                 <  2:225275/139‑1 (MQ=255)
                              cccAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGtt                                                                                               >  2:193666/1‑139 (MQ=255)
                                   cGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCAt                                                                                           >  1:752261/1‑138 (MQ=255)
                                          cccACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATc                                                                                   <  1:193666/139‑1 (MQ=255)
                                                           aaGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGAtgttg                                                                  >  1:223076/1‑139 (MQ=255)
                                                             gCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                                                <  2:408038/139‑1 (MQ=255)
                                                             gCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                                                <  2:752261/139‑1 (MQ=255)
                                                                 ttCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                                                    >  1:78256/1‑99 (MQ=255)
                                                                 ttCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATggtggt                                                                                                    <  2:78256/99‑1 (MQ=255)
                                                                  tcatcaGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCa                                                                                                                    >  2:554305/1‑82 (MQ=255)
                                                                  tcatcaGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCa                                                                                                                    <  1:554305/82‑1 (MQ=255)
                                                                  tcatcaGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTc                                                           >  2:243607/1‑139 (MQ=255)
                                                                         cAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACtt                                                                                                                         <  2:847982/70‑1 (MQ=255)
                                                                         cAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACtt                                                                                                                         >  1:847982/1‑70 (MQ=255)
                                                                             gagTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTAt                                                <  1:428288/139‑1 (MQ=255)
                                                                             gagTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTAt                                                <  1:243607/139‑1 (MQ=255)
                                                                              agTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTAt                                                >  2:461667/1‑138 (MQ=255)
                                                                                   ggCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCa                                                      >  1:835739/1‑127 (MQ=255)
                                                                                   ggCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCa                                                      <  2:835739/127‑1 (MQ=255)
                                                                                     cGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTcaca                                        <  1:814667/139‑1 (MQ=255)
                                                                                                  aTAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCaa                           >  1:301613/1‑139 (MQ=255)
                                                                                                  aTAAGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCa                            <  2:585235/138‑1 (MQ=255)
                                                                                                    aaGCGGTGATGGTCGCGAAACGTTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAAt                          >  1:700424/1‑139 (MQ=255)
                                                                                                    aaGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAAt                          >  2:547114/1‑138 (MQ=255)
                                                                                                    aaGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAAt                          >  1:305060/1‑138 (MQ=255)
                                                                                                    aaGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAAt                          >  2:658906/1‑138 (MQ=255)
                                                                                                    aaGCGGTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAAt                          >  1:569751/1‑138 (MQ=255)
                                                                                                     aGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                                                <  2:755346/99‑1 (MQ=255)
                                                                                                     aGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTtgtg                                                                >  1:755346/1‑99 (MQ=255)
                                                                                                         gTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTa                    >  1:638801/1‑139 (MQ=255)
                                                                                                         gTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTc                                                                                             >  1:676490/1‑66 (MQ=255)
                                                                                                         gTGATGGTCGCGAAACGTTTTTTTCCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTc                                                                                             <  2:676490/66‑1 (MQ=255)
                                                                                                          tGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTAt                   >  1:88303/1‑139 (MQ=255)
                                                                                                            aTGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTATc                  <  2:447528/138‑1 (MQ=255)
                                                                                                                 cgcgAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCt                                                          >  1:119190/1‑93 (MQ=255)
                                                                                                                 cgcgAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCt                                                          <  2:119190/93‑1 (MQ=255)
                                                                                                                    gAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTATCCGTagcagc         <  2:638801/139‑1 (MQ=255)
                                                                                                                    gAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTATCCGTagcag          >  2:544373/1‑138 (MQ=255)
                                                                                                                             tttttcCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTATCCGTAGCAGCACTGCCa  <  2:281834/138‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
AGGGTGGAGTCGTGCCCCTGAACCGTCATGCCCAGCGAATCGCCCACCAGCATGACGTTAAGCCCTTCATCAGCAAAGAGTTTGGCGAAGCTATAGTCATAAGCGGTGATGGTCGCGAAACGTTTTTTTTCCTGTTTGTACTTCTGCAGTAAGGAGATGGTGGTCGGTTTCATAACGTATCCTGATAAATTGATGTTGTGCTGTCTGGCATTTTATCAGTCACATTGGTGGGGGCAATGATTTATCCGTAGCAGCACTGCCA  >  NZ_CP009273/145916‑146177

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