Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,292,0170TC73.2% 51.4 / 33.1 41T521A (ACC→GCC) adhEbifunctional acetaldehyde‑CoA/alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base C (13/17);  total (18/23)
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.15e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGAGTTTCCGGATGTTCGTACATAACCCACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAGCGCACGTTTGTGGCCATCAGT  >  NZ_CP009273/1291883‑1292152
                                                                                                                                      |                                                                                                                                       
gtgtgtTTCCGGATGTTCGTACATAACCCACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAgtgccggg                                                                                                                                     <  1:646669/139‑1 (MQ=255)
gTGAGTTTCCGGATGTTCGTACATAACCCACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAgggccggg                                                                                                                                     <  1:312707/139‑1 (MQ=255)
    gTTTCCGGATGTTCGTACATAACCCACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCg                                                                                                                                 >  2:324974/1‑139 (MQ=255)
        ccGGATGTTCGTACATAACCCACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCtt                                                                                                                              <  1:522589/138‑1 (MQ=255)
                          cccACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACtt                                                                                                            >  2:829455/1‑138 (MQ=255)
                          cccACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACtt                                                                                                            >  2:466013/1‑138 (MQ=255)
                               aTGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGttt                                                                                                      >  1:38450/1‑139 (MQ=255)
                                            cgTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACC‑‑‑‑‑‑‑‑‑‑‑‑GCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACgg                                                                              >  1:274732/1‑138 (MQ=255)
                                             gTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTgctgct                                                                                        <  1:583513/139‑1 (MQ=255)
                                              tCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCtt                                                                                       >  1:856441/1‑139 (MQ=255)
                                                    ggggAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTa                                                                                 <  1:656774/139‑1 (MQ=255)
                                                           caccaccCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACC‑‑‑‑‑‑‑‑‑‑‑‑GCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGc                                                              <  2:274732/139‑1 (MQ=255)
                                                              caccCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgat                                                                       <  1:466013/139‑1 (MQ=255)
                                                                   aGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  <  1:780698/139‑1 (MQ=255)
                                                                    gcgcGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgatc                                                                 <  1:829455/139‑1 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTa                                                                                 >  2:683145/1‑122 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTa                                                                                 <  1:683145/122‑1 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGaa                                                                            >  2:753391/1‑127 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGaa                                                                            <  1:753391/127‑1 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  >  1:9556/1‑137 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  >  1:703082/1‑137 (MQ=255)
                                                                     cgcgATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  >  1:125726/1‑137 (MQ=255)
                                                                            aTCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAAc                                                         <  2:539819/139‑1 (MQ=255)
                                                                            aTCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAAc                                                         <  1:324974/139‑1 (MQ=255)
                                                                               aCGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCa                                                       <  1:139851/138‑1 (MQ=255)
                                                                                    tGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACgg                                                                              >  1:129588/1‑110 (MQ=255)
                                                                                    tGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACgg                                                                              <  2:129588/110‑1 (MQ=255)
                                                                                      gTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCgaag                                                                                                                   <  2:31172/71‑1 (MQ=255)
                                                                                      gTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCgaag                                                                                                                   >  1:31172/1‑71 (MQ=255)
                                                                                           aaGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAgg                                           >  2:757826/1‑138 (MQ=255)
                                                                                               aGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCg                                      <  2:694077/139‑1 (MQ=255)
                                                                                                 tttGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGgt                                    <  2:569688/139‑1 (MQ=255)
                                                                                                 tttACCAGTTCTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACAGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGgt                                    <  1:133203/139‑1 (MQ=255)
                                                                                                           cTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgat                                                                       <  2:692532/94‑1 (MQ=255)
                                                                                                           cTGCACCTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgat                                                                       >  1:692532/1‑94 (MQ=255)
                                                                                                              cACCTTTACGAACGATGCTCAGGGCCTGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  <  1:452283/96‑1 (MQ=255)
                                                                                                              cACCTTTACGAACGATGCTCAGGGCCTGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  >  2:452283/1‑96 (MQ=255)
                                                                                                                ccTTTACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAgcg                     <  2:76312/139‑1 (MQ=255)
                                                                                                                    tACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAGCGCa                   <  1:757826/137‑1 (MQ=255)
                                                                                                                    tACGAACGATGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAGCGCa                   <  2:399617/137‑1 (MQ=255)
                                                                                                                             tGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAGCGCACGTTTGTGGc         >  2:724468/1‑138 (MQ=255)
                                                                                                                             tGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  <  2:800002/81‑1 (MQ=255)
                                                                                                                             tGCTCAGGGCCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGtgatctgat                                                                  >  1:800002/1‑81 (MQ=255)
                                                                                                                                   gggtcgggtcCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAGCGCACGTTTGTGGCCATCAGt  <  1:724468/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GTGAGTTTCCGGATGTTCGTACATAACCCACATGATCTTCGCGGCGTCCATCGGGGAACCACCACCCAGCGCGATAATCACGTCTGGTTTGAAGGAGTTTGCCAGTTCTGCACCTTTACGAACGATGCTCAGGGTCGGGTCCGCTTCTACTTCGAAGAAGACTTCAGTTTCAACGCCTGCTGCTTTCAGTACGGAAGTGATCTGATCAGCATAACCATTGTTGAACAGGAAGCGGTCAGTCACGATGAGCGCACGTTTGTGGCCATCAGT  >  NZ_CP009273/1291883‑1292152

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