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

TTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTG  >  NC_000913/2386641‑2387092
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
ttCCGGGTATTGATTGCTCGGCCCCCCTGATTTATTGAAAACCCGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACattttat                                                                                                                                                                                                                          <  1:452633/239‑1 (MQ=255)
                                       aaCACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCTAAAAAAAGGAAATGCACATTTATTTTTCTCTCTACTGTTTTTAGAATGGAAACATTAAAtta                                                                                                                                                                                    >  2:226582/1‑238 (MQ=255)
                                                ttCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGaa                                                                                                                                                                         <  1:226582/239‑1 (MQ=255)
                                                        aTAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAaaataaat                                                                                                                                                                 <  1:161788/239‑1 (MQ=255)
                                                                                     taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATg                                                                                                                                                          <  1:503636/217‑1 (MQ=255)
                                                                                     taataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAACGACAAGTACAATTTTTTTTTTCTCATCTGTTTTTAAAAGAGGAAAATTAAAAATAACAGA‑GGGAAAAAAAATAGGGAAGAG‑TGGGGATTTTgttgcaaatg                                                                                                                                   >  2:261393/1‑230 (MQ=255)
                                                                                             aaTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGaa                                                                                                                           <  1:505926/239‑1 (MQ=255)
                                                                                                        cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCAt                                                                                                                                                                                                                <  1:22942/145‑1 (MQ=255)
                                                                                                                       gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGATCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGaaaa                                                                                                                         <  2:488508/215‑1 (MQ=255)
                                                                                                                       gacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGATCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGaaaa                                                                                                                         >  1:488508/1‑215 (MQ=255)
                                                                                                                            aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACattttatt                                                                                                                                                                                                                         <  1:162983/116‑1 (MQ=255)
                                                                                                                            aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgga‑gg                                                                                                                                                                           <  1:145597/161‑1 (MQ=255)
                                                                                                                                    aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTGGTTATACAACAGTAAATAATAGAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTg                                                                                   <  2:415861/239‑1 (MQ=255)
                                                                                                                                      cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAAAGAAAAGCCCATTTTTTTTTTCCCCATCTGTTTTTAGAATGGAAAAAATAAAATCAACAGA‑GGGaaaaaaaaa                                                                                                                                                                 >  2:145597/1‑160 (MQ=255)
                                                                                                                                          aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTTCCGAAAAAAAGACAAGCACAATTTATTTTTTCTCATCTGTTTTTAGATGTGAAAAATGAAAATCAACGGAGGGAAAAAAAAATGGGATAGATGGGGCTTTTTTAGA‑AAATTAACTAGAAAAAAATTAACAACTTGCTTGGGT‑TAC‑TGGAAAAcacctttttggt                                                                             >  2:547351/1‑228 (MQ=255)
                                                                                                                                           aaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGata                                                                            <  1:338935/239‑1 (MQ=255)
                                                                                                                                               gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTTCCGAAAAAAAGACATGCACATTTTATTTTTCCTCCTCTGTTTTTAAATGCGAAACATGAAAATTAAACGA‑GGGAAAATAAAAGGGGATGAT‑TGTGTATTTTTTCGAAAATGACTCAGAAAAAAATAAAACACAGCCTGAGTGTACTTTAAAAACCtttttatggaatat                                                                         >  2:624127/1‑228 (MQ=255)
                                                                                                                                                              ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTc                                                         <  2:450301/239‑1 (MQ=255)
                                                                                                                                                                tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCCGAAAAAAAGACATGGACATTTTATTTTTTCTCATCTGTTTTTATATGCGAAACATGAAAATCAACGGA‑GGGAAAAAAAATGGGGATGTT‑TTGCTATTTTTCAGAAAATTAATCAGAAAAGATTAAACAAATGATTGTGTGTAC‑TGAAAAACCATTGATGGATAATATACATAAAACGCTTTCac                                                       >  2:490246/1‑238 (MQ=255)
                                                                                                                                                                tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGa                                                                              <  1:658251/216‑1 (MQ=255)
                                                                                                                                                                                taCAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgga‑gg                                                                                                                                                                           <  2:168519/109‑1 (MQ=255)
                                                                                                                                                                                          aataatAAAAATTTCTCGAGGC‑GGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAAAATGAAAATTAAAGGA‑GGGaaaaaaaaa                                                                                                                                                                 >  1:168519/1‑108 (MQ=255)
                                                                                                                                                                                                                  ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTg       <  1:604020/238‑1 (MQ=255)
                                                                                                                                                                                                                    tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTggg     <  1:490246/238‑1 (MQ=255)
                                                                                                                                                                                                                          aaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTg  <  1:261393/235‑1 (MQ=255)
                                                                                                                                                                                                                              gCGACATGCACATTTTATCTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                                                                                                               <  1:709659/123‑1 (MQ=255)
                                                                                                                                                                                                                              gCGACATGCACATTTTATCTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCaa                                                                                                               >  2:709659/1‑123 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
TTCCCGGTATTGATTTCTCTGACCACCTGAATTATTGGCAACACGATATTCCCGCCATAATGATTACTGACACCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGC‑GGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGA‑GGAAAAATAAATGGGAATGAT‑TGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAC‑TGAAAAACCTTTGATGGATAATATACATGACACGCTTTCAGGATTACGCCGACTGGATATTGATAAACGCTGGGATTTTTTACACTTTGGGCTG  >  NC_000913/2386641‑2387092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑