Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,171,1070CA46.7% 97.3 / 32.4 60E55* (GAG→TAG) nudFADP‑ribose diphosphatase
Reads supporting (aligned to +/- strand):  ref base C (29/3);  new base A (0/28);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.36e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCC  >  NZ_CP009273/3171025‑3171191
                                                                                   |                                                                                    
gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt                                                                               >  2:1052965/1‑90 (MQ=255)
gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt                                                                               >  1:983144/1‑90 (MQ=255)
   gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc                                                                            >  1:790017/1‑90 (MQ=255)
   gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc                                                                            >  2:1898441/1‑90 (MQ=255)
      aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                         >  2:962206/1‑90 (MQ=255)
      aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGACGTGACCGCGCTCAAAAATTTCccgc                                                                         >  1:478912/1‑90 (MQ=255)
      aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGACGTGACCGCGCTCAAAAATTTCccgc                                                                         >  1:261168/1‑90 (MQ=255)
         cGAATCTGCTCAATCAGCACACCTTCACCCCGCCCGGGGCAAAGGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCcgccgc                                                                      <  1:1501867/90‑1 (MQ=255)
         cGAATCTGCTAAATCAGCACAACTCCACCCCGCACTGGGTCAAGGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTCCCcgccgc                                                                      <  1:307083/90‑1 (MQ=255)
            aTCTGCTCAATCAGCACAACTTCATCACGCACGGGGTAAAGGGGTAGCAAGACTGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCAcc                                                                   <  2:2118238/90‑1 (MQ=255)
            aTCTGCTCAATCAGCACAACTTCACCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCAcc                                                                   <  2:1993912/90‑1 (MQ=255)
                gCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCAt                                                               >  2:810341/1‑90 (MQ=255)
                  tCAATCAGCACACCTTCATCCCGCACGGGGCAAAGGGGTAGCAAGACTGGGGCGTGGCGGCGCTAAAAAATTTCCCGCCGCACCTCATGa                                                             <  2:49176/90‑1 (MQ=255)
                        aGCACAACTTCACCACGCACGGGGCCAAGGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCAtt                                                       <  1:510481/90‑1 (MQ=255)
                          cacaACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTg                                                     >  1:2014039/1‑90 (MQ=255)
                            caACTTCATCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc                                                   <  1:961725/90‑1 (MQ=255)
                            caACTCCACCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc                                                   <  1:757852/90‑1 (MQ=255)
                              aCTTCATCCCGAACGGGGTAAAAGGGAAGAAAGATGGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCg                                                 <  2:689612/90‑1 (MQ=255)
                               cTTCATCCCGCACTGGGCAAAGGGGTAGCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                <  1:1514768/90‑1 (MQ=255)
                               cTTCACCCCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                <  2:684869/90‑1 (MQ=255)
                               cTTCACCACGCACGGGGTCAAGGGGTAGCAAGCCTGCGGGGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGt                                                <  1:800722/90‑1 (MQ=255)
                                ttCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGtt                                               >  1:233345/1‑90 (MQ=255)
                                 tcatcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                              >  1:1578106/1‑90 (MQ=255)
                                 tcatcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCCCCGCACCTCATGACTCATTTGCCCGTTg                                              <  1:1451163/90‑1 (MQ=255)
                                      gcGCACTGGGCAAAAGGGAAGCAGGTCGGGGGGGGGCCCGCGTAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                         <  1:1687332/89‑1 (MQ=255)
                                      ccGCCCTGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                         <  2:1182056/89‑1 (MQ=255)
                                      aCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                         >  1:120544/1‑90 (MQ=255)
                                        gCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc                                       >  2:1017529/1‑90 (MQ=255)
                                        gCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc                                       >  1:1343672/1‑90 (MQ=255)
                                        ccACTGGGTCAAGGGGAACAAAGACCGCGGCGTGCCCGCGCAAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc                                       <  1:1160383/89‑1 (MQ=255)
                                         cACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                      >  1:384154/1‑90 (MQ=255)
                                         cACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                      >  2:800722/1‑90 (MQ=255)
                                          ccTGGGTCAAGGGGAAGCAAGCCGGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                     <  1:1398068/89‑1 (MQ=255)
                                          aCTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                     <  1:545768/90‑1 (MQ=255)
                                          aCTGGGTCAAAGGGAAGAAAGACTGCGACGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                     <  2:261168/90‑1 (MQ=255)
                                          aCTGGGTAAAAGGGTAGCAAGACTGCGACGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGa                                     <  2:478912/90‑1 (MQ=255)
                                               gTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                >  1:287539/1‑90 (MQ=255)
                                                taaaaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                               <  2:1721213/88‑1 (MQ=255)
                                                 cAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGac                              >  2:1687332/1‑89 (MQ=255)
                                                 cAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                              >  1:354992/1‑90 (MQ=255)
                                                 aaaggggTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaa                              <  2:695306/86‑1 (MQ=255)
                                                  aaggggTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  2:2156584/87‑1 (MQ=255)
                                                  aaaGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  2:1792756/90‑1 (MQ=255)
                                                  aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  2:384154/90‑1 (MQ=255)
                                                  aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  2:911057/90‑1 (MQ=255)
                                                  aaaGGGTAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  2:1301006/90‑1 (MQ=255)
                                                   aaGGGTTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  2:287539/90‑1 (MQ=255)
                                                   aaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt                            >  2:1060291/1‑90 (MQ=255)
                                                    aGGGTAGCAAGACTGCGGCGTGCCCGCGAAAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATc                           <  1:1795654/90‑1 (MQ=255)
                                                     gggTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                          >  2:2129208/1‑90 (MQ=255)
                                                     gggTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                          <  2:1613849/90‑1 (MQ=255)
                                                       gTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat                        >  2:997198/1‑90 (MQ=255)
                                                         tAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                       >  1:1225020/1‑90 (MQ=255)
                                                                  tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTaa              >  2:200923/1‑90 (MQ=255)
                                                                  tGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTaa              >  1:271497/1‑90 (MQ=255)
                                                                   gcggcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt             >  2:1174792/1‑90 (MQ=255)
                                                                        gTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             <  1:166661/69‑1 (MQ=255)
                                                                        gTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGaaa                             >  2:166661/1‑69 (MQ=255)
                                                                              gcgcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc  >  1:1558999/1‑90 (MQ=255)
                                                                              gcgcTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGcc  >  1:1948761/1‑90 (MQ=255)
                                                                                   |                                                                                    
GCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCC  >  NZ_CP009273/3171025‑3171191

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: