Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,510,3060AG72.7% 80.3 / 43.3 55K158E (AAA→GAA) patDaminobutyraldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (9/6);  new base G (24/16);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGG  >  NZ_CP009273/1510182‑1510439
                                                                                                                            |                                                                                                                                     
cTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGgcg                                                                                                                         >  1:461539/1‑139 (MQ=255)
     tGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTggc                                                                                                                    >  1:476870/1‑139 (MQ=255)
         tGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCa                                                                                                                <  2:491012/139‑1 (MQ=255)
             cTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGgcgc                                                                                                                        >  2:388137/1‑127 (MQ=255)
             cTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGgcgc                                                                                                                        <  1:388137/127‑1 (MQ=255)
             cTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAggg                                                                                                             >  1:791245/1‑138 (MQ=255)
             cTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAggg                                                                                                             >  1:282661/1‑138 (MQ=255)
                    cAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGt                                                                                                     <  1:479252/139‑1 (MQ=255)
                                 tGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCa                                                                                        >  2:840125/1‑139 (MQ=255)
                                 tGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCa                                                                                        >  1:255357/1‑139 (MQ=255)
                                      gTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAg                                                                                    >  1:466453/1‑138 (MQ=255)
                                        cATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGaaa                                                                                 >  2:410238/1‑139 (MQ=255)
                                         aTACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCa                                                                                                                >  2:88965/1‑107 (MQ=255)
                                         aTACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCa                                                                                                                <  1:88965/107‑1 (MQ=255)
                                         aTACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAAt                                                                                >  2:449111/1‑139 (MQ=255)
                                                    atCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGa                                                                     <  2:466453/139‑1 (MQ=255)
                                                       cGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCa                                                                                     >  1:455378/1‑120 (MQ=255)
                                                       cGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCa                                                                                     <  2:455378/120‑1 (MQ=255)
                                                       cGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACcg                                                                  >  1:661478/1‑139 (MQ=255)
                                                                   ttGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTgg                                                      >  1:220239/1‑139 (MQ=255)
                                                                             tGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCt                                               <  1:440458/136‑1 (MQ=255)
                                                                             tGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCt                                               >  2:440458/1‑136 (MQ=255)
                                                                                    tATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCa                                                                                                                <  2:115630/64‑1 (MQ=255)
                                                                                    tATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCa                                                                                                                >  1:115630/1‑64 (MQ=255)
                                                                                      tCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGata                                    <  1:840125/138‑1 (MQ=255)
                                                                                         cACCGTGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCtt                                >  2:829210/1‑139 (MQ=255)
                                                                                              tGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGTACTGCGTa                                                                                                    >  1:141897/1‑66 (MQ=255)
                                                                                              tGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTa                                                                                                    <  2:141897/66‑1 (MQ=255)
                                                                                              tGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCgg                           >  2:49471/1‑139 (MQ=255)
                                                                                              tGGAATTATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCgg                           >  1:160189/1‑139 (MQ=255)
                                                                                              tGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCgg                           >  1:767020/1‑139 (MQ=255)
                                                                                              tGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCgg                           >  2:302633/1‑139 (MQ=255)
                                                                                               ggAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCgg                           >  2:582831/1‑138 (MQ=255)
                                                                                                   ttATCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGc                      <  1:410238/139‑1 (MQ=255)
                                                                                                      tCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTa                                                                                                    <  2:29460/58‑1 (MQ=255)
                                                                                                      tCCGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTa                                                                                                    >  1:29460/1‑58 (MQ=255)
                                                                                                      tCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGt                    >  2:706451/1‑138 (MQ=255)
                                                                                                        cGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTggcg                                                                                                                   >  1:773254/1‑41 (MQ=255)
                                                                                                        cGCTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTggcg                                                                                                                   <  2:773254/41‑1 (MQ=255)
                                                                                                          cTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCa                                                    <  1:140883/102‑1 (MQ=255)
                                                                                                          cTGATGATGGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCa                                                    >  2:140883/1‑102 (MQ=255)
                                                                                                             atgatgGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCATGCGTGATTAACa            <  2:282661/139‑1 (MQ=255)
                                                                                                             atgatgGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACa            <  1:449111/139‑1 (MQ=255)
                                                                                                             atgatgGCCGCGTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACa            >  1:70442/1‑139 (MQ=255)
                                                                                                             atgatgGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACa            <  1:302633/139‑1 (MQ=255)
                                                                                                               gatgGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATa          <  2:661478/139‑1 (MQ=255)
                                                                                                                atgGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATa          <  1:706451/138‑1 (MQ=255)
                                                                                                                     cgcgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTgg                                             >  1:672781/1‑98 (MQ=255)
                                                                                                                     cgcgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTgg                                             <  2:672781/98‑1 (MQ=255)
                                                                                                                     cgcgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGt                    >  1:85265/1‑123 (MQ=255)
                                                                                                                     cgcgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGt                    <  2:85265/123‑1 (MQ=255)
                                                                                                                       cgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGAc                                                                    <  2:767306/73‑1 (MQ=255)
                                                                                                                       cgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGAc                                                                    >  1:767306/1‑73 (MQ=255)
                                                                                                                       cgTGGGAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTgg  >  1:207328/1‑139 (MQ=255)
                                                                                                                       cgTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTgg  >  2:729999/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                     
CTGTCTGAATGGTCTGGCGGCAGGTGAATATCTTGAAGGTCATACTTCGATGATCCGTCGCGATCCGTTGGGGGTCGTGGCTTCTATCGCACCGTGGAATTATCCGCTGATGATGGCCGCGTGGAAACTTGCTCCGGCGCTGGCGGCAGGGAACTGCGTAGTGCTTAAACCATCAGAAATTACCCCGCTGACCGCGTTGAAGTTGGCAGAGCTGGCGAAAGATATCTTCCCGGCAGGCGTGATTAACATACTGTTTGG  >  NZ_CP009273/1510182‑1510439

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: ‑