Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,326,7180TA43.5% 25.5 / 13.2 23Y44F (TAT→TTT) adiCarginine/agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (13/0);  new base A (0/10);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 8.74e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATC  >  NZ_CP009273/4326650‑4326798
                                                                    |                                                                                 
cGAGGAACGACTTTTTGGCGTATCCCACCGAGAGCCCCAGTGCACCGATAATCGCCACCACCCATCCAAAAATGG‑CAATCCCGCCAGTag                                                             <  1:342837/90‑1 (MQ=255)
 gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTaga                                                            >  1:359386/1‑90 (MQ=255)
      aCGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCa                                                       >  2:171850/1‑90 (MQ=255)
                          accGAGACCCCAAGTGCCCCAAAAATCGCCACAACCCACCCAAAAATGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaa                                   <  1:409950/88‑1 (MQ=255)
                          aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaa                                   >  1:346920/1‑90 (MQ=255)
                          aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaa                                   >  2:468156/1‑90 (MQ=255)
                           tCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaaa                                  >  2:500852/1‑90 (MQ=255)
                           ccGAGACCCCCAGTCCCCCGAAAATCGCCCCCACCCACCAAAAAAGGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaaaa                                  <  2:34516/89‑1 (MQ=255)
                               gagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc                              >  1:34516/1‑90 (MQ=255)
                                 cccccAGTGCCCCAATAATCGCCCCAAACCACCAAAAAAGGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                            <  2:342357/89‑1 (MQ=255)
                                 cccccAGTGCACCGATAATCGCCACCAACCACCAAAAAAGGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                            <  2:336724/89‑1 (MQ=255)
                                 cccccAGTCCACCG‑TAAACCTCCCCAACCCACCAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                            <  2:346920/89‑1 (MQ=255)
                                  ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGa                           >  1:538292/1‑90 (MQ=255)
                                   cccAGTGCACGAATAATCGCCCCCACCCACCAAAAAAGGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                          <  2:562279/90‑1 (MQ=255)
                                   cccAGTCCACCAAAAACCGCCCCCACCCACCAAAAAAGGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                          <  2:21176/90‑1 (MQ=255)
                                   cccAGGCCACCAAAAACCGCCCCCACCCACCAAAAAAGGC‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                          <  2:10837/90‑1 (MQ=255)
                                     cAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAccc                        >  1:383652/1‑90 (MQ=255)
                                           aCCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                  >  1:21176/1‑90 (MQ=255)
                                            ccGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataata                 >  1:562279/1‑90 (MQ=255)
                                            ccGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataata                 >  1:153263/1‑90 (MQ=255)
                                            ccGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataata                 >  2:59143/1‑90 (MQ=255)
                                                 aaTCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTccc            >  1:28612/1‑90 (MQ=255)
                                                     gccaccacCCACCAAAAAAGGC‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac        <  2:28612/88‑1 (MQ=255)
                                                           acccatccaaAAATGG‑CACCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATc  <  2:359386/88‑1 (MQ=255)
                                                                    |                                                                                 
CGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATC  >  NZ_CP009273/4326650‑4326798

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