Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606895,8560GC41.3% 52.6 / 38.2 46A26P (GCC→CCC) ECB_00836putative Regulatory protein
Reads supporting (aligned to +/- strand):  ref base G (21/6);  new base C (0/19);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 5.07e-08
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 121 total aligned reads displayed.

ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGCCAG  >  REL606/895781‑895930
                                                                            |                                                                          
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                            <  1:1163302/76‑1 (MQ=255)
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                            <  1:1714962/76‑1 (MQ=255)
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                            <  1:2360431/76‑1 (MQ=255)
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                            <  1:6311956/76‑1 (MQ=255)
aTGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAgcg                                                                            <  1:3493167/76‑1 (MQ=255)
  gACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgcc                                                                          <  1:1289371/76‑1 (MQ=255)
   aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCCCTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgcca                                                                         <  1:3549709/76‑1 (MQ=255)
   aCGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgcca                                                                         <  1:2436112/76‑1 (MQ=255)
    cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccag                                                                        <  2:4243853/76‑1 (MQ=255)
    cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccag                                                                        <  2:2991837/76‑1 (MQ=255)
    cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccag                                                                        <  2:5406554/76‑1 (MQ=255)
    cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccag                                                                        <  1:7041296/76‑1 (MQ=255)
    cgcgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTAAACAGCgccag                                                                        <  2:4377239/76‑1 (MQ=255)
      cgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccagcc                                                                      <  1:1862199/76‑1 (MQ=255)
      cgTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccagcc                                                                      <  2:7409569/76‑1 (MQ=255)
       gTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCgccagcca                                                                     <  2:4718477/76‑1 (MQ=255)
          cgcTGGCAGTAGTGGTGGCGCTGGCCCTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCacac                                                                  <  1:388817/76‑1 (MQ=255)
          cgcTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCacac                                                                  <  1:2562156/76‑1 (MQ=255)
          cgcTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCacac                                                                  <  1:4195697/76‑1 (MQ=255)
          cgcTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCacac                                                                  <  1:6265458/76‑1 (MQ=255)
             tGGCGGTGGGGGGGGCGCGGGCCCTCGTTGCGGGGGGCTGGCAGTCGGGGCGGTTTAACACCCCCAGCCACACCAt                                                               <  2:752882/76‑1 (MQ=255)
             gggCGGTGGTGGTGCCGCGGGCCCCCGTGGCGCGGGGCGGGCGGGCGTGGCGTCTAAACGCCCCCAGCCACACCAt                                                               <  2:2518064/75‑1 (MQ=255)
                 agtagtGGTGGCGCTGGGACTCGTTGCGCTGGGCTGGCAGGGGTGGGGGGCGTAACAACCGCCGCCCCACCAACGa                                                            >  2:5171395/1‑76 (MQ=255)
                    ggtggggGCGCTGCCCCTGGTGGGGCGGGGTGGGCGGTCGGGGCGCCTAAACACCCCCAGCCACACCATCGAAACg                                                        <  2:7116687/75‑1 (MQ=255)
                    ggtggcgGCGCCGGCACTCGTGTGTCGGGGCTGGCGGTCGTGCCGCCTTACCGGCCCCAGCCACACCATCGAAACg                                                        <  2:8588332/75‑1 (MQ=255)
                     gtggtgGCGTTGGCACGCGTTGCGCGGGGCTGGCAGTCGTGGCTGCTAAACAGCGCCAGCCACACCATCGAAACGc                                                       <  1:3488496/76‑1 (MQ=255)
                      gggtggCGCGGGCCCTGGTTGCGGGGGGGGGGCGGGGGGGGCGGCTAACCACCCCCAGCCACACCATCGAAACGCa                                                      <  2:5268768/75‑1 (MQ=255)
                          ggCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAGCGCCCGCCACCCCCTCGAAAAGCAacgc                                                  >  1:7935569/1‑76 (MQ=255)
                          ggCGCTGCCCCTCGTTGCGCTGGGCGGGCAGTGGTGGGGCTTTAACACCCCCAGCCACACCATCGAAACGCAgcgc                                                  <  2:8885819/76‑1 (MQ=255)
                           gcgcTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAAACGCCCCCGCCCCCCCCCCGAAAACCAgcgcc                                                 >  1:6089582/1‑75 (MQ=255)
                           gcgcGGGCCCTCGTTGCGTGGGGCTGCCGGTCGGGGGGGCTAAACCCCCCCAGCCACACCATCGAAACGCAgcgcg                                                 <  2:854896/76‑1 (MQ=255)
                              cTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAGCGCCCGCCCCCCCCTCCCAAAGCAACGGgcgg                                              >  1:7700836/1‑76 (MQ=255)
                               tGGCACTCGTTGCGCTTGGCTGGCAGTCGTGGCGGGTTTCCCTCGGCACGCACCCCCTCGAAACGCAGCGCgcggc                                             >  2:2848195/1‑76 (MQ=255)
                               tGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTTAAAGCGCCCGCCACACCCACCAAAAACCGCGCgcggg                                             >  2:6465200/1‑75 (MQ=255)
                               tGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAACGCCCCCCACACCCACCAAAAGCAGCGCgcggc                                             >  1:935790/1‑76 (MQ=255)
                               tGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAACGCCAACCCCACCACCCAAAAGCAGCGCgcggg                                             >  2:231663/1‑75 (MQ=255)
                                  cgcTCGTTGCGCGGGGCGGGCATGCGGGGCGAATAAAGAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCt                                          <  1:9009214/74‑1 (MQ=255)
                                  cACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGGTTGAAACCACCCAGCACCCCCCCGAAACGCAGCGCGCGgcgcg                                          >  2:1259023/1‑75 (MQ=255)
                                   ccTCGTGGCGCGGGGCGGGCAGTGGGGGCGGCTTAACAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTg                                         <  2:6745568/75‑1 (MQ=255)
                                     tCGTTGCGGGGGGCGGGCGGCGGGGGGGGCTTACCAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGaa                                       <  1:2229255/76‑1 (MQ=255)
                                         tGCGCTGGGCTGGCAGTCGTGGCGGGTTTAAAACGCCCGCCCCACCCCCCAAACGCAGCGCCCGGGGCCGgaaaaa                                   >  2:5345763/1‑74 (MQ=255)
                                         tGCGCTGGGCTGGCAGTCGTGGCGGGTTAACAGCCCCCACCACACCCTCGAAACGCAGGGCGCGGCCACGAAAAga                                   >  1:1497894/1‑75 (MQ=255)
                                         tGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCCCACACACCATCGAAACGCAGCGGGCGGCGCTGGCAAGt                                   >  1:2897044/1‑76 (MQ=255)
                                         gcggctgggtgggCGGTGGGGGCGGTTAAACGCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGt                                   <  2:8149108/73‑1 (MQ=255)
                                          gcgcGGGGCGGCCGGTGGGGGCGCCTTACCGCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTa                                  <  2:8276514/76‑1 (MQ=255)
                                           cgttgggcgggCGGTCGGGGCGGCTAAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTaa                                 <  2:712941/73‑1 (MQ=255)
                                           cgcTGGGCTGGCAGTCGTGGCGGGTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGaaa                                 >  1:3763481/1‑76 (MQ=255)
                                           ccgtgggctggCAGCGGTGCCGGCTTAACCCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTaa                                 <  1:4382446/73‑1 (MQ=255)
                                            ggggggctggCGGTGGGGGCGGCTTAACGCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTaaa                                <  2:2130140/73‑1 (MQ=255)
                                            ggggggcgggcGGCCGGGGCGTCTCAACGGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTaaa                                <  1:8176937/73‑1 (MQ=255)
                                            gctgggctggCAGTCGTGGCGGGTTAACAGCGCCAGCCACACCATCGAAACGCCGCGCGCCGCGCTGGAAAGTaaa                                >  1:6587413/1‑76 (MQ=255)
                                             gtgggcaggCGGGTGTGGCGTATAAACCGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAg                               <  2:3440598/75‑1 (MQ=255)
                                             cggggctggCAGTCGTGGCGGTTTAACAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAg                               <  2:8767351/74‑1 (MQ=255)
                                              ggggctggCAGTCGGGGCGTCTACACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgc                              <  1:4366404/75‑1 (MQ=255)
                                               ggggtggCCGTCGGCGCGACTAACCACCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAgcg                             <  1:7901450/76‑1 (MQ=255)
                                                  ggggCAGTCGTGGCGGCTAACCGCCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAg                          <  1:8150068/74‑1 (MQ=255)
                                                  ctggCAGTCGTGGCGGGTTTAAAGCGCCCGCCACACCCTCCAAACGCCGCGCGCGGGGCTGCAAAGTAAAgagcca                          >  2:2311998/1‑74 (MQ=255)
                                                   tggCAGTCGTGGCGGGTTTACAGCGCCCGCCACACCATCGAAAAGCAGCGCCCGGCGCCGGAAAGAAAAgagccag                         >  2:1834929/1‑73 (MQ=255)
                                                    ggCTGTGGTGGTGGTTAAGCAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCACAGGa                        <  2:1643124/76‑1 (MQ=255)
                                                     gCAGTCGTGGCGGGTTTAAAACCCCACCCACACCCACCAAACGCAGGGGCCGGGGGTGGAAAAGAAAACGCAGGaa                       >  2:5285728/1‑76 (MQ=255)
                                                      cAGTCGTGGCGGGTTTACACCGCCAGCCCCACCATCCAAACGCAGCGCGCGGCGCTGAAAAAGAAAACgaacgaac                      >  2:1172391/1‑76 (MQ=255)
                                                      cAGTCGTGGCGGGTTAACAGCGCCAGCCACACCATCCAAAAGCAACGCGCGGCGCTGGAAAGCAAAGCGCAGGGAc                      >  2:4538084/1‑76 (MQ=255)
                                                      cAGTCGTGGCGGCTTAACAGCGCCAGCCACACCCTCGAAACGCAGCGCGCGGGGCCGGAAAGGAAAACGCAGGGAc                      >  2:6214898/1‑76 (MQ=255)
                                                       aGTCGTGGCGGGTTAACAGCGCCCGCCACACCCTCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAACGCAGGGACt                     >  2:6167850/1‑76 (MQ=255)
                                                       aGTCGTGGCGGCTTAACAGCGCCCGCCCCACCAACCAAACGCAGCGCCCGGGGCTGGAAAGTAAAACGCAAGacat                     >  2:3550121/1‑73 (MQ=255)
                                                        ggggTGGCGGCACACCCGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTg                    <  1:1200479/73‑1 (MQ=255)
                                                        gTCGTGGCGGGTTTACAGCGCCCGCCACCCCATCCAAAAGCAGAGCGCGGCGCTGAAAAGAAAAGCGCAGGAACTg                    >  2:7799593/1‑76 (MQ=255)
                                                          cGTGGCGGGTTAACAGCGCCCGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCCGGAAATGAc                  >  2:1701281/1‑76 (MQ=255)
                                                          cGTGGCGGGTTAACAAAGCCCGCCACACCATCGAAACGCAGCGCGCGGCGCCGAAAAGGAAAGAGCAGGAACTGaa                  >  2:5623731/1‑75 (MQ=255)
                                                           gTGGCGGGTTTACAACGCCCGCCCCACCATCGAAACGCAGCGCGCGGCGCTGGAAAGTAAAGCGCAGGAACTGaag                 >  2:1374113/1‑74 (MQ=255)
                                                           gTGGCGGGTTTACAACGCCAGCCACACCATCCAAAAGCAGCGCGCGGCGCCGGAAAGAAAAAGGCAGGGACTGACg                 >  2:8954429/1‑76 (MQ=255)
                                                           gTGGCGGGTTAACAGCGCCCGCCCCACCACCCAAACGCAGCGCGCGGCGCCGGAAAAGAAAGAGCAGGAACTGACg                 >  2:1621174/1‑76 (MQ=255)
                                                            tGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACga                >  2:2728057/1‑76 (MQ=255)
                                                              ggcggcTTAACAGCGCCAGCCACACCATTGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa               >  1:3406390/1‑76 (MQ=255)
                                                              ggcggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTgaagaa               >  2:7755358/1‑76 (MQ=255)
                                                              ggcggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa               >  2:742480/1‑76 (MQ=255)
                                                              ggcggcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa               >  2:234849/1‑76 (MQ=255)
                                                              gccggcTTACCCGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaa               <  1:4005592/74‑1 (MQ=255)
                                                               ggagcTTACCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACgaag              <  2:1461517/73‑1 (MQ=255)
                                                                  gcTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaa           >  2:6148813/1‑76 (MQ=255)
                                                                   cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa          >  1:1505151/1‑76 (MQ=255)
                                                                   cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa          >  1:7370063/1‑76 (MQ=255)
                                                                   cTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaa          >  1:1701123/1‑76 (MQ=255)
                                                                    ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa         >  2:6285797/1‑76 (MQ=255)
                                                                    ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa         >  1:5671970/1‑76 (MQ=255)
                                                                    ttAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa         >  2:8151209/1‑76 (MQ=255)
                                                                    aaaaCAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGaaaaa         <  1:7599010/74‑1 (MQ=255)
                                                                     tAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAc        >  2:2662919/1‑76 (MQ=255)
                                                                      aaCCGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca       <  1:1951863/76‑1 (MQ=255)
                                                                      aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca       >  2:3483440/1‑76 (MQ=255)
                                                                      aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca       >  1:5171559/1‑76 (MQ=255)
                                                                      aaCAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAca       >  2:105906/1‑76 (MQ=255)
                                                                        cAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagc     >  2:8952227/1‑76 (MQ=255)
                                                                        cAGCCCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagc     <  1:6040239/76‑1 (MQ=255)
                                                                         aGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagcc    >  1:5012244/1‑76 (MQ=255)
                                                                         aGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagcc    >  1:6302495/1‑76 (MQ=255)
                                                                         aGCCCCACCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagcc    <  2:4757500/76‑1 (MQ=255)
                                                                           cgcCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagccag  <  1:4224690/76‑1 (MQ=255)
                                                                           cgcCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagccag  >  1:2730654/1‑76 (MQ=255)
                                                                           cgcCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAAcagccag  >  2:3406738/1‑76 (MQ=255)
                                                                            |                                                                          
ATGACGCGTGCGCTGGCAGTAGTGGTGGCGCTGGCACTCGTTGCGCTGGGCTGGCAGTCGTGGCGGCTTAACAGCGCCAGCCACACCATCGAAACGCAGCGCGCGGCGCTGAAAAGTAAAGCGCAGGAACTGACGAAGAAAAACAGCCAG  >  REL606/895781‑895930

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