Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,550,166 A→G T377T (ACT→ACC gntP ← fructuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,550,1660AG100.0% 209.6 / NA 58T377T (ACT→ACCgntPfructuronate transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/30);  total (28/30)

CAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAGCCATGCCATGATGTAGGGCGAGATATTGCCGTGGGACATCAACATG  >  NC_000913/4550023‑4550310
                                                                                                                                               |                                                                                                                                                
cAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTg                                                                                                                                           <  1:614740/151‑1 (MQ=255)
 aGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTgg                                                                                                                                          <  1:816271/151‑1 (MQ=255)
 aGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTgg                                                                                                                                          >  1:297146/1‑151 (MQ=255)
      ttCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCt                                                                                                                                      <  2:201521/150‑1 (MQ=255)
      ttCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCt                                                                                                                                      >  2:900763/1‑150 (MQ=255)
        cAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAg                                                                                                                                    <  2:297120/150‑1 (MQ=255)
        cAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAg                                                                                                                                    <  2:231874/150‑1 (MQ=255)
          aCGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGAt                                                                                                                                 <  2:491473/151‑1 (MQ=255)
          aCGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGAt                                                                                                                                 <  1:46413/151‑1 (MQ=255)
              gtCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTgc                                                                                                                              >  2:632007/1‑150 (MQ=255)
                cTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTgcagc                                                                                                                           >  1:407521/1‑151 (MQ=255)
                    tACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACt                                                                                                                        <  1:54805/150‑1 (MQ=255)
                     aCTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTg                                                                                                                       >  2:201632/1‑150 (MQ=255)
                     aCTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGa                                                                                                                      <  1:353529/151‑1 (MQ=255)
                      cTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGAt                                                                                                                     >  2:297761/1‑151 (MQ=255)
                       tGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGCTGCTGGATCTAGGATTGCAGCACTGAt                                                                                                                     >  2:847193/1‑150 (MQ=255)
                          cAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATc                                                                                                                 <  1:114803/151‑1 (MQ=255)
                               cAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCcggc                                                                                                            <  2:330907/151‑1 (MQ=255)
                                aaaGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCcggcg                                                                                                           >  1:203378/1‑151 (MQ=255)
                                   gTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCcggcggcg                                                                                                        <  2:2768/151‑1 (MQ=255)
                                     aaCCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCcggcggcgg                                                                                                       <  1:45480/150‑1 (MQ=255)
                                          ttGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCAt                                                                                                   <  1:6388/149‑1 (MQ=255)
                                            gAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGcc                                                                                               >  2:513063/1‑151 (MQ=255)
                                            gAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGcc                                                                                               >  2:150322/1‑151 (MQ=255)
                                             aacaacCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGcc                                                                                               <  2:691949/150‑1 (MQ=255)
                                                aacCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGaaa                                                                                           <  1:543411/151‑1 (MQ=255)
                                                    aGAAAGATGCATCGTTAATGTGGGTGAGGGCGTTGGAACCCGAAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAAGCCCGGCGGCGGTCATCGCCGAAACGa                                                                                        >  2:40339/1‑150 (MQ=255)
                                                          aTGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTg                                                                                 >  1:782761/1‑151 (MQ=255)
                                                          aTGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTg                                                                                 <  2:866030/151‑1 (MQ=255)
                                                            gCATCGGTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGAc                                                                               <  2:450116/151‑1 (MQ=255)
                                                                 gTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAGCGACACCCTGACCCGTc                                                                          >  2:551627/1‑151 (MQ=255)
                                                                   tAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCg                                                                         <  2:445514/150‑1 (MQ=255)
                                                                   tAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGc                                                                        >  1:730572/1‑151 (MQ=255)
                                                                             gAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAAtt                                                              <  2:525344/151‑1 (MQ=255)
                                                                               gggTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAg                                                            <  1:608585/151‑1 (MQ=255)
                                                                                ggTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGc                                                           >  1:629806/1‑151 (MQ=255)
                                                                                ggTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGc                                                           >  2:308163/1‑151 (MQ=255)
                                                                                   gtTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGcac                                                         >  2:431774/1‑150 (MQ=255)
                                                                                   gtTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGCCGCCAGACGAATTAGcac                                                         >  2:651901/1‑150 (MQ=255)
                                                                                    ttGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGcac                                                         >  1:739777/1‑149 (MQ=255)
                                                                                    ttGGAACCCGCAGCCGTAGCCAGTAGCAGCAGAGCCGGATGCACGCCAACCAGCCGACCGGTTGCTGGATCTAGGATAGCAACACTAACACTCCCGGCGGCGGTCAGTGCCGAATCGCTACCCTGACCCGTCGCCAGACGAATTAGca                                                          <  2:329275/148‑1 (MQ=255)
                                                                                         aCCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATc                                                  <  1:887082/151‑1 (MQ=255)
                                                                                          cccGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCa                                                 <  1:622641/151‑1 (MQ=255)
                                                                                           ccGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAg                                                <  2:186726/151‑1 (MQ=255)
                                                                                            cGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAgc                                               >  2:358230/1‑151 (MQ=255)
                                                                                                 ccGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAgccatgc                                          <  2:687024/151‑1 (MQ=255)
                                                                                                   gTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAgccatgcca                                        <  1:441186/151‑1 (MQ=255)
                                                                                                    tCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAgccatgccat                                       >  1:783504/1‑151 (MQ=255)
                                                                                                    tCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAgccatgccat                                       >  2:66660/1‑151 (MQ=255)
                                                                                                      gCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAg                                                            >  2:568610/1‑128 (MQ=255)
                                                                                                      gCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAg                                                            <  1:568660/128‑1 (MQ=255)
                                                                                                      gCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAgccatgccat                                       >  1:293179/1‑149 (MQ=255)
                                                                                                            aCCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGa                                                    >  1:449800/1‑130 (MQ=255)
                                                                                                            aCCAGCAGCGCCGGATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGa                                                    <  2:449753/130‑1 (MQ=255)
                                                                                                                        ggATTCACGCCAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAGCCATGCCATGATGTAGGGCGAGATAtt                     >  1:775489/1‑149 (MQ=255)
                                                                                                                                 ccAACCAGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAGCCATGCCATGATGTAGGGCGAGATATTGCCGTGGGACa          <  1:834411/151‑1 (MQ=255)
                                                                                                                                       aGCTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAGCCATGCCATGATGTAGGGCGAGATATTGCCGTGGGACATCAAc     >  2:21761/1‑150 (MQ=255)
                                                                                                                                         cTGACCGGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAGCCATGCCATGATGTAGGGCGAGATATTGCCGTGGGACATCAACATg  >  1:37191/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CAGGTTTTCAACGTGTCTTTTACTGACAGGTCAAAGTAACCTTTGAACAACCAGAAAGATGCATCGTTAATGTGGGTGAGGGTGTTGGAACCCGCAGCCGTCGCCAGTACCAGCAGCGCCGGATTCACGCCAACCAGCTGACCAGTTGCTGGATCTAGGATTGCAGCACTGATAATCCCGGCGGCGGTCATCGCCGAAACGACACCCTGACCCGTCGCCAGACGAATTAGCACAGTGATCAGCCATGCCATGATGTAGGGCGAGATATTGCCGTGGGACATCAACATG  >  NC_000913/4550023‑4550310

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