Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,144,2950AC36.0% 374.6 / 82.8 198G357G (GGT→GGGmrpantiporter inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (108/18);  new base C (0/71);  total (108/90)
Fisher's exact test for biased strand distribution p-value = 6.11e-37
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 226 total aligned reads displayed.

CGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGACGG  >  REL606/2144220‑2144370
                                                                            |                                                                           
cGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCa                                                                             <  1:5800952/76‑1 (MQ=255)
cGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCa                                                                             >  2:2236662/1‑76 (MQ=255)
cGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCa                                                                             >  1:2457134/1‑76 (MQ=255)
 gATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAc                                                                            >  2:3194502/1‑76 (MQ=255)
  aTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTGAcc                                                                           >  2:7730046/1‑76 (MQ=255)
  aTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAcc                                                                           >  2:3608485/1‑76 (MQ=255)
  aTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAca                                                                           <  2:7465069/76‑2 (MQ=255)
    gTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCt                                                                         <  2:2626338/76‑1 (MQ=255)
    gTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCt                                                                         >  1:469752/1‑76 (MQ=255)
     tATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTg                                                                        <  1:1679701/76‑1 (MQ=255)
     tATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAACGACTTCACCCTg                                                                        <  1:5203100/76‑1 (MQ=255)
      aTGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGACTGACTTCACCCCGc                                                                       >  2:2598428/1‑76 (MQ=255)
      aTGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGc                                                                       >  1:7219580/1‑76 (MQ=255)
      aTGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGc                                                                       >  1:5000812/1‑76 (MQ=255)
      aTGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGc                                                                       >  2:3638169/1‑76 (MQ=255)
        gaagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                                     >  1:2057487/1‑76 (MQ=255)
        gaagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                                     <  2:4744667/76‑1 (MQ=255)
        gaagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                                     >  2:561730/1‑76 (MQ=255)
        gaagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                                     <  1:3882035/76‑1 (MQ=255)
        gaagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                                     <  2:2824144/76‑1 (MQ=255)
         aagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                                    >  2:193221/1‑76 (MQ=255)
         aagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                                    >  2:5495235/1‑76 (MQ=255)
         aagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                                    <  2:6341201/76‑1 (MQ=255)
            aaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTag                                                                 <  2:4460404/76‑1 (MQ=255)
            aaTAACCAGAATGGTTAATAGGCGTGAGACCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTag                                                                 <  2:3721776/76‑1 (MQ=255)
            aaTAAACAGAATGGTAAATAGGCGTTATCCCGCGGGGAAGGAACTCTCGCCTGGAATGACTTCCCCCTGCCAGTag                                                                 <  2:5704978/76‑1 (MQ=255)
             aTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTaga                                                                <  1:1708089/76‑1 (MQ=255)
             aTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTaga                                                                >  2:1495067/1‑76 (MQ=255)
              tAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTagag                                                               >  1:5764465/1‑76 (MQ=255)
               aaCCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagc                                                              >  2:4758824/1‑76 (MQ=255)
               aaCCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagc                                                              >  1:6034854/1‑76 (MQ=255)
                aCCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAgagct                                                             <  1:239461/76‑1 (MQ=255)
                aCCAGAATGGTGAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGCCTTCCCCCTGCCAGTAgagct                                                             <  2:4087529/76‑1 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAgagctg                                                            <  1:5878008/76‑1 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctg                                                            >  1:2030135/1‑76 (MQ=255)
                 ccAGAAGGGTTAATAGGCGTTAGCCCGGGCGGAAGGAATTCCCGCCTGGAATGCCTTCCCCCTGCCAGTAgagctg                                                            <  2:2860476/76‑1 (MQ=255)
                   aGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctgag                                                          >  1:596724/1‑76 (MQ=255)
                      aTGGTTAATAGGCGTTAGCCCGCGGGGAAGAAAATCTCCCCTGGAAGGACTTCCCCCTGCCAGTAGagctgagctg                                                       <  2:5132235/76‑1 (MQ=255)
                        ggTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCa                                                     >  1:4386565/1‑76 (MQ=255)
                        ggTAAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCa                                                     <  2:4642029/76‑1 (MQ=255)
                         gTTAATGGGCGTTAGACCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCaa                                                    <  1:757045/76‑1 (MQ=255)
                         gTTAATAGGCGTTAGCCCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCaa                                                    <  2:5114516/76‑1 (MQ=255)
                         gTTAATAGGCGTTAGCCCGCGCGAAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCaa                                                    <  1:3136596/76‑1 (MQ=255)
                          ttAATGGGCGTTAGACGGCGGGGAAGGAAATCCGGCCTGAAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAAc                                                   <  2:2159512/76‑1 (MQ=255)
                          ttAATAGGGGTAAGACCGCGCGGAAGGAAATCTCCCCTGGAATGCCTTCCCCCTGCCAGTAGAGCTGAGCTGCAAc                                                   <  2:4962875/76‑1 (MQ=255)
                          ttAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAc                                                   >  2:6684483/1‑76 (MQ=255)
                          ttAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAc                                                   >  2:7214418/1‑76 (MQ=255)
                           tAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcg                                                  >  2:3968371/1‑76 (MQ=255)
                           tAATAGGCGTTAGACCGCGCGGAAGGAAATCCCGCCGGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAAcg                                                  <  2:1752273/76‑1 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                 <  2:4386445/76‑1 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                 >  1:4862940/1‑76 (MQ=255)
                              tAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCgg                                               >  2:5678197/1‑76 (MQ=255)
                               aGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                              >  1:6710614/1‑76 (MQ=255)
                               aGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                              <  1:4149886/76‑1 (MQ=255)
                               aGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                              >  1:547374/1‑76 (MQ=255)
                                ggCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                             <  2:5958250/76‑1 (MQ=255)
                                ggCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                             >  1:7331565/1‑76 (MQ=255)
                                ggCGTTAGACCGCGCGGAAGGAAATCTCCCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                             >  2:1526250/1‑76 (MQ=255)
                                  cGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcag                                           >  1:2875067/1‑76 (MQ=255)
                                    ttAGCCCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                         <  2:3119575/76‑1 (MQ=255)
                                    ttAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                         >  1:4926476/1‑76 (MQ=255)
                                    ttAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                         >  1:1150068/1‑76 (MQ=255)
                                     tAGCCCGCGCGGAAGGAAATCCCGCCGGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcca                                        <  2:1847163/76‑1 (MQ=255)
                                     tAGCCCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcca                                        <  1:2761495/76‑1 (MQ=255)
                                     tAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcca                                        >  2:5289000/1‑76 (MQ=255)
                                     tAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcca                                        >  2:1446101/1‑76 (MQ=255)
                                      aGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagccag                                       >  1:188247/1‑76 (MQ=255)
                                      aGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagccag                                       >  1:1759313/1‑76 (MQ=255)
                                       gcccGCGCGGAAGAAAATCCCGCCGGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagccagc                                      <  1:5186440/74‑1 (MQ=255)
                                        cccGCGCGGAAGAAAATCTCGCCGGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCt                                     <  1:789113/75‑1 (MQ=255)
                                         ccGCGCGGAAGAAAATCTCGCCTGGAAGGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                                    <  1:1387936/76‑1 (MQ=255)
                                           gcgcgAAAGAAAATCTGCCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGc                                  <  1:7006232/76‑1 (MQ=255)
                                            cggggaaggaaTTCGCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                                 <  1:7465241/73‑1 (MQ=255)
                                               ggaagaaaTTCTCGCCGGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                              <  2:5990627/76‑1 (MQ=255)
                                               ggaagaaaACCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                              <  2:2983228/76‑1 (MQ=255)
                                                gaaggaaCTCCCGCCGGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAg                             <  2:2388826/76‑1 (MQ=255)
                                                gaaggaaATATCGCCTGGAATGTCTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAg                             <  1:5666101/76‑1 (MQ=255)
                                                 aaggaaATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGa                            >  1:6366406/1‑76 (MQ=255)
                                                  aggaaATCTCGCCTGGATTGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                           <  1:1183952/76‑1 (MQ=255)
                                                  aggaaATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                           <  1:3564492/76‑1 (MQ=255)
                                                  aggaaATCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                           <  2:7604484/76‑1 (MQ=255)
                                                  aagaaATCCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGa                            <  2:6930063/74‑1 (MQ=255)
                                                   gaaaaTCTCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATa                          <  1:1020213/74‑1 (MQ=255)
                                                    gaaATCCTCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATa                          <  2:2721740/76‑1 (MQ=255)
                                                    gaaATCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                         <  2:3399863/76‑1 (MQ=255)
                                                    aaaaTCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                         <  1:6124960/75‑1 (MQ=255)
                                                     aactctcCCCGGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGc                        <  2:6528444/73‑1 (MQ=255)
                                                     aaaTCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGc                        >  2:2110807/1‑76 (MQ=255)
                                                     aaaTCCCGCCGGAAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGc                        <  1:6582825/76‑1 (MQ=255)
                                                       aTCCCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCgg                      <  2:7173841/76‑1 (MQ=255)
                                                          ctcGCCGGAAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTa                    <  2:3249338/76‑1 (MQ=255)
                                                          cccGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTa                    <  1:7495600/74‑1 (MQ=255)
                                                            cGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTaaa                  >  1:7304141/1‑76 (MQ=255)
                                                             gCCTGAAATACCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAAt                 <  1:644575/76‑1 (MQ=255)
                                                                  gAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGct            >  1:290248/1‑76 (MQ=255)
                                                                  gAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGct            <  2:3674937/76‑1 (MQ=255)
                                                                   aaTGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGctc           <  2:1673343/76‑1 (MQ=255)
                                                                      gACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTcc        >  2:479521/1‑76 (MQ=255)
                                                                        cTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCgg      >  1:6249259/1‑76 (MQ=255)
                                                                            aCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGACgg  >  1:901537/1‑76 (MQ=255)
                                                                            |                                                                           
CGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGACGG  >  REL606/2144220‑2144370

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