Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 428,188 G→T R181L (CGC→CTC)  secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913428,1880GT92.7% 264.8 / 7.6 82R181L (CGC→CTC) secDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (42/34);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 9.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGC  >  NC_000913/428092‑428281
                                                                                                |                                                                                             
aaaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGa                                                                                            <  1:268443/100‑1 (MQ=255)
 aaaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGAt                                                                                           <  1:1918649/100‑1 (MQ=255)
  aaTATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATg                                                                                          <  1:978023/100‑1 (MQ=255)
    tatCGATAGCCTACGCAGTGACCTGCACGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCt                                                                                        >  2:1183643/1‑100 (MQ=255)
      tCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaa                                                                                      >  1:965960/1‑100 (MQ=255)
      tCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaa                                                                                      >  2:165215/1‑100 (MQ=255)
       cGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTaaa                                                                                     >  2:840455/1‑100 (MQ=255)
        gATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAg                                                                                    >  2:94038/1‑100 (MQ=255)
        gATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAg                                                                                    >  2:1328801/1‑100 (MQ=255)
         aTAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGc                                                                                   <  1:881408/100‑1 (MQ=255)
          tAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCt                                                                                  <  2:1147516/100‑1 (MQ=255)
            gCCTCCGCAGTGACCTGCGCGAACAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTc                                                                                 >  1:1558603/1‑99 (MQ=255)
            gCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCg                                                                                <  1:878107/100‑1 (MQ=255)
             ccTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGt                                                                               >  2:634745/1‑100 (MQ=255)
                aCGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGTTCGtgat                                                                            <  2:1026774/100‑1 (MQ=255)
                 cGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatg                                                                           <  1:366900/100‑1 (MQ=255)
                  gCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGtgatga                                                                          >  2:195022/1‑100 (MQ=255)
                       gACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCa                                                                     >  2:393882/1‑100 (MQ=255)
                       gACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCa                                                                     >  2:1076421/1‑100 (MQ=255)
                          cTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTg                                                                  >  1:1009294/1‑100 (MQ=255)
                           tGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGc                                                                 >  1:1571718/1‑100 (MQ=255)
                           tGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGc                                                                 >  2:708413/1‑100 (MQ=255)
                             cgcgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGt                                                               >  1:1540422/1‑100 (MQ=255)
                             cgcgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGt                                                               >  1:941996/1‑100 (MQ=255)
                               cgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTAt                                                             >  2:510218/1‑100 (MQ=255)
                               cgAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTAt                                                             >  1:1094027/1‑100 (MQ=255)
                                gAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATc                                                            <  2:264838/100‑1 (MQ=255)
                                  aaaGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGTGATGATAAAGCTCGTGATGAAGCCATTGCGTATCTg                                                          >  1:1307395/1‑100 (MQ=255)
                                    aGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAg                                                        >  1:1911941/1‑100 (MQ=255)
                                         aTCCCGTATACCACTGTTCGTAAAGAAAACAACTGCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgc                                                   <  1:950551/100‑1 (MQ=255)
                                         aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCCCTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTCCGTATCTGAGCCAgc                                                   >  1:183320/1‑100 (MQ=255)
                                         aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgc                                                   <  2:798693/100‑1 (MQ=255)
                                         aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgc                                                   >  1:417335/1‑100 (MQ=255)
                                         aTCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgc                                                   >  1:1844152/1‑100 (MQ=255)
                                          tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTGAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                  <  1:1038198/100‑1 (MQ=255)
                                          tCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                  >  1:1267269/1‑100 (MQ=255)
                                           cccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcgc                                                 >  1:486985/1‑100 (MQ=255)
                                           cccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAgcg                                                  >  1:949866/1‑99 (MQ=255)
                                            ccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGcc                                                <  1:1231164/100‑1 (MQ=255)
                                            ccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGcc                                                >  1:1824614/1‑100 (MQ=255)
                                            ccGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGcc                                                <  2:860983/100‑1 (MQ=255)
                                                  aCCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCgg                                          <  1:773742/100‑1 (MQ=255)
                                                   ccACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGa                                         >  1:342145/1‑100 (MQ=255)
                                                     aCTGTTCGTAAAGAAAACAACGACGGCCTAAGCATCACTTTCCGTGATGTGAAAGCTCGTGATGAATCCATTGCGTATCTGAGCAAGCGCCATCCgg                                          <  1:671914/97‑1 (MQ=255)
                                                      cTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGAc                                        <  2:11271/98‑1 (MQ=255)
                                                      cTGGTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCt                                      <  2:514599/100‑1 (MQ=255)
                                                        gTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTgg                                    <  1:1716533/100‑1 (MQ=255)
                                                         ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                   <  1:220922/100‑1 (MQ=255)
                                                         ttCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGt                                   <  1:1317788/100‑1 (MQ=255)
                                                          tCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTg                                  >  1:611621/1‑100 (MQ=255)
                                                            gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                <  1:502880/100‑1 (MQ=255)
                                                            gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                <  2:592737/100‑1 (MQ=255)
                                                            gTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAt                                >  2:150109/1‑100 (MQ=255)
                                                             tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                               <  2:250855/100‑1 (MQ=255)
                                                             tAAAGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGAtt                               <  1:386152/100‑1 (MQ=255)
                                                              aaaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTa                              <  2:1056518/100‑1 (MQ=255)
                                                               aaGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTag                             >  1:1066290/1‑100 (MQ=255)
                                                                aGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCTATTGCGTATCTGAGCAAGCGCCATCCGGTCCTGGTGATTacc                            <  2:1316655/100‑3 (MQ=255)
                                                                aGAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagc                            >  1:315084/1‑100 (MQ=255)
                                                                 gAAAACAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagca                           >  1:71809/1‑100 (MQ=255)
                                                                  acaacaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcag                          <  1:1183641/98‑1 (MQ=255)
                                                                   aaaCAACTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTagcagc                         >  2:245922/1‑100 (MQ=255)
                                                                    aacaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGcc                        <  1:266710/100‑1 (MQ=255)
                                                                    aacaacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGcc                        >  1:1262006/1‑100 (MQ=255)
                                                                      caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                      >  2:1030786/1‑100 (MQ=255)
                                                                      caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                      <  2:6989/100‑1 (MQ=255)
                                                                      caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                      <  2:267660/100‑1 (MQ=255)
                                                                      caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                      <  1:1773806/100‑1 (MQ=255)
                                                                      caacTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAg                      >  1:1800153/1‑100 (MQ=255)
                                                                       aacTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAgg                     <  2:185573/100‑1 (MQ=255)
                                                                        acTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAggg                    >  1:816646/1‑100 (MQ=255)
                                                                         cTACGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgc                   >  2:98053/1‑100 (MQ=255)
                                                                           aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                 <  1:694072/100‑1 (MQ=255)
                                                                           aCGGCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcag                 <  1:305148/100‑1 (MQ=255)
                                                                             ggCCTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGgcagca               <  2:846704/100‑1 (MQ=255)
                                                                                cTGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAAcc            >  2:1603834/1‑100 (MQ=255)
                                                                                 tGAGCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCa           <  2:981122/100‑1 (MQ=255)
                                                                                    gCATCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCt        >  2:853070/1‑100 (MQ=255)
                                                                                       tCACTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcg     <  2:1066291/100‑1 (MQ=255)
                                                                                         aCTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtg   >  1:897113/1‑100 (MQ=255)
                                                                                          cTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCtgcgtgc  <  1:165214/100‑1 (MQ=255)
                                                                                          cTTTCCTCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATCAGCAGCCAGGGCAGCAACCAGCtgcgtgc  >  1:1624641/1‑100 (MQ=255)
                                                                                                |                                                                                             
AAAATATCGATAGCCTACGCAGTGACCTGCGCGAAAAGGGCATCCCGTATACCACTGTTCGTAAAGAAAACAACTACGGCCTGAGCATCACTTTCCGCGATGCTAAAGCTCGTGATGAAGCCATTGCGTATCTGAGCAAGCGCCATCCGGACCTGGTGATTAGCAGCCAGGGCAGCAACCAGCTGCGTGC  >  NC_000913/428092‑428281

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: