Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,080,0570TG53.8% 7.1 / 12.7 13V409G (GTG→GGG) efeBiron uptake transporter deferrochelatase/peroxidase subunit
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base G (7/0);  total (8/5)
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.13e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGTTTTTAGGC  >  NZ_CP009273/1079980‑1080114
                                                                              |                                                         
aaaaaGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGGGAAGGACGCGaa                                               >  2:698022/1‑90 (MQ=255)
     ggCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGatt                                          <  1:939151/90‑1 (MQ=255)
     ggCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGGGAAGGACGCGAACGatt                                          >  1:527005/1‑90 (MQ=255)
     ggCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑‑TTTTTTGCGCTGCCGGGGGGGAAGGACGCGAACGatta                                         >  2:326685/1‑90 (MQ=255)
               cGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAg                                <  1:33721/90‑1 (MQ=255)
                        ggaggaATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGGGAAGGACGCGAACGATTATTTCGGAAGCGCGTTAtt                       >  1:118142/1‑90 (MQ=255)
                             aaTACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGTCCGGGGGGGAGGACGCGGACGATTATTTCGGAAGCGCGTTATTGCggg                  >  1:652638/1‑90 (MQ=255)
                             aaTACGTTAAACCTATCGGCGGCGGGTA‑TTTTTTTGCGCTGCCGGGGGGGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCggg                  >  1:451139/1‑90 (MQ=255)
                                   ttAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGc                     >  1:663805/1‑81 (MQ=255)
                                   ttAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGc                     <  2:663805/81‑1 (MQ=255)
                                           aTCGGCGGCGGTTATTTTTTTTGCGCTGCCGGGGGGGAAGGACCCGACCGATTATTTGGGAAGGGCGTTTTTGGGGGTTTAATGTTTTTa     >  2:372007/1‑90 (MQ=255)
                                           aTCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGTTTTTAg    <  2:527005/90‑1 (MQ=255)
                                             gggcggggtTTT‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGTTTTTAGGc  <  1:372007/89‑1 (MQ=255)
                                                                              |                                                         
AAAAAGGCTCAATGGCGAAGCGCTGGAGGAATACGTTAAACCTATCGGCGGCGGTTA‑TTTTTTTGCGCTGCCGGGGGTGAAGGACGCGAACGATTATTTCGGAAGCGCGTTATTGCGGGTTTAATGTTTTTAGGC  >  NZ_CP009273/1079980‑1080114

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