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

AAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGG  >  NZ_CP009273/4326629‑4326802
                                                                                         |                                                                                    
aagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                      >  2:676075/1‑90 (MQ=255)
     ccaccaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATg                                                                                 >  1:544262/1‑90 (MQ=255)
       accaGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGc                                                                               >  2:609829/1‑90 (MQ=255)
           ggACTTGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGAGCCCAAGTGCACCGATAATCGTCACCAACCATCCAAAAATGGCAATc                                                                           <  1:40015/90‑1 (MQ=255)
               ttGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATgg                                                                                <  2:443957/81‑1 (MQ=255)
               ttGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATgg                                                                                >  1:443957/1‑81 (MQ=255)
                tGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGcc                                                                      >  2:704390/1‑90 (MQ=255)
                          aaCGACTTTTTGGCGTATCCCACGGAGACCCCCAGTCCCCCGATAATCGCCCCCACCCACCCATAAATGGCAATCCCGCCAGTAGAGGcc                                                            <  2:178136/90‑1 (MQ=255)
                                  ttgggCGTACCCCAGCGAGCCCCCCAGTGCCCCGAAAATCGCCCCCACCCATCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGc                                                    <  2:456800/87‑1 (MQ=255)
                                                  gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAaca                                    >  1:150712/1‑90 (MQ=255)
                                                    gagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc                                  >  1:325017/1‑90 (MQ=255)
                                                    gagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc                                  >  2:106565/1‑90 (MQ=255)
                                                      gCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                >  1:240847/1‑90 (MQ=255)
                                                      ccccaAGTGCACCGAAAACCTCCACCAACCACAAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                <  2:439995/89‑1 (MQ=255)
                                                       ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGa                               >  2:57885/1‑90 (MQ=255)
                                                        cccAGTGCCCCGATAACCGCCACCACCCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                              <  1:488981/90‑1 (MQ=255)
                                                        cccAGTGCCCCGAAAACCGTCCCCACCCATCCAAAAATGCCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                              <  2:447483/90‑1 (MQ=255)
                                                        cccAGTGCACCGATAATCGTCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                              <  1:57027/90‑1 (MQ=255)
                                                        cccAGTGCACCGATAATCGTCACCAACCACCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                              <  1:148052/90‑1 (MQ=255)
                                                          aaGTGCACCGATAATCGTCACCAACCATCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAccc                            <  1:422727/89‑1 (MQ=255)
                                                             cccACCGATAACCGCCCCCACCCACCAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCat                         <  2:693382/88‑1 (MQ=255)
                                                                aCCGATAATCGTCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                      <  2:478397/90‑1 (MQ=255)
                                                                      aaTCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTccc                >  1:70744/1‑90 (MQ=255)
                                                                          gccaccacCCATCAAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac            <  1:212012/88‑1 (MQ=255)
                                                                           tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                              >  2:520967/1‑71 (MQ=255)
                                                                           tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAc                              <  1:520967/71‑1 (MQ=255)
                                                                            caccacCCACCAAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac          <  2:240847/90‑1 (MQ=255)
                                                                            cACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac          >  1:569606/1‑90 (MQ=255)
                                                                            cACCAACCATCCATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGacac          <  1:676075/90‑1 (MQ=255)
                                                                                  ccatcaaaAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg    <  2:48745/90‑1 (MQ=255)
                                                                                    atccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg  >  2:683091/1‑90 (MQ=255)
                                                                                         |                                                                                    
AAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGG  >  NZ_CP009273/4326629‑4326802

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