Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,302,867 A→T D100E (GAT→GAA phnD ← phosphonate/organophosphate ester transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,302,8670AT99.2% 487.3 / NA 134D100E (GAT→GAAphnDphosphonate/organophosphate ester transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (90/43);  minor base G (1/0);  total (91/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

CGATCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGATCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCTGGATA  >  REL606/4302794‑4302939
                                                                         |                                                                        
cGATCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTc                                                                        >  1:1867631/1‑76 (MQ=255)
 gATCACGACGCCCCAGTAACCTGGCGACCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTcc                                                                       <  1:2545545/76‑1 (MQ=255)
 gATCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTcc                                                                       >  2:1605755/1‑76 (MQ=255)
  aTCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCa                                                                      >  2:4543644/1‑76 (MQ=255)
  aTCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCa                                                                      >  1:3899628/1‑76 (MQ=255)
  aTCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCa                                                                      >  2:2407588/1‑76 (MQ=255)
  aTCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGACGTTGGAGCGTTCCa                                                                      >  2:5556791/1‑76 (MQ=255)
   tCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCAc                                                                     >  2:2415240/1‑76 (MQ=255)
     aaCACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGGTCCACCg                                                                   >  2:1440113/1‑76 (MQ=255)
        gcGCTCCAGTAACCGGCCGACCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCtt                                                                <  2:2624865/75‑1 (MQ=255)
        ccGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCtt                                                                <  1:4708853/75‑1 (MQ=255)
        acGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCtt                                                                >  2:775797/1‑76 (MQ=255)
         cGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTc                                                               >  1:4420845/1‑76 (MQ=255)
         cGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTc                                                               >  1:6010368/1‑76 (MQ=255)
          gCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTcc                                                              >  1:4956027/1‑76 (MQ=255)
            tCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCAt                                                            >  1:3005802/1‑76 (MQ=255)
            cccAGTACCCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCAt                                                            <  2:605039/75‑1 (MQ=255)
             ccAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCAtt                                                           >  1:3481665/1‑76 (MQ=255)
             ccAGTAACCGGGCGACCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCAtt                                                           <  1:5703253/76‑1 (MQ=255)
              cAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTg                                                          >  2:3491058/1‑76 (MQ=255)
               aGTAGCCGGGCGATCCACCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGc                                                         <  1:4191235/76‑1 (MQ=255)
               aGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTACACCGCTTCCCTTGc                                                         >  1:2694035/1‑76 (MQ=255)
                gTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGcc                                                        >  2:4230963/1‑76 (MQ=255)
                gTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGcc                                                        >  1:4298945/1‑76 (MQ=255)
                 tAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCg                                                       >  1:4965384/1‑76 (MQ=255)
                 tAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCg                                                       >  2:274956/1‑76 (MQ=255)
                 tAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCg                                                       >  1:1078327/1‑76 (MQ=255)
                 tAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCAACGCTTCCATTGCCg                                                       >  1:4055840/1‑76 (MQ=255)
                  aaCCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGa                                                      >  2:2749016/1‑76 (MQ=255)
                  aaCCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGa                                                      >  2:3549315/1‑76 (MQ=255)
                    ccgggCGATCCATCCGCCGCGACCGTCTGGGCGAAGGCCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACa                                                    <  1:2688872/73‑1 (MQ=255)
                    ccTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGGCa                                                    >  2:3363938/1‑76 (MQ=255)
                    ccTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACa                                                    >  1:5593142/1‑76 (MQ=255)
                    ccTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGACGACa                                                    >  1:3794449/1‑76 (MQ=255)
                     cTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAg                                                   >  1:4344784/1‑76 (MQ=255)
                     cTGGCGACCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAg                                                   <  2:2659754/76‑1 (MQ=255)
                       ggCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGt                                                 >  2:235390/1‑76 (MQ=255)
                          gATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGc                                              >  2:1716845/1‑76 (MQ=255)
                           tgccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGcc                                             <  1:4903369/74‑1 (MQ=255)
                            cccacccGCGGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCg                                            <  2:1078327/75‑1 (MQ=255)
                             ccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCgg                                           >  1:2147186/1‑75 (MQ=255)
                             ccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGt                                           >  1:3009575/1‑76 (MQ=255)
                             ccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGt                                           >  2:3410820/1‑76 (MQ=255)
                             ccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGt                                           >  1:3702046/1‑76 (MQ=255)
                             ccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGt                                           >  2:3841738/1‑76 (MQ=255)
                             ccatccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGGTGCCGt                                           >  1:466497/1‑76 (MQ=255)
                               atccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTAc                                         >  1:2074634/1‑76 (MQ=255)
                               atccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTAc                                         >  1:6181431/1‑76 (MQ=255)
                               atccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTAc                                         >  1:5635243/1‑76 (MQ=255)
                               atccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTAc                                         >  1:715585/1‑76 (MQ=255)
                               atccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTAc                                         >  2:112766/1‑76 (MQ=255)
                                tccGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTAcc                                        <  1:5466463/76‑1 (MQ=255)
                                 ccgccgCGCCCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCa                                       <  1:5905585/76‑1 (MQ=255)
                                 ccgccgCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCa                                       >  1:1454500/1‑76 (MQ=255)
                                 ccgccgCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCa                                       <  2:4764947/76‑1 (MQ=255)
                                 ccgccgCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCa                                       <  1:4923473/76‑1 (MQ=255)
                                 ccgccgCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCa                                       <  1:6246374/76‑1 (MQ=255)
                                  cgccgcGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAg                                      >  1:3300960/1‑76 (MQ=255)
                                   gccgcGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAgg                                     <  2:656879/76‑1 (MQ=255)
                                    ccgcGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCCGGc                                    >  2:1475442/1‑76 (MQ=255)
                                    ccgcGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGc                                    >  1:5565847/1‑76 (MQ=255)
                                     cgcgCCCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCg                                   <  2:2313042/76‑1 (MQ=255)
                                        gcccGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGata                                <  2:6267463/74‑1 (MQ=255)
                                        gACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGTGatt                                >  1:4094301/1‑75 (MQ=255)
                                        gACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGata                                >  1:2355561/1‑76 (MQ=255)
                                           cGTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATcc                             >  1:5667388/1‑76 (MQ=255)
                                            gTCTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCa                            >  1:2877261/1‑76 (MQ=255)
                                            gTCGGGGCGAAGCCCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCa                            <  1:775797/76‑1 (MQ=255)
                                             tcggggCGAAGCCCTGGCCGTTGGCGCGTTCCCCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCAc                           <  2:5593142/73‑1 (MQ=255)
                                              cTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACt                          >  2:2105489/1‑76 (MQ=255)
                                              cTGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACt                          >  2:6278353/1‑76 (MQ=255)
                                               tGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACtt                         >  2:4798978/1‑76 (MQ=255)
                                               tGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACtt                         <  2:2196153/76‑1 (MQ=255)
                                               tGGGCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACtt                         >  1:2637697/1‑76 (MQ=255)
                                                gggCGAAGCCCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACttt                        <  1:5405311/76‑1 (MQ=255)
                                                 ggCAAAGACCTGGCGGTTGGCGCGTTCCCCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTa                       <  1:955150/76‑1 (MQ=255)
                                                  gCGAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTAt                      >  1:2931772/1‑76 (MQ=255)
                                                    gAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTg                    >  2:3813066/1‑76 (MQ=255)
                                                    gAAGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTg                    >  1:3398029/1‑76 (MQ=255)
                                                     aaGACCTGGCCGTTGGCGCGTTCCGCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGa                   >  2:3353103/1‑76 (MQ=255)
                                                     aaGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGa                   >  2:1159734/1‑76 (MQ=255)
                                                     aaGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGa                   >  2:5232621/1‑76 (MQ=255)
                                                     aaGACCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGa                   >  2:135073/1‑76 (MQ=255)
                                                      aggccTGGCCGTTGGCGCGTTCCACCGCTCCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGaa                  <  2:2944234/73‑1 (MQ=255)
                                                       gcccTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAg                 <  2:1852943/74‑1 (MQ=255)
                                                        aCCTGGTCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAgc                >  1:3479382/1‑76 (MQ=255)
                                                        aCCTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAgc                <  1:2304473/76‑1 (MQ=255)
                                                          cTGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAgcgc              >  1:1678477/1‑76 (MQ=255)
                                                           tGGCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCa             >  2:1664180/1‑76 (MQ=255)
                                                            ggCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCAt            >  1:848797/1‑76 (MQ=255)
                                                             gCCGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCAtt           <  2:658318/76‑1 (MQ=255)
                                                              ccGTTGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTc          <  2:3556637/76‑1 (MQ=255)
                                                               cGTTGGCGCTTTCCCCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTcc         <  1:5399280/76‑1 (MQ=255)
                                                                gtgggCGCGTTCCCCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTccc        <  2:3719610/73‑1 (MQ=255)
                                                                 ttGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTCCCGGGCGATATCCACTTTATTGAAGCGCATTCCCt       <  2:5461155/76‑1 (MQ=255)
                                                                 tgggCGCGTTCCCCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCt       <  2:1885839/74‑1 (MQ=255)
                                                                  tGGCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCTg      >  1:1514467/1‑76 (MQ=255)
                                                                  gggCGCGTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCTg      <  2:3389433/75‑1 (MQ=255)
                                                                      gcgTTCCGCCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCTGGATa  <  1:1099261/76‑1 (MQ=255)
                                                                      gcgTTCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCTGGATa  <  2:5974834/76‑1 (MQ=255)
                                                                         |                                                                        
CGATCAACACGCTCCAGTAACCTGGCGATCCATCCGCCGCGACCGTCTGGGCGAAGACCTGGCCGTTGGCGCGATCCACCGCTTCCATTGCCGACAGGTTGCCGTACCAGGCGATATCCACTTTATTGAAGCGCATTCCCTGGATA  >  REL606/4302794‑4302939

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