Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA30.8% 38.8 / 11.4 26intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (1/17);  new base A (8/0);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 5.76e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCCCAGGG‑TTTATTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACT  >  NC_000913/2386628‑2387084
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                               
aCCCCAGGGTTTTATCCCGGGTTTTGATTCCTCTGACCACCTAAATTATGGAAAACACAATTTTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGa                                                                                                                                                                                                                              <  1:455203/239‑1 (MQ=255)
     aGGG‑TTTATTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGCGTTATAACTTTTTTATACAACAGTAAATAATAAAAATTTCTCCAGGCGGGTTCACAAAAAACGagaact                                                                                                                                                                                                                         >  1:460716/1‑233 (MQ=255)
                                         tGATTTATGGGAAACACAATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaa                                                                                                                                                                                     <  1:570670/239‑1 (MQ=255)
                                                                       tAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCg                                                                                                                                                                                       <  1:81851/207‑1 (MQ=255)
                                                                                                   taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                                                                                                                                                                         <  1:394063/145‑1 (MQ=255)
                                                                                                           aaTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGATATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACat                                                                                                                                                                                                                     <  1:816394/141‑1 (MQ=255)
                                                                                                                     gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGAAAAGAACaatttttttttt                                                                                                                                                                                                           >  2:816394/1‑139 (MQ=255)
                                                                                                                     gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAACGAAAAGCACAATTTATTTTTTCTCATCTGTTTTTAGATGTGAAACAAA‑AAATTCAACAGGGGGAAAATAAAAAGGAAGGGTTGGCTATTTTTTAGAAATTAAATCAGAAAAAATTaaaa                                                                                                        >  2:352727/1‑237 (MQ=255)
                                                                                                                              tACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttat                                                                                                                                                                                                                <  1:337697/127‑1 (MQ=255)
                                                                                                                                          aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATg                                                                                                                                                                                                                          <  1:195392/105‑1 (MQ=255)
                                                                                                                                                        aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                                                                         <  1:17422/235‑1 (MQ=255)
                                                                                                                                                        aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                                                                         <  1:745439/235‑1 (MQ=255)
                                                                                                                                                        aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGAt                                                                         <  1:352727/235‑1 (MQ=255)
                                                                                                                                                                                         tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAaaata                                                                                                                                                          <  2:370293/121‑1 (MQ=255)
                                                                                                                                                                                         tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCt                                                                                                 <  1:795175/178‑1 (MQ=255)
                                                                                                                                                                                            tataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACt                                  <  2:545655/238‑1 (MQ=255)
                                                                                                                                                                                                   caGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa                                                                                   <  2:384325/182‑1 (MQ=255)
                                                                                                                                                                                                   caGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACATGca                                                                                                                                                                                                                        >  1:370293/1‑50 (MQ=255)
                                                                                                                                                                                                   caGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGCACATTTTTTTTTTTCTCTTCTTTTTTTATATTCGAAAAAAG‑AAAATTAAAGGAGGAAAAAAAAAATGGAATTATTTGCTATTTTTTAGAAAATTACTCAAAAAAAATTAAACAAATGCTTTAGTCttat                                                                                       >  2:795175/1‑175 (MQ=255)
                                                                                                                                                                                                     gTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAAATGGGAAAAAAG‑AAAATTAAAGGAGGGAAAAAAAAAGGGGAAGATTGGGTATTTTTTAGAAAATGAAACAGAAAAAAATAAACAAATGCTTGAGTTTAATGAAAAAACTTTTATGGATAAAATAAAATAAAAGCTTTTAGGAATAAACCGACTgggaatttat                        >  1:545655/1‑231 (MQ=255)
                                                                                                                                                                                                      tAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGata                             <  1:420129/233‑1 (MQ=255)
                                                                                                                                                                                                             taaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAGATGCGAAAAAAG‑AAAATTAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTTCAGAAAATGAATCAGAAAAAAATAAACAACTGCTTGAGTCTACTGAAAAAACTTTTAt                                                                         >  1:384325/1‑182 (MQ=255)
                                                                                                                                                                                                                  aTTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGGACAATTTTTTTTTCCTCCACTGTTTTTAAAAGGGAAACAAGAAAATTAAAGGGGGGAAAAAAAAAAGGGAA‑GATTGGCTATTTTTGAGAAAATTAATCCGAAAAAAATAAAAAAATGCTTTAGGTTACTGAAAAAACTTTTATGGATAATATAAAATAAACGCTTTCAGGATTACGCCGACTGGATAATTAAAAAAGCTGGGAtt           >  1:863394/1‑239 (MQ=255)
                                                                                                                                                                                                                             cGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACt  <  1:202497/237‑1 (MQ=255)
                                                                                                                                                                                                                                          aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                                                                                                      <  2:772416/124‑1 (MQ=255)
                                                                                                                                                                                                                                          aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                                                                                                      >  1:772416/1‑124 (MQ=255)
                                                                                                                                                                                                                                           gCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGAtt                                                                                                                                             <  2:385023/84‑1 (MQ=255)
                                                                                                                                                                                                                                           gCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGAtt                                                                                                                                             >  1:385023/1‑84 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                               
ACCCCAGGG‑TTTATTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACT  >  NC_000913/2386628‑2387084

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