Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3570AT60.0% 9.2 / 32.6 30intergenic (+47/‑56)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base A (2/10);  new base T (18/0);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 2.20e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT  >  NZ_CP009273/2380276‑2380446
                                                                                 |                                                                                          
aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGGCAAGGCCATTTTTTTTTTCCt                                                                                    >  2:71136/1‑90 (MQ=255)
  ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGACAAGCACAATTTTTTTTTCCTcc                                                                                  >  2:512611/1‑89 (MQ=255)
  ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCACAAAAAAAGAAAAACACAATTTTTTTTTCCTcc                                                                                  >  2:1149770/1‑89 (MQ=255)
    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                <  2:900614/90‑1 (MQ=255)
    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                <  2:785688/90‑1 (MQ=255)
    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                <  2:286857/90‑1 (MQ=255)
    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                <  1:877374/90‑1 (MQ=255)
      ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAAACATTTTATTTTTTCTCTTCTg                                                                              >  1:135378/1‑90 (MQ=255)
           aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTa                                                                        <  2:378808/90‑1 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCtttttttta                                                                      >  2:877374/1‑87 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTGtttttaa                                                                      >  2:375838/1‑87 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCCTCTGgttttta                                                                      >  1:467908/1‑87 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCATCtttttttaa                                                                      >  1:93876/1‑87 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAAAACATTTTTTTTTTTCTCATCTGtttttaa                                                                      >  1:726668/1‑87 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTTTTTTTTCTCCTCTGtttttaa                                                                      >  2:505331/1‑87 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGAAAATTTTATTTTTTCTCATCTTTTTTTAg                                                                      >  1:892796/1‑90 (MQ=255)
             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAAAATTCACATTTTTTTTTTTCTCCTCtttttttat                                                                      >  2:419146/1‑87 (MQ=255)
              ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCtttttttaaa                                                                     >  2:209156/1‑86 (MQ=255)
                   ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa                                                                <  1:512611/90‑1 (MQ=255)
                        acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTGTTTTTAAATGCGaaaaaa                                                           >  1:143381/1‑87 (MQ=255)
                        acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCCTCTTTTTTTAAATGAGaaaaaa                                                           >  2:5206/1‑87 (MQ=255)
                           gTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGAACATTTTTTTTTTTCTCCTCTTTTTTTATATGGGAAAAAAGaa                                                        >  1:1045793/1‑90 (MQ=255)
                           gTAAATAATAAAAATTTCTCGAGGCGGGGCCAGGAAAAACGCCAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATGCGGAAAAAGaa                                                        >  2:12632/1‑90 (MQ=255)
                                   taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg                                                <  1:797346/90‑1 (MQ=255)
                                   taaAAATTTCTCGAGGCGGGTACCGGAAAAAAGACAAGAACATTTTTTTTTTCCTCTTCTGTTTTTAAAATGGAAACAAGAAATTaaaag                                                >  1:79874/1‑88 (MQ=255)
                                     aaaaTTTCTCGAGGCGGGGGCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAAATTGGAAAAAAGAAATTTAAAgga                                              >  2:1007607/1‑90 (MQ=255)
                                     aaaaTTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAGATGGGAAACATGAAATTTAACgga                                              >  2:996180/1‑90 (MQ=255)
                                               gAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataa                                    <  2:1045793/90‑1 (MQ=255)
                                                      gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                             <  1:592347/90‑1 (MQ=255)
                                                      gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa                             <  1:419146/90‑1 (MQ=255)
                                                                             ttttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAtt      <  1:186982/90‑1 (MQ=255)
                                                                              tttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt                            <  2:320125/67‑1 (MQ=255)
                                                                              tttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt                            >  1:320125/1‑67 (MQ=255)
                                                                               ttattttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGa    >  2:54164/1‑90 (MQ=255)
                                                                                 attttCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACt  <  2:1085989/90‑1 (MQ=255)
                                                                                 |                                                                                          
ATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACT  >  NZ_CP009273/2380276‑2380446

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