Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,500,2850AC53.8% 12.0 / 16.2 13V18G (GTA→GGA) ydcOBenE family transporter YdcO
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base C (0/7);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAG  >  NZ_CP009273/1500205‑1500364
                                                                                |                                                                                
cagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGaa                                                                         >  2:480044/1‑90 (MQ=255)
                       ccGGCAACAATCGCTGCTTGCCAGATTATTGCCGCTGAACGGGCGTAGCAAATTAATCCCGCCAGAAACCCCGCCAGTAGCGTGGGTgga                                                 <  1:466564/90‑1 (MQ=255)
                        cGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggag                                                >  1:40849/1‑90 (MQ=255)
                                         tGCCAGTTTATGGCCGCTGACCGGGGGTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGc                               <  2:114365/90‑1 (MQ=255)
                                         tGCCAGATTATTGCCGCTGAACGGGCGTAGCCAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGc                               <  1:230374/90‑1 (MQ=255)
                                              gATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGt                          >  2:203127/1‑90 (MQ=255)
                                              gATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGt                          >  2:543172/1‑90 (MQ=255)
                                                   ttGCCGCTGAACGGGCGTAGCCAATTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                     <  1:143224/90‑1 (MQ=255)
                                                     gCCGCTGACGGGGCGTAGCAAATAAAGCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                   <  2:112904/90‑1 (MQ=255)
                                                             aggggCGTACCAAATAAACCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat           <  2:51216/87‑1 (MQ=255)
                                                                ggCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGc        >  2:10433/1‑90 (MQ=255)
                                                                     aGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACa   >  1:560973/1‑90 (MQ=255)
                                                                       ccAAATAAACCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAg  <  2:189042/90‑1 (MQ=255)
                                                                                |                                                                                
CAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAG  >  NZ_CP009273/1500205‑1500364

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