Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273575,6140AG63.2% 6.5 / 17.7 19intergenic (+72/+93)BW25113_RS02885/ylcIhypothetical protein/YlcI/YnfO family protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (7/5);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCG  >  NZ_CP009273/575528‑575700
                                                                                      |                                                                                      
cTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAg                                                                                     >  1:124987/1‑90 (MQ=255)
  ggtggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTg                                                                                   >  2:35212/1‑90 (MQ=255)
     ggtTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTa                                                                                >  2:264378/1‑90 (MQ=255)
               tCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAAc                                                                      <  2:124158/90‑1 (MQ=255)
                                  ttCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCt                                                   <  1:18501/90‑1 (MQ=255)
                                        ctactTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAAt                                             >  1:191738/1‑90 (MQ=255)
                                              aaCCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                                       <  1:114683/90‑1 (MQ=255)
                                                ccTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGaa                                     >  2:224791/1‑90 (MQ=255)
                                                     gACAGAGCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                                <  2:316087/90‑1 (MQ=255)
                                                     gACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTAc                                <  1:35212/90‑1 (MQ=255)
                                                         gagCCAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAata                            >  2:211739/1‑90 (MQ=255)
                                                             cAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                                       >  2:107497/1‑75 (MQ=255)
                                                             cAACACGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTg                                       <  1:107497/75‑1 (MQ=255)
                                                               acacGCAGTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTcc                      >  2:215/1‑90 (MQ=255)
                                                                      gTCTGTCACTGTCAGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCgg               <  2:191738/90‑1 (MQ=255)
                                                                            cactgtcaGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGgaa         <  1:215/90‑1 (MQ=255)
                                                                                gtcaGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAc     >  1:205188/1‑90 (MQ=255)
                                                                                  caGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAcgc   >  2:162932/1‑90 (MQ=255)
                                                                                   aGGGAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAAcgcg  <  2:205188/90‑1 (MQ=255)
                                                                                      |                                                                                      
CTGGTGGTTGACTGATCACCATAACTGCTAATCATTCAAACTACTTAACCTGTGACAGAGCCAACACGCAGTCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATTAAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCG  >  NZ_CP009273/575528‑575700

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