Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,144,2090AG64.5% 26.3 / 30.9 31I189I (ATT→ATC) argEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (7/4);  new base G (11/9);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCGGCAAAATAACGCGCTCCGGCCATACTGGTTT  >  NZ_CP009273/4144072‑4144342
                                                                                                                                         |                                                                                                                                     
ggCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGaa                                                                                                                                      >  2:134898/1‑139 (MQ=255)
ggCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGaa                                                                                                                                      >  1:553220/1‑139 (MQ=255)
     ggTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACg                                                                                                                                  >  2:813484/1‑138 (MQ=255)
     ggTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACg                                                                                                                                  >  2:797956/1‑138 (MQ=255)
     ggTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACg                                                                                                                                  >  2:780702/1‑138 (MQ=255)
     ggTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACg                                                                                                                                  >  2:575457/1‑138 (MQ=255)
       tAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGaa                                                                                                                                      >  2:94617/1‑132 (MQ=255)
       tAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGaa                                                                                                                                      <  1:94617/132‑1 (MQ=255)
       tAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGAt                                                                                                                               >  2:409889/1‑139 (MQ=255)
       tAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGAt                                                                                                                               >  2:405828/1‑139 (MQ=255)
                cGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTa                                                                                                                        >  1:428598/1‑137 (MQ=255)
                    gTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGata                                                                                                                  <  2:218695/139‑1 (MQ=255)
                    gTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGata                                                                                                                  <  2:549753/139‑1 (MQ=255)
                    gTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGata                                                                                                                  <  1:55669/139‑1 (MQ=255)
                     tGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGata                                                                                                                  <  1:405828/138‑1 (MQ=255)
                      gATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGatat                                                                                                                 >  1:316189/1‑138 (MQ=255)
                       aTAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGAt                                                                                                                                     >  2:52620/1‑117 (MQ=255)
                       aTAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGAt                                                                                                                                     <  1:52620/117‑1 (MQ=255)
                                      gTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGcgcg                                                                                                 >  2:387971/1‑138 (MQ=255)
                                         aTCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTAc                                                                                              <  1:41861/138‑1 (MQ=255)
                                                                           tGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGgcgc                                                           <  2:141616/139‑1 (MQ=255)
                                                                           tGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGgcgc                                                           <  1:81892/139‑1 (MQ=255)
                                                                               ttAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATc                                                       <  1:813484/139‑1 (MQ=255)
                                                                                  gTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCgg                                                    >  1:246902/1‑139 (MQ=255)
                                                                                   ttCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCggg                                                   >  1:437042/1‑139 (MQ=255)
                                                                                                     gcgTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTa                                                                                                                        >  2:216401/1‑52 (MQ=255)
                                                                                                     gcgTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTa                                                                                                                        <  1:216401/52‑1 (MQ=255)
                                                                                                     gcgTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGt                                     >  1:802476/1‑135 (MQ=255)
                                                                                                     gcgTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGt                                     <  2:802476/135‑1 (MQ=255)
                                                                                                      cgTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCgg                                >  2:58752/1‑139 (MQ=255)
                                                                                                      cgTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCgg                                >  2:838331/1‑139 (MQ=255)
                                                                                                                 gCTGGAGTGCCCCGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCGGCAAAATAAcgc                     <  1:58752/139‑1 (MQ=255)
                                                                                                                          ccccGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCGGCAAAATAACGCGCTCCGGc             >  1:96502/1‑138 (MQ=255)
                                                                                                                             cGACTGGCCCTGGATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCGGCAAAATAACGCGCTCCGGCCATa         <  1:387971/139‑1 (MQ=255)
                                                                                                                                    cccTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCGGCAAAATAACGCGCTCCGGCCATACTGGttt  >  2:815530/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GGCACGGTAAACGCTTCGTAGTGATAACGTTCTTTCAGGTTATCGCGCAATTGCAAAATATGCCCGATGGCGTCGTGCATTAGTTCGATAGCGTTAACTCCGCGTGCTGGATCGCTGGAGTGCCCCGACTGGCCCTGAATACGGATGGCGTTAGAGATATGACCTTTATGTGCGCGTACCGGTTGTAGTGACGTCGGTTCGCCAATGATGGCGCAATCCGGGCGCAGGGCGGTAGTTTCGGCAAAATAACGCGCTCCGGCCATACTGGTTT  >  NZ_CP009273/4144072‑4144342

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