Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3620CT52.9% ‑3.4 / 17.9 17intergenic (+52/‑51)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base C (1/7);  new base T (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTC  >  NZ_CP009273/2380275‑2380447
                                                                                       |                                                                                      
acaggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                                      <  1:1249905/87‑1 (MQ=255)
   ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTcc                                                                                   >  1:465528/1‑89 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                 <  1:1088514/90‑1 (MQ=255)
     tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                 <  2:572957/90‑1 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGGAAAAACGGCATGCACATTTTATTTTTCCTCCTCTTGTTTTag                                                                        >  1:93444/1‑89 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTTCTCATCtttttttta                                                                       >  1:161136/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACAATTTATTTTTCCTCCTCTGtttttaa                                                                       >  1:182097/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTCCTCtttttttaa                                                                       >  2:255106/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTATTTTTCCTCCTCTGtttttaa                                                                       >  1:562353/1‑87 (MQ=255)
              tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAAAAATTTTTTTTTTCTCATCTGttttttt                                                                       >  2:723996/1‑87 (MQ=255)
                         acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGaaaaaa                                                            >  1:485571/1‑87 (MQ=255)
                            gTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGAAAAATGaa                                                         >  1:1208035/1‑90 (MQ=255)
                                   ataaAAATTTCTCGAGGCGGATCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                  >  1:615948/1‑90 (MQ=255)
                                                             aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTg                        <  1:629494/90‑1 (MQ=255)
                                                               aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc                      <  1:723996/90‑1 (MQ=255)
                                                                       aTGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGc              <  2:1549343/90‑1 (MQ=255)
                                                                                  attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt   <  2:768125/90‑1 (MQ=255)
                                                                                   ttttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc  <  1:686432/90‑1 (MQ=255)
                                                                                       |                                                                                      
GATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTC  >  NZ_CP009273/2380275‑2380447

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