Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273890,3860AC53.8% 11.5 / 14.5 13intergenic (+34/‑61)potF/potGspermidine/putrescine ABC transporter substrate‑binding protein PotF/putrescine ABC transporter ATP‑binding subunit PotG
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base C (0/7);  total (5/8)
Fisher's exact test for biased strand distribution p-value = 4.66e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAG  >  NZ_CP009273/890305‑890473
                                                                                 |                                                                                       
gACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCgg                                                                                 >  1:706847/1‑90 (MQ=255)
        gcgcCCCCGCGGGTGGACAAAAGGGAAGAGCGGAAAATAACCCGCAGTCGTAGAGGCCGGAGGGGCGCACCACCCCCGCCGGCAATTCGc                                                                         <  1:482787/90‑1 (MQ=255)
        gcgcACCGGGGCGAGGAAAAAGGGGAAGGGGGGAAAAAAATCCGAAGTCGTGGATGCGGGGGGGGCGCACCACACCCGCCGGCAATTCGc                                                                         <  1:372646/90‑1 (MQ=255)
          cccccccCGCGTGGACAAAAGTGAAGGGCGGAAAATAATCCGCAGTTTAAGATGCCGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCAc                                                                       <  2:763405/87‑1 (MQ=255)
                gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACCCCCGCCGGCAATTCGCACCATTAt                                                                 <  2:127511/90‑1 (MQ=255)
                gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt                                                                 >  2:365186/1‑90 (MQ=255)
                         caaaaGTGAAGAGCGAAAAATAACCCGCAGTCGTAGATGCGGGAGGGGCCCACCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCtt                                                        <  2:536986/88‑1 (MQ=255)
                                 aaGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACAt                                                >  1:556609/1‑90 (MQ=255)
                                   gagCGGAAAATAATCCGCAGTGGAAGAGGCGGGGGGGGCCCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTc                                              <  1:31241/90‑1 (MQ=255)
                                                cccGCAGTCGTAGATGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagc                                 <  2:556609/89‑1 (MQ=255)
                                                   gCAGTCGTAGATGCCGGAGGGGCGCACCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagcagc                              <  2:539890/90‑1 (MQ=255)
                                                                ccGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCt                 >  2:774732/1‑90 (MQ=255)
                                                                      aagcgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCt           >  2:310536/3‑90 (MQ=255)
                                                                               acacCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAg  >  2:427247/1‑90 (MQ=255)
                                                                                 |                                                                                       
GACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCAG  >  NZ_CP009273/890305‑890473

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