Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,309 A→G I84T (ATC→ACC)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,3090AG99.5% 902.4 / ‑5.1 216I84T (ATC→ACC) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (110/104);  total (110/106)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 221 total aligned reads displayed.

GATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGG  >  REL606/4202236‑4202381
                                                                         |                                                                        
gATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGtt                                                                        <  2:799788/76‑1 (MQ=255)
gATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGtt                                                                        <  2:2437539/76‑1 (MQ=255)
gATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGtt                                                                        <  2:6132492/76‑1 (MQ=255)
 aTAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTg                                                                       >  1:1207430/1‑76 (MQ=255)
 aTAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTg                                                                       <  2:735892/76‑1 (MQ=255)
 aTAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTg                                                                       <  1:4188209/76‑1 (MQ=255)
  tAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGc                                                                      >  1:4446810/1‑76 (MQ=255)
  tAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGc                                                                      >  1:3548877/1‑76 (MQ=255)
   aGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGcc                                                                     >  1:6943029/1‑76 (MQ=255)
   aGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGcc                                                                     <  1:5319279/76‑1 (MQ=255)
   aGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGcc                                                                     >  1:1855083/1‑76 (MQ=255)
   aGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGcc                                                                     >  2:2901820/1‑76 (MQ=255)
     ggTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCa                                                                   <  1:4290310/76‑1 (MQ=255)
     ggTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCa                                                                   <  2:517173/76‑1 (MQ=255)
     ggTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGGTGCCCa                                                                   >  1:5920084/1‑76 (MQ=255)
      gTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCaa                                                                  >  2:6679472/1‑76 (MQ=255)
       tGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAAt                                                                 <  2:3975706/76‑1 (MQ=255)
       tGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAAt                                                                 <  2:127230/76‑1 (MQ=255)
       tGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAAt                                                                 >  1:1649129/1‑76 (MQ=255)
       tGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAAt                                                                 <  2:3488729/76‑1 (MQ=255)
         aaCAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGt                                                               <  2:4868964/76‑1 (MQ=255)
           cAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTcc                                                             >  2:1062402/1‑76 (MQ=255)
           cAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTcc                                                             >  1:3438313/1‑76 (MQ=255)
             aTCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCa                                                           >  1:6878785/1‑76 (MQ=255)
             aTCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCa                                                           <  2:2967067/76‑1 (MQ=255)
              tCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAg                                                          >  1:2746312/1‑76 (MQ=255)
              tCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAg                                                          <  1:132544/76‑1 (MQ=255)
               cGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAgc                                                         <  1:500794/76‑2 (MQ=255)
               cGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGt                                                         <  2:1977217/76‑1 (MQ=255)
                gCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGtt                                                        <  1:3236587/76‑1 (MQ=255)
                gCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGtt                                                        <  2:1126992/76‑1 (MQ=255)
                 cTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTc                                                       >  1:3793943/1‑76 (MQ=255)
                 cTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTc                                                       >  2:6736692/1‑76 (MQ=255)
                  tAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCg                                                      <  2:1793295/76‑1 (MQ=255)
                  tAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCg                                                      >  1:105461/1‑76 (MQ=255)
                  tAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCg                                                      <  2:3196449/76‑1 (MQ=255)
                    aCAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGcc                                                    <  1:1124112/76‑1 (MQ=255)
                     cAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCg                                                   <  2:91484/76‑1 (MQ=255)
                     cAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCg                                                   >  1:3891666/1‑76 (MQ=255)
                      aaaTTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGa                                                  <  2:2971460/76‑1 (MQ=255)
                        aTTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGAcc                                                >  1:2554528/1‑76 (MQ=255)
                        aTTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGAcc                                                <  2:2765784/76‑1 (MQ=255)
                        aTTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGAcc                                                >  1:4395262/1‑76 (MQ=255)
                        aTTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGAcc                                                >  1:4650229/1‑76 (MQ=255)
                         ttACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCt                                               <  1:5793430/76‑1 (MQ=255)
                         ttACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCt                                               <  2:5102180/76‑1 (MQ=255)
                          tACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTg                                              <  2:5859659/76‑1 (MQ=255)
                          tACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTg                                              <  2:4683997/76‑1 (MQ=255)
                          tACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTg                                              <  2:4376716/76‑1 (MQ=255)
                          tACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTg                                              <  2:4314665/76‑1 (MQ=255)
                          tACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTg                                              >  2:413737/1‑76 (MQ=255)
                          tACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTg                                              >  2:1546381/1‑76 (MQ=255)
                             ggCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACg                                           >  2:3424352/1‑76 (MQ=255)
                             ggCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACg                                           >  1:2192451/1‑76 (MQ=255)
                             gcctctGGAGAAATCTGCTGCCGACCATAAAGGCGTGTGCGCCGTTTGCCCAATGTCCCAGTTCGCCGACCTGACg                                           <  1:1643384/74‑1 (MQ=255)
                               ctctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCa                                         <  2:6584349/76‑1 (MQ=255)
                                tctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCAc                                        >  1:1476853/1‑76 (MQ=255)
                                tctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCAc                                        >  1:6391919/1‑76 (MQ=255)
                                tctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCAc                                        >  2:3808618/1‑76 (MQ=255)
                                 ctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACg                                       <  1:1010609/76‑1 (MQ=255)
                                 ctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACg                                       <  1:4103226/76‑1 (MQ=255)
                                 ctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACg                                       <  2:6005586/76‑1 (MQ=255)
                                 ctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACg                                       >  2:4790506/1‑76 (MQ=255)
                                 ctGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACg                                       <  1:5820699/76‑1 (MQ=255)
                                   ggAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGaa                                     <  1:887265/76‑1 (MQ=255)
                                     agaAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAAc                                   >  2:6420876/1‑76 (MQ=255)
                                       aaaGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAAccc                                 >  1:4743560/1‑76 (MQ=255)
                                       aaaGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCAGACCTGACGCCCTAAAccc                                 >  2:530756/1‑76 (MQ=255)
                                        aaGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCct                                <  2:703319/76‑1 (MQ=255)
                                        aaGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACAAAACCct                                <  1:1906556/76‑1 (MQ=255)
                                        aaGCTCCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCct                                <  2:4393154/76‑1 (MQ=255)
                                         aGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCctg                               <  2:2452366/76‑1 (MQ=255)
                                         aGCTGCTGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCctg                               <  1:778010/76‑1 (MQ=255)
                                          gctgctGCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCctgc                              >  1:2332655/1‑76 (MQ=255)
                                           ctgctgCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCctgct                             >  2:5117696/1‑76 (MQ=255)
                                           ctgctgCCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCctgct                             >  1:6623059/1‑76 (MQ=255)
                                             gctgcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGt                           >  2:6171635/1‑76 (MQ=255)
                                             gctgcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGt                           <  2:5060360/76‑1 (MQ=255)
                                              ctgcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGtt                          >  2:5918862/1‑76 (MQ=255)
                                              ctgcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGtt                          <  1:3991736/76‑1 (MQ=255)
                                               tgcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTg                         <  2:2152118/76‑1 (MQ=255)
                                               ggcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTg                         <  1:3489241/75‑1 (MQ=255)
                                                gcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGc                        >  2:2552050/1‑76 (MQ=255)
                                                gcCGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGc                        >  2:3277446/1‑76 (MQ=255)
                                                 ccGACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCa                       <  2:6289095/76‑1 (MQ=255)
                                                   gACCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCAtt                     >  1:5381597/1‑76 (MQ=255)
                                                    aCCATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATtg                    <  2:2885491/76‑1 (MQ=255)
                                                     ccATAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATtgt                   >  1:5128357/1‑76 (MQ=255)
                                                       aTAAAGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTgg                 <  1:5637258/76‑1 (MQ=255)
                                                         aaaGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGtt               >  1:6754011/1‑76 (MQ=255)
                                                          aaGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTa              >  2:1621191/1‑76 (MQ=255)
                                                           aagCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAg             <  2:2282324/74‑1 (MQ=255)
                                                           aGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAg             <  1:2358018/76‑1 (MQ=255)
                                                           aGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAg             <  1:3966954/76‑1 (MQ=255)
                                                           aGGCATGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAg             <  1:4883663/76‑1 (MQ=255)
                                                               aTGTGCGCGGGTTGCCCAATGTCCCAGTAAGCTGGCCTGGCGCGCGAAACCCTGCTGTTGCATTGTGGTTAGCAgg         <  1:338267/76‑1 (MQ=255)
                                                               aTGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAgg         >  2:4595193/1‑76 (MQ=255)
                                                               aTGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAgg         >  1:6254519/1‑76 (MQ=255)
                                                               aTGTGCGCCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAgg         >  2:2133569/1‑76 (MQ=255)
                                                                     gcCGGTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTgg   >  2:1143924/1‑76 (MQ=255)
                                                                         |                                                                        
GATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGG  >  REL606/4202236‑4202381

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