Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,810,759 G→A intergenic (+80/‑11) ygaH → / → mprA putative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,810,7590GA94.2% 265.6 / 5.3 86intergenic (+80/‑11)ygaH/mprAputative L‑valine exporter, norvaline resistance protein/transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (35/46);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810663‑2810854
                                                                                                |                                                                                               
tGATGGCGATTATATAACCCACAAGAATCCTTTTTCTAAAACAATACATTTACTTTATTTGTCACTTTCGTTACTATATCGGCTGAAATTAATGAGATCa                                                                                              <  2:1382231/100‑1 (MQ=255)
tGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCa                                                                                              <  2:144826/100‑1 (MQ=255)
 gATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCAt                                                                                             <  1:1371967/100‑1 (MQ=255)
  aTGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACGGTCGTTACTATATTGGCTGAAATTAATGAGGTCATa                                                                                            >  1:1621595/1‑100 (MQ=255)
   tGGCGATTATAAAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAc                                                                                           >  2:1208470/1‑100 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTTTTTTTTTTTTCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCAt                                                                                             >  2:1705691/1‑97 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAcc                                                                                          <  2:866454/100‑1 (MQ=255)
    ggCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAcc                                                                                          <  2:1300679/100‑1 (MQ=255)
     cgatctatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAccc                                                                                         <  1:1375628/95‑1 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACCTTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCa                                                                                        <  2:765455/100‑1 (MQ=255)
      cGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATAcccn                                                                                        <  2:1448085/100‑2 (MQ=255)
      cGAGTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCa                                                                                        <  1:1086182/100‑1 (MQ=255)
        aTTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCaaa                                                                                      >  1:833840/1‑100 (MQ=255)
          tatataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATg                                                                                    <  2:2275100/100‑1 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATgg                                                                                   >  2:131086/1‑100 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATAGGCTGAAATTAATGAGATCATACCCagatcg                                                                                   >  2:1375625/1‑95 (MQ=255)
           atataACCCACAAGAATCATTTTTCTAAAACAAAACATTTACTTTATTTTTCACCGTCGTTACTATATCGGAGGAACTTAATGAGATCATACCCAAATgg                                                                                   >  1:1312786/1‑100 (MQ=255)
             ataACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATAAAGTCATACCCAAATGGAt                                                                                 <  1:1463951/100‑1 (MQ=255)
               aaCCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAg                                                                               <  1:1404436/100‑1 (MQ=255)
                  ccACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTc                                                                            >  1:469294/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCg                                                                           <  1:441508/100‑1 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCg                                                                           >  2:1795828/1‑100 (MQ=255)
                   cacaAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCg                                                                           >  2:1989129/1‑100 (MQ=255)
                     caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTGTATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGtt                                                                         <  1:846972/100‑1 (MQ=255)
                     caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGtt                                                                         <  1:342897/100‑1 (MQ=255)
                     caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGtt                                                                         <  2:833838/100‑1 (MQ=255)
                     caAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATACCGGCTGAAAGTAATGAGATCATACCCAAATGGATAGTTCGtt                                                                         >  2:1618739/1‑100 (MQ=255)
                      aaGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGttt                                                                        <  2:117364/100‑1 (MQ=255)
                       aGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTa                                                                       <  1:1705696/100‑1 (MQ=255)
                          aTCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGc                                                                    >  2:732513/1‑100 (MQ=255)
                           tCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGcc                                                                   >  2:23288/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGccc                                                                  >  1:1031340/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGccc                                                                  >  1:262602/1‑100 (MQ=255)
                            cATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGccc                                                                  <  2:353540/100‑1 (MQ=255)
                                tttCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCAttn                                                              <  2:1013748/100‑2 (MQ=255)
                                  tCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGaa                                                            >  2:2162641/1‑100 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           <  2:70132/100‑1 (MQ=255)
                                   cTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAAc                                                           <  2:1564263/100‑1 (MQ=255)
                                     aaaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaa                                                         >  1:2190335/1‑100 (MQ=255)
                                      aaaCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACaaa                                                        >  1:987096/1‑100 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  1:704206/100‑1 (MQ=255)
                                        aCAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATg                                                      <  1:2280377/100‑1 (MQ=255)
                                         cAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTGCTGAGATCATACCCGCATTGATAGTTCGTTTACGCCCATTGAACAAATGc                                                     <  2:1312783/100‑1 (MQ=255)
                                          aaTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCt                                                    <  1:131085/100‑1 (MQ=255)
                                           aTACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTa                                                   >  1:2201454/1‑100 (MQ=255)
                                             aCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaa                                                 <  2:1894835/100‑1 (MQ=255)
                                             aCATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaa                                                 >  1:1701886/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  2:1252709/1‑100 (MQ=255)
                                              cATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTaaaa                                                >  1:2083634/1‑100 (MQ=255)
                                                  tACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTc                                            >  1:889609/1‑100 (MQ=255)
                                                   aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                           <  2:276537/100‑1 (MQ=255)
                                                   aCTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcg                                           <  2:1287817/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAGTGGATAGTTCGTTTACGCCCATTGAACAAATGATAAAATTTcg                                           >  2:1951933/1‑99 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          <  2:262603/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          <  2:889607/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          <  2:1153808/100‑1 (MQ=255)
                                                    cTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgc                                          >  2:1246780/1‑100 (MQ=255)
                                                     tttATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTcgcg                                         <  1:1012823/100‑1 (MQ=255)
                                                       tATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGcc                                       >  1:916387/1‑100 (MQ=255)
                                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                     <  2:916385/100‑1 (MQ=255)
                                                         tttGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAg                                     <  2:1031338/100‑1 (MQ=255)
                                                            gTCACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccg                                  >  2:24591/1‑100 (MQ=255)
                                                              cACTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAgccgcc                                >  1:1568780/1‑100 (MQ=255)
                                                               aCTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCa                               <  1:1964771/100‑1 (MQ=255)
                                                               aCTGTCGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCa                               >  1:2276129/1‑100 (MQ=255)
                                                                    cGTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAg                          >  1:834840/1‑100 (MQ=255)
                                                                     gTTACTATATCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGa                         >  2:231271/1‑100 (MQ=255)
                                                                          tataTCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                    <  1:1460723/100‑1 (MQ=255)
                                                                          tataTCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                    >  2:101197/1‑100 (MQ=255)
                                                                          tataTCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTc                    >  1:588284/1‑100 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  2:2190325/100‑1 (MQ=255)
                                                                           atatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTcc                   <  2:1917774/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  1:732515/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  1:794038/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  1:24591/100‑1 (MQ=255)
                                                                            tatCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCt                  <  2:834838/100‑1 (MQ=255)
                                                                             atCGGCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCtt                 <  1:19768/100‑1 (MQ=255)
                                                                                ggCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTAt               >  1:478814/1‑99 (MQ=255)
                                                                                ggCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc              >  2:1722314/1‑100 (MQ=255)
                                                                                ggCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc              >  1:631901/1‑100 (MQ=255)
                                                                                ggCTGAAATTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATc              <  2:1774527/100‑1 (MQ=255)
                                                                                     aaaTTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgag         <  1:2039528/100‑1 (MQ=255)
                                                                                      aaTTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGgaga        <  2:588285/100‑1 (MQ=255)
                                                                                       aTTAATGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGAt       <  1:1252712/100‑1 (MQ=255)
                                                                                          aaTGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCt    <  2:1560221/100‑1 (MQ=255)
                                                                                            tGAGATCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCACCACGAAGATTTTCCTTATCAGGAGATCCTTc  >  2:1618040/1‑100 (MQ=255)
                                                                                                |                                                                                               
TGATGGCGATTATATAACCCACAAGAATCATTTTTCTAAAACAATACATTTACTTTATTTGTCACTGTCGTTACTATATCGGCTGAAATTAATGAGGTCATACCCAAATGGATAGTTCGTTTACGCCCATTGAACAAATGCTAAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTC  >  NC_000913/2810663‑2810854

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