Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 942,738 C→T R594W (CGG→TGG)  dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913942,7380CT100.0% 215.6 / NA 56R594W (CGG→TGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/27);  total (29/27)

GGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGCGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGCAGGCAAACC  >  NC_000913/942592‑942882
                                                                                                                                                  |                                                                                                                                                
ggAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCa                                                                                                                                              <  1:1043923/151‑1 (MQ=255)
      tGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGt                                                                                                                                         >  2:1164594/1‑150 (MQ=255)
        tCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATg                                                                                                                                      >  2:821356/1‑151 (MQ=255)
                gATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTcg                                                                                                                              >  1:32001/1‑151 (MQ=255)
                   tttCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCgg                                                                                                                           >  1:1043005/1‑151 (MQ=255)
                    ttCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCggg                                                                                                                          >  2:697930/1‑151 (MQ=255)
                      cAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGa                                                                                                                         <  2:454139/150‑1 (MQ=255)
                       aaCGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGa                                                                                                                         >  1:455963/1‑149 (MQ=255)
                       aaCGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAg                                                                                                                       <  2:860425/151‑1 (MQ=255)
                        aCGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGc                                                                                                                      >  1:337560/1‑151 (MQ=255)
                           aTCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGAt                                                                                                                   <  2:1093302/151‑1 (MQ=255)
                                  gATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGaa                                                                                                            <  2:415488/151‑1 (MQ=255)
                                     tAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTg                                                                                                         >  1:850424/1‑151 (MQ=255)
                                         ccGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCa                                                                                                      >  1:64440/1‑150 (MQ=255)
                                          cgcgCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAAc                                                                                                    <  2:278925/151‑1 (MQ=255)
                                           gcgcTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACg                                                                                                   >  2:764770/1‑151 (MQ=255)
                                               tttGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGttt                                                                                                >  2:417561/1‑150 (MQ=255)
                                                  gAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGaa                                                                                             <  1:71605/150‑1 (MQ=255)
                                                   aaTGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAaga                                                                                           <  1:603862/151‑1 (MQ=255)
                                                   aaTGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAaga                                                                                           <  2:235394/151‑1 (MQ=255)
                                                      gTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGtt                                                                                        <  2:222113/151‑1 (MQ=255)
                                                        aaGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTcc                                                                                      >  1:202068/1‑151 (MQ=255)
                                                             cATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCa                                                                                   >  2:741442/1‑149 (MQ=255)
                                                             cATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAg                                                                                 <  1:821932/151‑1 (MQ=255)
                                                                    gAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGa                                                                          <  1:27545/151‑1 (MQ=255)
                                                                         gACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCttt                                                                     >  2:945881/1‑151 (MQ=255)
                                                                         gACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCttt                                                                     <  1:491790/151‑1 (MQ=255)
                                                                                  aCTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGGTCTTTAAAAAgcg                                                             >  1:1061972/1‑150 (MQ=255)
                                                                                  aCTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAgcgc                                                            >  1:400736/1‑151 (MQ=255)
                                                                                  aCTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAgcgc                                                            <  1:82740/151‑1 (MQ=255)
                                                                                    tGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGcgcg                                                           >  2:538804/1‑150 (MQ=255)
                                                                                        aaaaCGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCca                                                      <  2:444222/151‑1 (MQ=255)
                                                                                         aaaCGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCGCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCcac                                                     <  1:789485/151‑1 (MQ=255)
                                                                                          aaCGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCcaca                                                    <  2:845661/151‑1 (MQ=255)
                                                                                            cGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCcaca                                                    >  2:506087/1‑149 (MQ=255)
                                                                                            cGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAg                                                  <  1:741841/151‑1 (MQ=255)
                                                                                              tCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAggg                                                >  1:31374/1‑151 (MQ=255)
                                                                                                  ggCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATc                                            <  2:1159992/151‑1 (MQ=255)
                                                                                                  ggCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATc                                            <  2:619606/151‑1 (MQ=255)
                                                                                                       tGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGtt                                       >  2:1026465/1‑151 (MQ=255)
                                                                                                       tGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAACGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGt                                        <  1:673275/150‑1 (MQ=255)
                                                                                                            aaCAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCtt                                   >  2:788469/1‑150 (MQ=255)
                                                                                                            aaCAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTta                                  <  2:1146276/151‑1 (MQ=255)
                                                                                                             aCAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTtat                                 <  2:569388/151‑1 (MQ=255)
                                                                                                             aCAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTtat                                 >  2:129281/1‑151 (MQ=255)
                                                                                                               aGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATaa                               <  1:414722/151‑1 (MQ=255)
                                                                                                                 tttACTGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATaa                               >  2:79787/1‑149 (MQ=255)
                                                                                                                      tGAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc                        >  1:280923/1‑151 (MQ=255)
                                                                                                                       gAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc                        <  1:1061992/150‑1 (MQ=255)
                                                                                                                       gAAGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTc                        <  1:1049967/150‑1 (MQ=255)
                                                                                                                         aGGCCGTACCCAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGt                     <  2:198598/151‑1 (MQ=255)
                                                                                                                                   cAGGAAGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGc           >  2:482085/1‑151 (MQ=255)
                                                                                                                                       aaGAGTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGCAGGc       >  2:219667/1‑151 (MQ=255)
                                                                                                                                          agTGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGCAGGCaa     >  1:658201/1‑150 (MQ=255)
                                                                                                                                            tGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGCAGGCAAAcc  >  2:436080/1‑151 (MQ=255)
                                                                                                                                            tGGATGTGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGCAGGCAAAcc  >  1:858065/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGAATATGTCTTACGTGATTTTCAACGATCAGGTGATTAAACCGCGCTTTGAATGTAAGACCATCTATGAAATGACCAGCGAACTGGCAAAACGTCTTGGCGTTGAGCAACAGTTTACTGAAGGCCGTACCCAGGAAGAGTGGATGCGGCATCTGTATGCCCAGTCGCGGGAAGCGATTCCTGAACTGCCAACGTTTGAAGAGTTCCGCAAGCAGGGGATCTTTAAAAAGCGCGACCCACAAGGGCATCACGTTGCTTATAAAGCCTTCCGTGAAGATCCGCAGGCAAACC  >  NC_000913/942592‑942882

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