Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606925,6390TG36.2% 216.9 / 33.4 116K121N (AAA→AACltaEL‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (22/52);  new base G (42/0);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 1.47e-15
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 165 total aligned reads displayed.

ACGCGGCAGCATTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAATACTG  >  REL606/925564‑925714
                                                                           |                                                                           
  gcgGCAGCATTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAAAGGATATCGTCGGGgtt                                                                           >  2:849672/1‑76 (MQ=255)
   cgGCAGCATTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTg                                                                          <  2:2071169/76‑1 (MQ=255)
   cgGCAGCATTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTg                                                                          <  1:2966061/76‑1 (MQ=255)
    ggCAGCATTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGa                                                                         <  2:6525268/76‑1 (MQ=255)
      cagcaTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGAtt                                                                       <  2:1847855/76‑1 (MQ=255)
       agcaTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGAttt                                                                      <  1:6137411/76‑1 (MQ=255)
       agcaTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGGttt                                                                      >  1:417422/1‑76 (MQ=255)
        gcaTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATTGATATCGTCGGGTTTGAtttt                                                                     <  2:7794748/76‑1 (MQ=255)
        gcaTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGAtttt                                                                     >  1:7594724/1‑76 (MQ=255)
         caTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGGttttt                                                                    >  1:7438816/1‑75 (MQ=255)
          aTTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGAttttta                                                                   >  2:4613418/1‑74 (MQ=255)
            tttGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATc                                                                 <  2:2743131/76‑1 (MQ=255)
            tttGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATc                                                                 <  2:5657065/76‑1 (MQ=255)
             ttGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGGTTTTTATCg                                                                >  2:3500799/1‑76 (MQ=255)
               gCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGGTTTTCATCTcc                                                              >  2:448445/1‑76 (MQ=255)
               gCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTCTCGcc                                                              >  2:3254142/1‑76 (MQ=255)
                ccGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCa                                                             <  1:6965521/76‑1 (MQ=255)
                ccGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGccc                                                             >  2:4357986/1‑75 (MQ=255)
                ccGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGccc                                                             >  2:7092937/1‑75 (MQ=255)
                 cGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATTGCCAc                                                            >  2:2979498/1‑76 (MQ=255)
                  gTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGGTTTTTTTCTCCACt                                                           >  1:4254163/1‑76 (MQ=255)
                  gTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCACt                                                           >  1:3276580/1‑76 (MQ=255)
                   ttGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCCCtt                                                          >  2:7694718/1‑76 (MQ=255)
                   ttGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTTATTTTCCTCGCCACtt                                                          >  2:681756/1‑76 (MQ=255)
                   ttGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATTGGCACtt                                                          >  2:6038937/1‑76 (MQ=255)
                   ttGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATTGCCACtt                                                          >  2:355315/1‑76 (MQ=255)
                    tgtgGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACttt                                                         <  1:6642718/76‑1 (MQ=255)
                    tgtgGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCACttt                                                         >  2:3290365/1‑76 (MQ=255)
                    tgtgGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGCCCCttt                                                         >  2:1222871/1‑76 (MQ=255)
                     gtgGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTa                                                        <  2:6975891/76‑1 (MQ=255)
                      tgGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGTTTTTTATTGCCACTTTAt                                                       >  1:454582/1‑76 (MQ=255)
                      tgGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGGCCCTCTAt                                                       >  2:6306856/1‑76 (MQ=255)
                       gggTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATc                                                      <  2:5054395/76‑1 (MQ=255)
                       gggTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGCCACTTTATc                                                      >  1:2265191/1‑76 (MQ=255)
                          tgtTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGTCACTTTATCCAg                                                   <  1:1237281/76‑1 (MQ=255)
                           gtTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGc                                                  <  2:3932536/76‑1 (MQ=255)
                            ttttCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCg                                                 <  2:7215165/76‑1 (MQ=255)
                             tttCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCACTTTATCCAGCgg                                                >  1:5343162/1‑76 (MQ=255)
                              ttCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGt                                               <  1:5946462/76‑1 (MQ=255)
                              ttCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGt                                               <  2:1720491/76‑1 (MQ=255)
                              ttCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGt                                               <  1:3664355/76‑1 (MQ=255)
                               tCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCggga                                              >  1:3813665/1‑74 (MQ=255)
                               tCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGCCCCTTTATCCAGCggga                                              >  1:7196800/1‑74 (MQ=255)
                                 cAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCACTTTATCCAGCGGTaac                                            >  2:1570224/1‑74 (MQ=255)
                                  aGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTTGTTTTCATCGCCACTTTATCCAGCGGTAGCg                                           >  1:953527/1‑76 (MQ=255)
                                  aGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATTGCCACTTTATCCAGCCGGAGCg                                           >  2:7035100/1‑76 (MQ=255)
                                   gACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGt                                          <  2:7511338/76‑1 (MQ=255)
                                     cTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGGTTTTTATCGCCACTTTATCCAGCGGGAGGGTGc                                        >  1:5340179/1‑76 (MQ=255)
                                     cTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGCCCCTTTATCCAGCGGGAGCGTGc                                        >  1:1116843/1‑76 (MQ=255)
                                      tGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGcc                                       >  2:1299716/1‑76 (MQ=255)
                                       gAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCCCTTTTTCCAACGGGAGGGTGGCg                                      >  1:4433111/1‑76 (MQ=255)
                                        aGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTTATCGGCACTTTATCCAGCCGGCGCGTTCCgc                                     >  2:7520887/1‑75 (MQ=255)
                                          tAATTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGCTTTTTATCGCCACTTTATCCAGCGGGAGCGTGCCGTCg                                   >  2:6731879/1‑76 (MQ=255)
                                           aaTTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCgg                                  <  1:3455493/76‑1 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGGTTTTCATCGCCACTTTATCCAGCGGTAGGGTGCCCGCGGc                                 >  2:6550280/1‑76 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGc                                 <  2:5197097/76‑1 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGGTTTATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCggg                                 >  2:4092882/1‑75 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCCCTTTATCCAGCGGTAGCGTGCCGTCggg                                 >  1:688236/1‑75 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCACTTTATTCAGCGGTAGCGTGCCGTCggg                                 >  2:6147469/1‑75 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCACTTTATCCAGCGGTAGCGTGCCCTCggg                                 >  2:5392498/1‑75 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCCTCGCCAATTTCTTCAACGGTAGCGCGCCGTCGGc                                 >  2:7353533/1‑76 (MQ=255)
                                            aTTTGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTCATCGCCACTTTATCCAGCGGGAGCGTGCCGTCggg                                 >  1:5294440/1‑75 (MQ=255)
                                             tttGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCa                                <  2:7859309/76‑1 (MQ=255)
                                              ttGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACCTTATCCAGCGGGAGCGTGCCGTCGGcgg                               >  1:5784548/1‑74 (MQ=255)
                                               tGGTGCGGGCGAAATGGATATCGTCGGGGTTGATTTTTATCGCCCCTTTATCCAGCGGTAACGTGTCGTCGGGAGc                              >  2:4390234/1‑76 (MQ=255)
                                                   gCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCcgcg                          <  1:5811575/76‑1 (MQ=255)
                                                   gCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCcgcg                          <  1:1644310/76‑1 (MQ=255)
                                                     gggCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTc                        <  1:647959/76‑1 (MQ=255)
                                                     gggCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTc                        <  1:4666488/76‑1 (MQ=255)
                                                      ggCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCg                       >  1:5621542/1‑76 (MQ=255)
                                                       gagAAATGGATATCGTCGGTTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGa                      <  2:4423456/74‑1 (MQ=255)
                                                        ggAAAGGAATATCGTCGGGTTTGATTTCCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGAt                     <  2:743105/75‑1 (MQ=255)
                                                        ctaaaTGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGAt                     <  1:3821869/74‑1 (MQ=255)
                                                         gAAATGGATATCGTCGGTTTTAATTTTCATCGCCCCTTTACCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATg                    <  2:407861/76‑1 (MQ=255)
                                                         gAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATg                    <  1:1906997/76‑1 (MQ=255)
                                                         gAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATg                    <  1:4157063/76‑1 (MQ=255)
                                                          aaaTGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATgg                   <  1:2040590/76‑1 (MQ=255)
                                                           aaTGGATATCGTCGGGTTTGATTTTTATCGCCACTTTATCCAGCGGTAGCGGGCCGTCGGGAGGCGCCGCGGTggg                  >  2:7306841/1‑76 (MQ=255)
                                                            aTGGATATCGTCGGTTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGt                 <  1:4420585/76‑1 (MQ=255)
                                                             tGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGtt                <  1:849672/76‑1 (MQ=255)
                                                              gtatatCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTg               <  1:3409903/74‑1 (MQ=255)
                                                               gttatCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGc              <  1:5043956/74‑1 (MQ=255)
                                                               gATATCGTCGGTTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGc              <  1:7190349/76‑1 (MQ=255)
                                                                   tcgtcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGGCGGTGGGGTGCgggt          >  1:5019314/1‑74 (MQ=255)
                                                                    cgtcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGGAGGCCCGGCGGtggggtgggggcg         >  2:1092210/1‑73 (MQ=255)
                                                                    cgtcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGCCGATGGGTTGCGGcgg         >  2:1903077/1‑74 (MQ=255)
                                                                    cgtcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCCCGTCGCTGGGGTGCGGCtt         >  2:4669566/1‑75 (MQ=255)
                                                                      tcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGTCGCCGCGTCGCTGGGTGGCGGCTgga       >  1:281287/1‑74 (MQ=255)
                                                                      tcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGaa       <  1:100225/76‑1 (MQ=255)
                                                                       ggtgTTTGATTTTCATCGCCCCTTTACCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAAt      <  2:2002635/73‑1 (MQ=255)
                                                                       ggggTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAAt      <  2:5731450/75‑1 (MQ=255)
                                                                       cgtgtttgtttTCCATCGCCGCTATTTCCAGTGTGACCGTGCAGTCGCCAGCCGCGTCGATGGGTTGCGGCTGAAt      <  2:2247693/73‑1 (MQ=255)
                                                                       cgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGGGGCTGAAt      >  1:7337972/1‑76 (MQ=255)
                                                                       cgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAAt      <  1:3930450/76‑1 (MQ=255)
                                                                       cgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGAGTTGCGGCTGaaa      >  1:5830008/1‑75 (MQ=255)
                                                                        gggTTTGTTTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAATa     <  2:6329422/76‑1 (MQ=255)
                                                                          gtttgtttTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAATACt   <  2:2403541/76‑1 (MQ=255)
                                                                          gTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGAAGGGGTGCGGCTGGAAACt   >  1:155445/1‑76 (MQ=255)
                                                                           tttGTTTTCCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAATACTg  <  2:6507663/76‑1 (MQ=255)
                                                                           tttGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGGATACTg  >  1:458875/1‑76 (MQ=255)
                                                                           |                                                                           
ACGCGGCAGCATTTTGCCGTTGTGGGTGTTTTCCAGACTGAGTAATTTGGTGCGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCGATGGGTTGCGGCTGAATACTG  >  REL606/925564‑925714

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