Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273693,5420TG67.5% 45.1 / 41.3 40A364A (GCA→GCC) asnBasparagine synthase B
Reads supporting (aligned to +/- strand):  ref base T (7/6);  new base G (15/12);  total (22/18)
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.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGAATTTTTTGTCGAGGAACGGAACGCGTGCTTCCACGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGTCACATCA  >  NZ_CP009273/693405‑693677
                                                                                                                                         |                                                                                                                                       
aGGAATTTTTTGTCGAGGAACGGAACGCGTGCTTCCACGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTg                                                                                                                                        >  2:137263/1‑139 (MQ=255)
           gTCGAGGAACGGAACGCGTGCTTCCACGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGaa                                                                                                                              <  1:137263/138‑1 (MQ=255)
                               cTTCCACGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAacac                                                                                                         <  2:420536/139‑1 (MQ=255)
                                  ccACGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACt                                                                                                        <  2:773967/137‑1 (MQ=255)
                                     cGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTtcatc                                                                                                   <  2:805329/139‑1 (MQ=255)
                                       ccccAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTtcat                                                                                                    >  2:280539/1‑136 (MQ=255)
                                                 aCATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCa                                                                                         <  1:338082/137‑1 (MQ=255)
                                                  cATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTtgtg                                                                                                                                 <  1:318535/96‑1 (MQ=255)
                                                  cATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTtgtg                                                                                                                                 >  2:318535/1‑96 (MQ=255)
                                                  cATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATaa                                                                                                                    <  1:19299/109‑1 (MQ=255)
                                                  cATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATaa                                                                                                                    >  2:19299/1‑109 (MQ=255)
                                                                  cgcgcgCAGTCATACATATGCAGGGCCAGCATTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTaa                                                                      <  2:50368/139‑1 (MQ=255)
                                                                                   atGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGcc                                                           <  2:709424/133‑1 (MQ=255)
                                                                                   atGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGcc                                                           >  1:709424/1‑133 (MQ=255)
                                                                                     gCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCt                                                    >  2:332497/1‑138 (MQ=255)
                                                                                     gCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCCGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCt                                                    >  1:712657/1‑138 (MQ=255)
                                                                                     gCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCt                                                    >  2:483788/1‑138 (MQ=255)
                                                                                     gCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCt                                                    >  1:139672/1‑138 (MQ=255)
                                                                                     gCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCt                                                    >  1:659084/1‑138 (MQ=255)
                                                                                         ggCCAGCAGTTTACGCACCGTCTCTTCATGCGTTTCTTTGGCATTCGGGGCTTTTTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCa                                                                               <  2:396608/107‑1 (MQ=255)
                                                                                         ggCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTAcgc                                               >  2:648736/1‑139 (MQ=255)
                                                                                         ggCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCa                                                                               >  1:396608/1‑107 (MQ=255)
                                                                                         ggCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTAcgc                                               >  2:611925/1‑139 (MQ=255)
                                                                                           ccAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGa                                             >  1:295321/1‑139 (MQ=255)
                                                                                              gcagTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCa                                                                >  2:448816/1‑117 (MQ=255)
                                                                                              gcagTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCa                                                                <  1:448816/117‑1 (MQ=255)
                                                                                                  tttACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTaaa                                      <  1:483788/139‑1 (MQ=255)
                                                                                                    tACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATa                                    >  2:414981/1‑139 (MQ=255)
                                                                                                            gTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGgtg                            >  1:16039/1‑139 (MQ=255)
                                                                                                            gTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGgtg                            >  2:778322/1‑139 (MQ=255)
                                                                                                            gTCTCTTCATGCAGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGgtg                            >  1:319381/1‑139 (MQ=255)
                                                                                                                        aGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATAc                                   >  1:71200/1‑120 (MQ=255)
                                                                                                                        aGTTCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATAc                                   <  2:71200/120‑1 (MQ=255)
                                                                                                                          ttCTTTGGCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGaa                >  2:198891/1‑137 (MQ=255)
                                                                                                                                ggCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATGCATCGGTGTTGAAGCGCGAATAGTGGt         <  2:712657/138‑1 (MQ=255)
                                                                                                                                ggCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGt         <  2:383788/138‑1 (MQ=255)
                                                                                                                                ggCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGt         <  1:611925/138‑1 (MQ=255)
                                                                                                                                 gCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGTca       <  1:648736/139‑1 (MQ=255)
                                                                                                                                 gCATTCGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGTca       <  2:659084/139‑1 (MQ=255)
                                                                                                                                      cGGGGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGTCAcatca  >  2:823961/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                       
AGGAATTTTTTGTCGAGGAACGGAACGCGTGCTTCCACGCCCCAGGCTGACATCGCTTTGTTGGCACGCGCGCAGTCATACATATGCAGGGCCAGCAGTTTACGCACCGTCTCTTCATGCAGTTCTTTGGCATTCGGTGCTTTGTGGAAGTAAAGATAACCGCCGAACACTTCATCAGAACCTTCACCGGACAGCACCATTTTAATGCCCATCGCCTTGATCTTACGCGACATTAAATACATCGGTGTTGAAGCGCGAATAGTGGTCACATCA  >  NZ_CP009273/693405‑693677

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