Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 547,835 +G coding (255/261 nt) b4507 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913547,8311.G96.0% 400.8 / 10.1 154coding (251/261 nt)b4507hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (5/1);  new base G (81/65);  total (86/68)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 171 total aligned reads displayed.

AAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCC  >  NC_000913/547756‑547900
                                                                            |                                                                      
aaGGGGCGTGCGGGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCa                                                                         <  2:444818/76‑1 (MQ=255)
  tgggCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGACGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                      <  2:2655414/75‑1 (MQ=255)
  ggggCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                      <  1:1178630/76‑1 (MQ=255)
  ggggCGGGCGTGTGCGCATGTCGGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                      <  1:1182386/76‑1 (MQ=255)
   gggCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                                    <  2:4768269/76‑1 (MQ=255)
    ggCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                                   >  1:135219/1‑76 (MQ=255)
      cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATCCTGGAGCAggggga                                                                 <  2:3183814/76‑3 (MQ=255)
       gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                                                >  2:3651781/1‑73 (MQ=255)
       gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                                                >  2:3885179/1‑73 (MQ=255)
       gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                                                <  2:2217906/76‑4 (MQ=255)
          cgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCCCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTg                                                              <  1:2254009/76‑1 (MQ=255)
           gtgGCCGCATGACTGTTCGAAGGTTGGGCGAACGATGAATCGCCGGCGCTGGGAATACTGGAGCAGGGGGAAGTGa                                                             <  2:5128140/76‑1 (MQ=255)
            ggggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                            <  1:648406/75‑1 (MQ=255)
             gggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAc                                                           >  1:368950/1‑76 (MQ=255)
              ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACCGGAGCAGGGGGAAGAGAACt                                                          >  2:5391111/1‑76 (MQ=255)
                cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGTAGTGAACGTc                                                        >  1:997094/1‑76 (MQ=255)
                cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                                        >  2:2404711/1‑76 (MQ=255)
                  cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGGACTTCAt                                                      >  1:3189834/1‑76 (MQ=255)
                  cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAt                                                      >  1:2838048/1‑76 (MQ=255)
                   aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAATGAACTTCAtt                                                     >  2:766367/1‑76 (MQ=255)
                   aTGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGTGGAAGTGAACTTCAtt                                                     >  1:3920176/1‑76 (MQ=255)
                     gTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTcc                                                   <  1:430743/76‑1 (MQ=255)
                      tCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCt                                                  >  2:3418250/1‑76 (MQ=255)
                       cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                                 >  2:1772339/1‑76 (MQ=255)
                       cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                                 <  2:748504/76‑1 (MQ=255)
                       cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                                 <  1:5014206/76‑1 (MQ=255)
                       cTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                                 <  1:4784893/76‑1 (MQ=255)
                        tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                                >  2:1655864/1‑76 (MQ=255)
                         gTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGt                                               <  1:3797336/76‑1 (MQ=255)
                          ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTc                                              >  1:1461195/1‑76 (MQ=255)
                           tCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                             <  2:4213023/76‑1 (MQ=255)
                           tCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCa                                             >  1:5021888/1‑76 (MQ=255)
                            cGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAc                                            >  1:5729510/1‑76 (MQ=255)
                             gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                           >  2:5657173/1‑76 (MQ=255)
                             gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCAcc                                           <  2:2698106/76‑1 (MQ=255)
                              aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAGTTCATTCCTTGTCACCa                                          >  1:4901470/1‑76 (MQ=255)
                              aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                          >  1:3345017/1‑76 (MQ=255)
                              aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCa                                          <  1:1585756/76‑1 (MQ=255)
                               aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAAGTTCATTCCTTGTCACCAt                                         <  2:4901469/76‑1 (MQ=255)
                               aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                         <  1:4780001/76‑1 (MQ=255)
                               aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                         >  1:28743/1‑76 (MQ=255)
                               aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                         <  1:1737652/76‑1 (MQ=255)
                               aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGCGAACTTCATTCCTTGTCAccct                                         >  2:4983210/1‑74 (MQ=255)
                                ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                        >  2:576222/1‑76 (MQ=255)
                                ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                        >  2:4532353/1‑76 (MQ=255)
                                 gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                       >  2:3962610/1‑76 (MQ=255)
                                 gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGCGAACTTCATTCCTTGTCACCCtgt                                       >  2:5620770/1‑76 (MQ=255)
                                   tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGGAGTGAACTTCATTCCTTGTCACCATGTGa                                     >  2:1075544/1‑76 (MQ=255)
                                    gggCGAAAGATGAAGCGCAGGCGCTGGCGATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGaa                                    <  1:3120320/76‑1 (MQ=255)
                                     ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                                   >  1:3431171/1‑76 (MQ=255)
                                     ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAt                                   <  1:199306/76‑1 (MQ=255)
                                      gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATg                                  <  2:2672925/76‑1 (MQ=255)
                                      gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGAGGACGTGAACTTCATTCCTTGTCACCATGTGaagg                                  >  1:4427431/1‑74 (MQ=255)
                                        gagaGATGAAGCGCAGGCGCTGGCAATACTGGGGCATGGGGAAGTGACCTTCATTCCTTGTCACCATGTGAATGcc                                <  1:5695337/73‑1 (MQ=255)
                                         aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                               <  2:2989888/76‑1 (MQ=255)
                                          aataTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                              >  1:4034365/4‑76 (MQ=255)
                                          aaGATGAGGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                              <  2:1894097/76‑1 (MQ=255)
                                          aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                                <  1:3591328/58‑1 (MQ=255)
                                          aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                                >  2:3591328/1‑58 (MQ=255)
                                          aaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGt                              <  1:3049937/76‑1 (MQ=255)
                                           aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                             >  2:1785221/1‑76 (MQ=255)
                                             aTGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg                           <  2:5619450/76‑1 (MQ=255)
                                               gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                         >  1:2733262/1‑76 (MQ=255)
                                               gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                         <  2:1966322/76‑1 (MQ=255)
                                                aaGCGCAGGCGCTGGAAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCGCCATGTGAATGCCGTCGGGcc                        <  1:4549728/76‑1 (MQ=255)
                                                  gcgcAGGCGCTGGCCATACTGGAGCAGGGGGAAGCGCACTTCATTCCTTTCCCCCATGTGAAACCCCTCGCGCCac                      >  2:55253/1‑75 (MQ=255)
                                                  gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                      >  1:126688/1‑76 (MQ=255)
                                                   cgcAGGCGCTGGCAATACACGAGCAGGGGGCCGTGAACTTCATTCCTTGTCACCATGTGAATCCCGTCGGGcccag                     >  2:4270378/1‑73 (MQ=255)
                                                    gcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg                    <  2:2520288/76‑1 (MQ=255)
                                                    gcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg                    <  1:4459437/76‑1 (MQ=255)
                                                    gcAGGCGCGGGCAATACGGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATg                    <  2:8546/76‑1 (MQ=255)
                                                     cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATgg                   >  2:2270242/1‑76 (MQ=255)
                                                      aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg                  <  2:3550075/76‑1 (MQ=255)
                                                      aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg                  >  2:3791155/1‑76 (MQ=255)
                                                       ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc                 >  1:239860/1‑76 (MQ=255)
                                                       ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc                 >  2:4559093/1‑76 (MQ=255)
                                                        gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCCATGGGCg                >  2:1123482/1‑76 (MQ=255)
                                                        gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg                <  1:1772339/76‑1 (MQ=255)
                                                         cgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg               <  2:2306833/76‑1 (MQ=255)
                                                          gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt              >  2:4342564/1‑76 (MQ=255)
                                                          gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt              >  1:5463414/1‑76 (MQ=255)
                                                           cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCggga             >  2:1232683/1‑74 (MQ=255)
                                                            tGGCAATACTGGAGCAGGGGGAAGTGAGCTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt            <  2:4377013/76‑1 (MQ=255)
                                                            tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGCAt            >  2:1989518/1‑76 (MQ=255)
                                                            tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCAGTCGGGCCAATGGGCGGTAt            >  1:2585510/1‑76 (MQ=255)
                                                            tGGCAATACGGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt            <  1:1330311/76‑1 (MQ=255)
                                                                aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt        <  1:5254438/76‑1 (MQ=255)
                                                                 aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt       <  2:1340691/76‑1 (MQ=255)
                                                                 aTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt       >  2:2479617/1‑76 (MQ=255)
                                                                  tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc      <  2:1578674/76‑1 (MQ=255)
                                                                  tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc      >  1:3294497/1‑76 (MQ=255)
                                                                   aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCt     <  2:4704374/76‑1 (MQ=255)
                                                                   aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCt     >  1:1445893/1‑76 (MQ=255)
                                                                    cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTg    >  2:4430959/1‑76 (MQ=255)
                                                                    cTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTg    >  2:2543879/1‑76 (MQ=255)
                                                                    cTGGAGCAGGGGGAAGTGAACTTCATGCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTg    >  1:2897806/1‑76 (MQ=255)
                                                                     tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGc   >  2:2148724/1‑76 (MQ=255)
                                                                     tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGc   <  2:4043543/76‑1 (MQ=255)
                                                                     tGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAAAGCCGTCGGGCCAATGGGCGGTATTACTTCTGc   >  1:867689/1‑76 (MQ=255)
                                                                      ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGcc  >  1:1966187/1‑76 (MQ=255)
                                                                            |                                                                      
AAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCC  >  NC_000913/547756‑547900

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