Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273890,3810AC52.4% 21.5 / 25.4 21intergenic (+29/‑66)potF/potGspermidine/putrescine ABC transporter substrate‑binding protein PotF/putrescine ABC transporter ATP‑binding subunit PotG
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/11);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTC  >  NZ_CP009273/890311‑890465
                                                                      |                                                                                    
gtgCGCCCCCGCGCGGGGACCAAAGTGAAGAGCGAAAAATAACCCGCAGTCGTAGATGCCGGGGGGGCGCCCCACACCCGCCGGCAATTc                                                                   <  1:395194/90‑1 (MQ=255)
  gcgcCCCCGCGCGTGGACCAAAGTAAAGAGCGGAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGc                                                                 <  2:113970/90‑1 (MQ=255)
  gcgcCCCCGCGCGGGGACAAAAGTGAAGAGCGGAAAATAACCCGCAGTCGTAGAGGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGc                                                                 <  2:206207/90‑1 (MQ=255)
     cccccGCGCGTGGACCAAAGTTAAGAGCGAAAAATAACCCCCAGTCGAAGAGGCGGGAGGGGCGCCCCACCCCCGCCGGCAATTCGcccc                                                              <  2:591977/88‑1 (MQ=255)
          gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                         >  2:563366/1‑90 (MQ=255)
          gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                         >  2:520333/1‑90 (MQ=255)
          gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                         >  2:25038/1‑90 (MQ=255)
             cgTGGACAAAAGTGAAGAGCGGAAAAAAACCCGCAGTCGTAGATGCCGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGt                                                      <  2:508273/90‑1 (MQ=255)
              gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTg                                                     >  2:116794/1‑90 (MQ=255)
               tGGACAAAAGGGAAGAGGGAAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTgc                                                    <  2:294626/90‑1 (MQ=255)
               gggACAAAAGTGAAGAGCGAAAAATAACCCGCAGTCGTAGAGGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTgc                                                    <  2:608719/89‑1 (MQ=255)
                 gACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgcgc                                                  >  1:655183/1‑90 (MQ=255)
                           aaGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACAt                                        >  1:563229/1‑90 (MQ=255)
                             gggCGAAAAAAAATCCCCAGTCGAAGAGGCGGGGGGGCCGCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTc                                      <  2:655183/88‑1 (MQ=255)
                             gagCGGAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTc                                      <  2:146416/90‑1 (MQ=255)
                                cgaaaaaTAACCCGCAGCCGTAGAGGCGGGGGGGGCCCCCCCCACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAAt                                   <  2:650955/87‑1 (MQ=255)
                                     aataacCCGCAGTGGTAGAGGCGGGGGGGGCCCACCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCgg                              <  2:573007/90‑1 (MQ=255)
                                         atCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGagag                          >  1:425167/1‑90 (MQ=255)
                                          tCCGCAGTCGTAGATGCGGGGGGGGCGCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCATTGCCGGAGagc                         <  1:373200/90‑1 (MQ=255)
                                                        tGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACg           >  1:701468/1‑90 (MQ=255)
                                                                 ggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTc  >  1:628281/1‑90 (MQ=255)
                                                                      |                                                                                    
GTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTC  >  NZ_CP009273/890311‑890465

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