Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,144,2950AC43.2% 309.7 / 122.3 220G357G (GGT→GGGmrpantiporter inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (101/24);  new base C (0/95);  total (101/119)
Fisher's exact test for biased strand distribution p-value = 1.10e-39
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 259 total aligned reads displayed.

CGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGACG  >  REL606/2144220‑2144369
                                                                           |                                                                          
cGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCa                                                                            <  2:1386762/76‑1 (MQ=255)
cGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCa                                                                            <  2:5229736/76‑1 (MQ=255)
cGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCa                                                                            >  2:2344501/1‑76 (MQ=255)
  aTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAcc                                                                          >  2:6951221/1‑76 (MQ=255)
  aTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAcc                                                                          >  1:1626640/1‑76 (MQ=255)
    gTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCt                                                                        <  1:9153728/76‑1 (MQ=255)
    gTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCt                                                                        >  1:852934/1‑76 (MQ=255)
     tATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTg                                                                       <  1:1646125/76‑1 (MQ=255)
      aTGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGc                                                                      >  1:2127050/1‑76 (MQ=255)
      aTGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGc                                                                      >  1:367831/1‑76 (MQ=255)
       tGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGcc                                                                     <  1:8854038/76‑1 (MQ=255)
       tGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGcc                                                                     <  1:7687652/76‑1 (MQ=255)
        gaagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                                    <  1:7514890/76‑1 (MQ=255)
         aagaaTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                                   <  2:2978347/76‑1 (MQ=255)
             aTAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTaga                                                               <  1:5871736/76‑1 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAgagctg                                                           <  2:3679163/76‑1 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctg                                                           >  2:4966642/1‑76 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctg                                                           >  2:9514434/1‑76 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctg                                                           >  2:5836964/1‑76 (MQ=255)
                 ccAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCGGGAATGACTTCCCCCTGCCAGTAgagctg                                                           <  2:3426798/76‑1 (MQ=255)
                   aGAATGGTTAATAGGCGTTAGCCCGCGCGGAAGGAAATCTCGCCTGGAATGCCTTCCCCCTGCCAGTAgagctgag                                                         <  2:5680554/76‑1 (MQ=255)
                   aGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAACTCTCGCCTGGAATGCCTTCCCCCTGCCAGTAgagctgag                                                         <  2:1759636/76‑1 (MQ=255)
                   aGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctgag                                                         >  1:9339336/1‑76 (MQ=255)
                   aGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctgag                                                         >  1:6105193/1‑76 (MQ=255)
                   aGAATGGTTAATAGGCGTTAGACCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAgagctgag                                                         <  2:9485034/76‑1 (MQ=255)
                     aaTGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAgagctgagct                                                       <  2:1569138/76‑1 (MQ=255)
                      aTGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGagctgagctg                                                      <  1:4500799/76‑1 (MQ=255)
                      aTGGTTAATAGGCGTAAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGagctgagctg                                                      <  1:2491012/76‑1 (MQ=255)
                       tGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGc                                                     >  2:5337176/1‑76 (MQ=255)
                          ttAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCt                                                                        <  1:9382112/54‑1 (MQ=255)
                          ttAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAc                                                  >  2:321838/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:598043/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:9187124/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:8073446/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:6430437/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:2373554/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:3263897/1‑76 (MQ=255)
                            aaTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgc                                                >  2:4259812/1‑76 (MQ=255)
                             aTAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAcgcg                                               >  1:4102860/1‑76 (MQ=255)
                              tAGGCGTCAGCCCGCGCGGAAGGAAAGCTCCCCTGGAATGTCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCgg                                              <  2:1927907/76‑1 (MQ=255)
                               aGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                             >  2:2568820/1‑76 (MQ=255)
                               aGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                             >  1:8375448/1‑76 (MQ=255)
                               aGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                             >  1:7661391/1‑76 (MQ=255)
                                ggCGTTAGCCCGCGCGGAAGGAAATCTCGCCTGGAATGCCTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                            <  1:6370414/76‑1 (MQ=255)
                                ggCGTTAGCCCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                            <  1:5373972/76‑1 (MQ=255)
                                ggCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                            <  1:5552390/76‑1 (MQ=255)
                                ggCGTTAGACCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTc                                            <  1:1936034/76‑1 (MQ=255)
                                 gCGTTAGGCCGCGCGGAAGAAAACCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTca                                           <  2:2140083/76‑1 (MQ=255)
                                  cGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCCGTATAGCTGAGCTGCAACGCGGTcag                                          >  1:10508/1‑76 (MQ=255)
                                  cGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcag                                          >  1:5299582/1‑76 (MQ=255)
                                   gTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagc                                         >  1:4521122/1‑76 (MQ=255)
                                    ttAGGCCGCGGGGAAGAAAATCTTCCCGGTAATGCTTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                        <  1:181293/76‑1 (MQ=255)
                                    ttAGCCCGCGCGGAAGAAAATCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                        <  2:1554643/76‑1 (MQ=255)
                                    ttAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                        <  2:2695013/76‑1 (MQ=255)
                                    ttAGACCGCGCGGAAGGAAATCTCGCCGGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                        <  2:865521/76‑1 (MQ=255)
                                    ttAGACCGCGCGGAAGAAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                                        <  1:6924862/76‑1 (MQ=255)
                                     tAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcca                                       >  2:5102407/1‑76 (MQ=255)
                                     tAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcca                                       >  2:7207745/1‑76 (MQ=255)
                                      aGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagccag                                      >  1:3072433/1‑76 (MQ=255)
                                         ccGGGCGGAAGAAAATCTCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                                   <  2:1734504/76‑1 (MQ=255)
                                           gcgggGAGGGAAATCACGCCTGGAAGGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGc                                 <  2:8847233/76‑1 (MQ=255)
                                           gcgcgGAAGGAAATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGc                                 <  2:2056966/76‑1 (MQ=255)
                                            ggcgAAAGGAAGCCTCACCTGGAATGCTTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                                <  2:859815/75‑1 (MQ=255)
                                            cgcgGAGGAAAATCTCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                                <  2:6196065/76‑1 (MQ=255)
                                            cgcgGAAGGAAATCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                                <  2:5900694/76‑1 (MQ=255)
                                            cgcgGAAGAAAACCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                                <  2:747205/76‑1 (MQ=255)
                                            cgcgAAAGGAAATCTCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                                <  2:5296311/76‑1 (MQ=255)
                                             gcgGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCgg                               >  2:3478295/1‑76 (MQ=255)
                                              ggaaagaaaACCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGt                              <  2:1034359/73‑1 (MQ=255)
                                               ggaaggaaATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                             >  2:7963105/1‑76 (MQ=255)
                                                gaagaaaATCGCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAg                            <  1:4416445/76‑1 (MQ=255)
                                                 aaggaaATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGa                           >  1:1223990/1‑76 (MQ=255)
                                                 aaggaaATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGa                           >  1:7416069/1‑76 (MQ=255)
                                                 aagaaaACCTCCCCTGAAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGa                           <  1:3989204/76‑1 (MQ=255)
                                                  aggaaATCTCGCCGGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                          <  2:7817136/76‑1 (MQ=255)
                                                  agaaaaCCTCCCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                          <  2:6878026/73‑1 (MQ=255)
                                                   ggaaATCTCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATa                         <  1:4772002/76‑1 (MQ=255)
                                                    gaaTTCTCGCGTGGAATGCCTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:4002354/76‑1 (MQ=255)
                                                    gaaGTCGCGCCTAAAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:6056672/76‑1 (MQ=255)
                                                    gaaGTCCCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:4832368/76‑1 (MQ=255)
                                                    gaaATCTCGCCTGGAATGCCTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:92618/76‑1 (MQ=255)
                                                    gaaATCTCGCCTGGAATGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:8776225/76‑1 (MQ=255)
                                                    gaaACCGCGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:4725343/76‑1 (MQ=255)
                                                    aaaaTCTCGCCTGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:7452890/75‑1 (MQ=255)
                                                    aaaaTCTCGCCGGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg                        <  2:9403973/75‑1 (MQ=255)
                                                      aacctcGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCg                      <  2:7100235/73‑1 (MQ=255)
                                                      aaTCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCg                      >  1:6618875/1‑76 (MQ=255)
                                                      aaTCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCg                      >  1:3614328/1‑76 (MQ=255)
                                                           cgcgTGGAATCCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTaaa                 <  2:6105193/76‑1 (MQ=255)
                                                           cGCCTGGAATGACTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTaaa                 <  2:9339336/76‑1 (MQ=255)
                                                             ccTGGAATGCCTTCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAAtt               <  2:7913705/76‑1 (MQ=255)
                                                               tGGAATGCCTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCg             <  1:5067105/76‑1 (MQ=255)
                                                                ggAATGACTTCACCCTGCCAGTAGGGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGc            >  2:8521506/1‑76 (MQ=255)
                                                                ggAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGc            >  2:4641784/1‑76 (MQ=255)
                                                                ggAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGc            >  1:2868196/1‑76 (MQ=255)
                                                                   aTGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGctct         >  2:4397686/1‑76 (MQ=255)
                                                                   aTGACTCCCCCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGctct         <  2:4439703/76‑1 (MQ=255)
                                                                     gACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTcc       >  1:2720760/1‑76 (MQ=255)
                                                                         tCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGAc   >  2:7695930/1‑76 (MQ=255)
                                                                         tCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGAc   >  1:3289622/1‑76 (MQ=255)
                                                                           |                                                                          
CGATGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCCGGACG  >  REL606/2144220‑2144369

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