Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,291,0630AC62.5% 3.6 / 12.4 16G145G (GGT→GGGyraQpermease
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base C (0/10);  total (5/11)
Fisher's exact test for biased strand distribution p-value = 1.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGCCGCAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACAC  >  NZ_CP009273/3290974‑3291126
                                                                                          |                                                               
aTCGCCGCAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCa                                                                 >  1:238504/1‑90 (MQ=255)
     cgcAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGcc                                                            <  2:525044/90‑1 (MQ=255)
                 gCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTg                                                >  2:283395/1‑90 (MQ=255)
                           aaGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGc                                      >  1:497091/1‑90 (MQ=255)
                                 aTAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTccc                                >  1:261167/1‑90 (MQ=255)
                                          aGCGTCGCCGGGTTTAACCCCGGATTGCCCATCCAGAGTGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGAcc                       <  1:370421/90‑1 (MQ=255)
                                          aGCGTCGCCGGGTTTAACCCCGGATTGCCCATCCAGAATGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGAcc                       <  1:91796/90‑1 (MQ=255)
                                           gCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCg                      >  2:617882/1‑90 (MQ=255)
                                               ccccGGGTTTAACCCCGGATTGCCCACCCAAAATGCCCGCCCCCCCCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGc                  <  1:548143/88‑1 (MQ=255)
                                                 gcGGGGTTAACCCCGGAATCCCCCACCAAGAAGGCCGCCCCCCCCCCCCCCGACCCCTTTTGGCGACGCTTCCCCGCCGCGACCGGAGcc                 <  2:456734/89‑1 (MQ=255)
                                                   ggttttAACCCCGGATCCCCCACCCAGAAGGCCAGCCCCCCCCCCATCGACCCCTGTTGGCGCCGCATTCCCGCCGCGACCGGAGCCGca              <  2:497091/87‑1 (MQ=255)
                                                     ttttACCCCCGGTTTGCCCACCAAGAATGCCAGCCCCCCCCCAATCGACACCTGTTGGCGACGATTTCCCGCCGCGACCGGAGCCGcaca            <  1:617882/89‑1 (MQ=255)
                                                     gTTTACCACCGGATTGCCCATCCAGAATGCCAGCCCCCCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGcaca            <  1:498698/90‑1 (MQ=255)
                                                       ttAACACCGGATTGCCCACCCAGAATGCCAGCGCCCCGCCCATCGACCCCTTTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCAcagc          <  1:323781/90‑1 (MQ=255)
                                                             ccGGTTTGCCCATCCAGAATGCCAGCGCCCCGCCCATCGACCCCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTac    <  1:422603/90‑1 (MQ=255)
                                                             ccGGTTTCCCCACCCAGAATCCCCGCCCCCCCCCCACCGACCCCTGTTGGCGCCGCTTTCCCGCCGCGACCGGAGCCGCACAGCAGGTac    <  2:330410/90‑1 (MQ=255)
                                                               ggATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTacac  >  2:228834/1‑90 (MQ=255)
                                                                                          |                                                               
ATCGCCGCAAAACCCCAGCCGAGGACAAAGCCCATAAACACCAGCGTCGCCGGGTTTAACACCGGATTGCCCATCCAGAATGCCAGCGCACCGCCCATCGACACCTGTTGGCGACGCATTCCCGCCGCGACCGGAGCCGCACAGCAGGTACAC  >  NZ_CP009273/3290974‑3291126

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