Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,326,7280TC55.6% 4.7 / 13.8 18I41V (ATT→GTT) adiCarginine/agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base C (0/10);  total (5/13)
Fisher's exact test for biased strand distribution p-value = 6.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATT  >  NZ_CP009273/4326645‑4326813
                                                                                   |                                                                                     
tGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGcc                                                                                 >  1:239633/1‑90 (MQ=255)
          aacaacATTTGGGCGTACCCCACCGAGCCCCCCAGTCCCCCGATAATCGCCCCCACCCACCAAAAAATGGCAACCCCGCCAGTAGAGGcc                                                                       <  2:407517/90‑1 (MQ=255)
          aaCGACTTTTTGGCGTATCCCACCGAGACCCCCAGTCCCCCGATAACCGCCACCACCCATCCAAAAATGGAAATCCCGCCAGTAGAGGcc                                                                       <  2:26137/90‑1 (MQ=255)
                         taccccAGCGAGCCCCCCAGTCCCCCGATAACGCCCCCCAACCACCAAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAAc                                                        <  2:311705/86‑1 (MQ=255)
                               aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaa                                                  >  1:382756/1‑90 (MQ=255)
                                     acccccAGGGCACCAATAACCGCCCCCACCCACCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCt                                            <  2:382756/88‑1 (MQ=255)
                                      gCCCCAGTGCCCCAATAATCGCCCCCAACCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                           <  1:140597/90‑1 (MQ=255)
                                      cccccAGTCCACGAAAAACCGCCCCCACCCACCCAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                           <  2:113991/89‑1 (MQ=255)
                                      cccccAGGGCAACGATAACCGCCCCCACCCACCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                           <  2:205320/89‑1 (MQ=255)
                                       ccccAGTCCACCGATAACCGCCCCCAACCACCAAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGa                                          <  2:41985/90‑1 (MQ=255)
                                        cccAGTGCCCCGAAAATCGCCCCCAACCCTCCAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                                         <  1:201497/90‑1 (MQ=255)
                                        cccAGTGCACCGAAAATCGTCACCAACCATCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                                         <  2:46234/90‑1 (MQ=255)
                                             gcccccGATAACCGCAACCACCCACCCAAAAACGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCat                                    <  1:304775/86‑1 (MQ=255)
                                                          gccACCAACCATCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                       <  1:357761/88‑1 (MQ=255)
                                                           acACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGaca                      >  2:368497/2‑90 (MQ=255)
                                                                            aTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTa                                                          <  2:101605/37‑1 (MQ=255)
                                                                            aTGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTa                                                          >  1:101605/1‑37 (MQ=255)
                                                                               gCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt  >  1:168398/1‑90 (MQ=255)
                                                                                   |                                                                                     
TGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATT  >  NZ_CP009273/4326645‑4326813

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: