Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA61.0% 12.9 / 63.5 59intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (5/18);  new base A (36/0);  total (41/18)
Fisher's exact test for biased strand distribution p-value = 5.20e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-10

ATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTA  >  NC_000913/2386728‑2386997
                                                                                                                                       |                                                                                                                                               
ataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatttta                                                                                                                              <  1:716743/151‑1 (MQ=255)
  aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAACACaatttttt                                                                                                                            >  1:475058/1‑151 (MQ=255)
    aCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatttt                                                                                                                          <  2:654771/151‑1 (MQ=255)
     cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATg                                                                                                                                         <  1:415682/137‑1 (MQ=255)
     cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACattttatt                                                                                                                            <  1:888673/148‑1 (MQ=255)
     cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAATCACatttttttttt                                                                                                                         >  2:953640/1‑150 (MQ=255)
     cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCCCGGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTACAAAAAAAAGGAAATGCACattttattttt                                                                                                                         >  2:888980/1‑150 (MQ=255)
       aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAACaatttatttttt                                                                                                                        >  2:154439/1‑148 (MQ=255)
             cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGGACCATTTTTTTTTCCTCAt                                                                                                                   >  1:246654/1‑149 (MQ=255)
                   cGGTGATACCGCAGACAGATTGAATTATCAGAAAATTGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAGGCAAATTTTATTTTCCCTCATCTGTTTTa                                                                                                        <  2:7836/151‑1 (MQ=255)
                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCTGAAAAAATGACAAGCACCCTTTTTTTTTCCTCTTCTCTTTTaaa                                                                                                      >  2:266929/1‑150 (MQ=255)
                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGCACCATTTTTTTTTCCTCCTCTGTTTTTAg                                                                                                      >  2:300641/1‑151 (MQ=255)
                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCAAAATTTTTTTTTCCTCATCTTTTTTaaa                                                                                                      >  2:391661/1‑150 (MQ=255)
                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATAAAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGAAAAGCACATTTTATTTTTTCTCATCTGTTTTata                                                                                                      >  1:761550/1‑149 (MQ=255)
                     gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATAAAACAGTAAATAATAAAAATTTCTCGAGGCGGGATCCTAAAAAAAGAGATGTACAATTTTTTTTTCCttttttttttttaa                                                                                                         >  2:415240/1‑138 (MQ=255)
                               agacagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAAAAAAAAAGAAAAGAACATTTTTTTTTTCCTCCTCTTTTTTTAGAAGCGgaaaa                                                                                           >  2:80135/1‑149 (MQ=255)
                                 acagaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAAc                                                                                            <  2:918640/148‑1 (MQ=255)
                                     aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAt                                                                                     <  1:241833/151‑1 (MQ=255)
                                            tATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgg                                                                              <  2:8600/151‑1 (MQ=255)
                                                      tGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaata                                                                   <  1:501637/151‑1 (MQ=255)
                                                       ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCAAATTTTTTTTTTTCTCCTCTGGTTTTAGATGAGAAAAAAGAAAATTAAAAGAGGGAaaaaaa                                                                  >  1:313685/1‑151 (MQ=255)
                                                        gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaataa                                                                  <  2:278462/150‑1 (MQ=255)
                                                           cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAAACACATTTTTTTTTTTCTCCTCTGTTTTTTTATTCGGAAAATGAAAATTAACGGAGGAAaaataaaatg                                                              >  1:821706/1‑148 (MQ=255)
                                                            aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAACGAAAAGCAAATATTTTTTTTCCTCACCTGTTTTTAGAAGGTTAACATGAAAATCTAACGAGGAaaaaaaaaaagg                                                             >  2:918759/1‑147 (MQ=255)
                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAAAAAAGCACACTTTTTTTTTCCTCTTCTGTTTTTAAAAAGGGAAAAAGAAAATTAACAGAGGGAaaataaaatggg                                                            >  1:855056/1‑146 (MQ=255)
                                                             ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGACAAGCAACATTTTTTTTTCCTTATCCGGTTTTAGATGCGGAAAAAAAAATTTAAAGGGGGGAAAAAAAAATggg                                                             >  1:28333/1‑151 (MQ=255)
                                                              gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGCACAATTTATTTTTTCTCATCTGTTTTAAGATGCGAAACAAGAAAATCAACAGAGGGAAAATAAATGGGaa                                                           >  1:691110/1‑151 (MQ=255)
                                                                 gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAAGAAAATTTTATTTTTCCTCATCTGTTTTTAAATGAGAAAAAAAAAATTCAACAGGGGGAAAATAAAAGGGGATGa                                                        >  1:652822/1‑151 (MQ=255)
                                                                   ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGACAAACAAATTTTATTTTTTCTCCTCTGTTTTTAAATGAGAAAAAAAAAAATTAAAGGAGGGAAAATAAATAGGGAAGAtt                                                      >  1:305600/1‑151 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAACGACAAGCCCAATTTTTTTTTCCCCCTCTTTTTTTAAAATGGAAAAAAGAAAATTAAAAGAGGGAAAAAAAAATGGGAATGTTTggtt                                                  >  2:166650/1‑149 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAAGAAAATTTTTTTTTTCCTCATCTGTTTTTAGAAGCGAAAAAAGAAAATTAAACGAGGGAAAATAAATTGGGATGAATGGCTa                                                 >  1:403992/1‑151 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAAAACACAATTTTTTTTTCCCCTTCTTTTTTTAAAATCGGAAAATAAAAATTAAACGAGGGAAAAAAAATTGGGAATgttgggta                                                  >  2:1042354/1‑144 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTTTTTTTAAAAACGAAAAAAGAAAATTAAAGGGGGAAAAAAAAAAAGGAATGATTGGGTa                                                 >  1:344950/1‑151 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                  <  2:424895/150‑1 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                  <  1:673827/150‑1 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGAACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGTAACAAAAAAATTAAAGGAGGAAAAATAAATAGGGATGATTTGCtt                                                 >  2:329827/1‑150 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTACAGAAAAAAAGACATGCACAATTTATTTTTCCCCATCTTTTTTTAGAAGCGAAAAAAGAAAATTAACGGAGGAAAAATAAATGGGAATTATTGGCTa                                                 >  2:363618/1‑151 (MQ=255)
                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAAACAAGCACATTTTTTTTTTCCCCCTCTGTTTTTAGATGCGGAAAAAAAAATTAAAAGGCGGGAAAATAAATTGGAAAGaattggta                                                 >  1:546963/1‑144 (MQ=255)
                                                                         tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa                                                 <  1:920445/150‑1 (MQ=255)
                                                                          gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAAACAAAAAGAAAATTTTTTTTTTCCCATCTGTTTTTAAAAGCGGAAAAATAAAATTAAACGGGGGAAAAATAAAGGGCAAAAATTGGcatt                                                >  1:591579/1‑148 (MQ=255)
                                                                                   ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCCAAAAAAGAACAACCCAATTTTTTTTTTCCTCCTCTTGTTTTAGATGGGGAACATGAAAATTAACGGGGGGAAAATAAAATGGGATGGATTGGTTTTTTTCAGaaa                                      >  2:511788/1‑151 (MQ=255)
                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGGACAATTTTTTTTTCCTCCTTTTTTTTTAAAAGGGGAAAAAGAAAATAAACGGAGGAAAAAAAAAATGGGATTATTGGGTATTTTTTAGaaaat                                    >  2:544933/1‑149 (MQ=255)
                                                                                      ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGAAGCGAAAAATGAAAATCAAAGGAGGGAAAAAAAAAAGGAATGATTTGCTATTTTTTAGAAAATg                                   >  2:116850/1‑151 (MQ=255)
                                                                                        ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGGGAAGCAAAATTTATTTTTCCTCATCTCGTTTTAAATGGGGAAAAAGAAAATTAAAAGGGGGAAAAAAAATAGGAATAATTGGTTCTTTTTCAGaaaataaa                                 >  2:129691/1‑146 (MQ=255)
                                                                                        ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc                                 <  2:30262/151‑1 (MQ=255)
                                                                                             acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCAAATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAAAAAGAAAATCAAAGGAGGGAAAAAAAATAGGGATGATTTGCTATTTTTGAGAAAATGAATCAGa                            >  1:857572/1‑151 (MQ=255)
                                                                                             acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGGAAATTTTATTTTTCCTCATCTGTTTTTAGAAGCGAAAAAAGAAATTTAACCGGGGAAAAATAAAAGGGGATGATTTGGTATTTTTTAGAAAATTAATCAGa                            >  1:822664/1‑151 (MQ=255)
                                                                                             acaGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATTGGAAACATGAAAATTAAAAGAGGGAAAATAAATAGGGATGATTTGCTATTTTTGAGAAAATGACTCAGa                            >  1:930860/1‑151 (MQ=255)
                                                                                                gTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAACGGCAAGGCCCTTTTTTTTTTCCCCCTTTTTTTTTAAATGTGGAACAAAAAAATTAACGGAGGAAAAAAAAATTGGAATGATTTTTTTTTTTTCAAAAAATGAATCaaaaaa                         >  2:349249/1‑151 (MQ=255)
                                                                                                     taataaAAATTTCTCGAGGCGGGTCCCGGAAAAAAGACATTGACAATTTTTTTTTCCTCCTCCGTTTTTAAAAGCGAAAAAAGAAAATTAAAGGGGGGAAAAAAAATGGGGAAGATTTGGTTTTTTTTCGAAAATTACTTAGaaaaaaata                    >  2:645978/1‑146 (MQ=255)
                                                                                                        taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAAcc                 <  1:154403/151‑1 (MQ=255)
                                                                                                         aaaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCa                <  1:895699/151‑1 (MQ=255)
                                                                                                                     aGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGt     <  2:857863/150‑1 (MQ=255)
                                                                                                                       gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc    <  2:657541/149‑1 (MQ=255)
                                                                                                                          ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  2:837266/1‑148 (MQ=255)
                                                                                                                          ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  2:88652/1‑148 (MQ=255)
                                                                                                                          ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  1:905706/1‑148 (MQ=255)
                                                                                                                          ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  1:414426/1‑148 (MQ=255)
                                                                                                                          ggTCCAGAAAAAGCGACATGCACATTTTATTTCCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTa  >  2:785704/1‑148 (MQ=255)
                                                                                                                                       |                                                                                                                                               
ATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTA  >  NC_000913/2386728‑2386997

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