Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA61.9% ‑6.3 / 41.3 42intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/16);  new base A (26/0);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 6.01e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-08

GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCT  >  NC_000913/2386725‑2387008
                                                                                                                                         |                                                                                                                                                  
gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACATACACatt                                                                                                                                       >  1:601005/1‑151 (MQ=255)
gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGAAAAGCACatt                                                                                                                                       >  2:469574/1‑151 (MQ=255)
gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAACACatt                                                                                                                                       >  1:303870/1‑151 (MQ=255)
     aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatt                                                                                                                                  <  2:328532/151‑1 (MQ=255)
     aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACattttttt                                                                                                                                  >  2:262313/1‑151 (MQ=255)
      aaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttattt                                                                                                                                 <  2:495401/151‑1 (MQ=255)
        cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACattttattttt                                                                                                                               >  2:214243/1‑150 (MQ=255)
        cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAACTCaatttattttt                                                                                                                               >  1:244666/1‑150 (MQ=255)
          aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAAGACattttttttttc                                                                                                                              >  2:596742/1‑148 (MQ=255)
          aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                                                                             <  1:360109/151‑1 (MQ=255)
                    cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGttt                                                                                                                   >  2:390668/1‑151 (MQ=255)
                         tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                              <  1:214243/151‑1 (MQ=255)
                         tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                              <  1:337715/151‑1 (MQ=255)
                         tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                              <  2:121584/151‑1 (MQ=255)
                             aCCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCCCGAGGCGGGTCCAGAAAAAACGACAAGCCCATTTTATTTTTCCTCATCTGGTTTTAGATGCg                                                                                                          >  2:136087/1‑151 (MQ=255)
                                             aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                                                          <  2:103016/151‑1 (MQ=255)
                                                aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTTTTTTTAAATGCGGAAAAAGAAAATTAACggg                                                                                       >  2:427451/1‑150 (MQ=255)
                                                aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACCTTTTTTTTTTCCTCCTCTGTTTTTAGATGCGGAAAATTAAAATCAACggg                                                                                       >  2:343542/1‑150 (MQ=255)
                                                             tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATg                                                                          <  1:596736/151‑1 (MQ=255)
                                                                ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAATCACAATTTATTTTTCCTCATCTTTTTTTAGAATCGAAAAAAGAAAATTAAAGGGGGAAAAAAAAATtggg                                                                       >  1:525545/1‑150 (MQ=255)
                                                                    gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTATTTTTCCTCATCTGTTTTTAAATGCGGAAAAAGAAAATCAACGGAGGGAAAATAAATGGGAATGa                                                                   >  2:348290/1‑151 (MQ=255)
                                                                           gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCATCTGTTTTTAGAAGCGAAAAATGAAATTTAACGGAGGAAAAATAAATTGGGATTATTggtgt                                                            >  1:659890/1‑148 (MQ=255)
                                                                           gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCACATTTTATTTTTCCCCCTCTGTTTTTAAATGCGAAAAAAGAAAATTAAAGGAGGGAAAAAAAAAGGGGATGATTGGcgt                                                            >  1:782730/1‑149 (MQ=255)
                                                                           gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGCACATTTTTTTTTTCCTCCTTTGTTTTTAAATTCGGAAAATGAAATTTAAAGGGGGGAAAATAAATAGGGATGATTGGCTa                                                            >  2:252570/1‑151 (MQ=255)
                                                                           gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGCACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAACAAGAAAATTAAAGGGGGAAAAATAAATTGGGATGATTTGCtt                                                            >  2:396107/1‑150 (MQ=255)
                                                                            tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACATGAAAATTTTATTTTTTCTCCTCTGTTTTTAAAAGCGGAACAAGAAAATTTACGGAAGAAAAATAAATTGGGATGATTGGCTAt                                                           >  1:214831/1‑151 (MQ=255)
                                                                             gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCACAAAAAACGGAACACACAATTTTTTTTTCCACTTCTTTTTTTAAATGCGAAAAAAAAAATTTAAAGGAGGGAAAATAAATGGGGATGGTTggttttt                                                          >  2:43896/1‑146 (MQ=255)
                                                                              ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAttt                                                         <  2:413849/151‑1 (MQ=255)
                                                                                   aCTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAg                                                    <  2:600442/151‑1 (MQ=255)
                                                                                        gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTATTTTTCCTCATCTTTTTTTAGATGCGAAACAAGAAATTCAACGGAGGGAAAAAAAATTGGGATGATTGGCTATTTTGCAGAAAtt                                               >  2:643602/1‑151 (MQ=255)
                                                                                           ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGAAGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTTGCTATTTTGCAGAAAATGaa                                            >  2:377396/1‑150 (MQ=255)
                                                                                              caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGGACATTTTATTTTTCCTCATCTGTTTTTAGATGCGAAACATGAAATTTAACGGAGGAAAAATAAATGGGGATGATTGGCTATTTTTCAGAAATTGAATCa                                         >  1:811837/1‑151 (MQ=255)
                                                                                              caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAAAAAAATTTATTTTTTCTTATCTGGTTTTAGAAGCGAAACAAGAAAATTAACAGAGGAAAAATAAATGGGGATGAATTGGTATTTTTCAGAAATTGaaata                                         >  1:74094/1‑147 (MQ=255)
                                                                                              caacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACATTTTATTTTTTCTCATCTGTTTTTAGATGCGAAAAATGAAATTCAAAGGAGGGAAAATAAATGGGGATGATTGGCCATTTTTCAGAAAATGACTCa                                         >  1:57054/1‑151 (MQ=255)
                                                                                                acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAAAAAGAAATTCAACGGAGGAAAAAAAAATGGGGATGATTGGCTATTTTGTAGAAAATGACTCAGa                                       >  2:398884/1‑151 (MQ=255)
                                                                                                       ataataAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTATTTTTTCTCATCTGTTTTTAGATGCGAAAAATGAAATTTAACGGGGGAAAAATAAATGGGGATGATTGGCTATTTTTTAGAAAAAGACTCAGaaaaaaat                                >  1:215123/1‑147 (MQ=255)
                                                                                                                      cGAGGCTGGTCCAGACAAAGCGACATGCCCATTTCATTTTCCCTCATCTGTTTTAAGATGCTAAACATGACATTCACCTGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTACTTGAg                 <  1:186553/151‑1 (MQ=255)
                                                                                                                               tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa         <  2:399282/150‑1 (MQ=255)
                                                                                                                                    aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAcc   <  2:36796/151‑1 (MQ=255)
                                                                                                                                    aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc    <  1:54184/150‑1 (MQ=255)
                                                                                                                                    aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAGTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAATGAAAAAc    <  1:416053/150‑1 (MQ=255)
                                                                                                                                     aaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCt  <  2:20196/151‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCT  >  NC_000913/2386725‑2387008

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