Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,624,8290TC63.3% 23.1 / 31.0 30E114E (GAA→GAG) ydfImannitol dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base C (9/10);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGG  >  NZ_CP009273/1624708‑1624941
                                                                                                                         |                                                                                                                
gCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCaa                                                                                                 <  1:22390/139‑1 (MQ=255)
    gtggGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACg                                                                                             >  2:212779/1‑139 (MQ=255)
                cACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCa                                                                                 <  2:19362/139‑1 (MQ=255)
                   gTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCatct                                                                              <  2:176896/139‑1 (MQ=255)
                                 tCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGTTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGc                                                                  >  2:483625/1‑137 (MQ=255)
                                        aTCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCtttt                                                         >  2:164988/1‑139 (MQ=255)
                                        aTCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGc                                                                  <  2:268493/130‑1 (MQ=255)
                                        aTCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGc                                                                  >  1:268493/1‑130 (MQ=255)
                                        aTCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCtttt                                                         >  1:1020/1‑139 (MQ=255)
                                        aTCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCtttt                                                         >  2:75655/1‑139 (MQ=255)
                                            aGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTaa                                                     <  2:334165/139‑1 (MQ=255)
                                              cATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAac                                                    >  1:349945/1‑138 (MQ=255)
                                               aTTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTa                                                                                      <  2:103303/103‑1 (MQ=255)
                                               aTTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTa                                                                                      >  1:103303/1‑103 (MQ=255)
                                                       tCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACg                                                                                             >  2:441890/1‑88 (MQ=255)
                                                       tCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACg                                                                                             <  1:441890/88‑1 (MQ=255)
                                                              gCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAg                                   <  1:164988/139‑1 (MQ=255)
                                                              gCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGa                                    <  1:75655/138‑1 (MQ=255)
                                                               ccGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGc                                  <  2:1020/139‑1 (MQ=255)
                                                               ccGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGc                                  <  1:252877/139‑1 (MQ=255)
                                                               ccGGAGAGTGGAAATACCCTTTTTCGGTGAATGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGc                                  <  1:368320/139‑1 (MQ=255)
                                                                        ggAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCcaca                          >  2:22432/1‑138 (MQ=255)
                                                                               cccTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGcc                  >  1:128807/1‑139 (MQ=255)
                                                                               cccTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGcc                  >  1:142909/1‑139 (MQ=255)
                                                                                 cTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGc                                  >  2:251769/1‑121 (MQ=255)
                                                                                 cTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGc                                  <  1:251769/121‑1 (MQ=255)
                                                                                    tttCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCa                                             >  1:101333/1‑107 (MQ=255)
                                                                                    tttCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCa                                             <  2:101333/107‑1 (MQ=255)
                                                                                            aTTGTCAGAGAGACAATCGCGATTTGCGGCTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCaa     <  2:387895/139‑1 (MQ=255)
                                                                                               gTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACgg  >  2:198793/1‑139 (MQ=255)
                                                                                                                         |                                                                                                                
GCTGGTGGGGATTTTGCACGTCGGCAGCTACCATCGGGTGATCGAGCATTAACTGTCCGGTCGCCGGAGAGTGGAAATACCCTTTTTCGGTGATTGTCAGAGAGACAATCGCGATTTGCGGTTCACACATCGCTGCCAACACGGTTTCTAAGCCATCTATCTGTACGTGCAAGGCTTTTTTAACGACGCCAACGACGCGAGCCGTCCACACATCGGCCGACATTTCCGCAACGG  >  NZ_CP009273/1624708‑1624941

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