Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,072,060 G→A L34L (CTG→TTG)  rutC ← putative aminoacrylate deaminase, reactive intermediate detoxification, weak enamine/imine deaminase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,072,0600GA98.6% 233.4 / NA 72L34L (CTG→TTG) rutCputative aminoacrylate deaminase, reactive intermediate detoxification, weak enamine/imine deaminase activity
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/35);  minor base T (0/1);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAGCGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGAC  >  NC_000913/1071929‑1072192
                                                                                                                                   |                                                                                                                                    
gCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCTTCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               >  2:969123/1‑139 (MQ=255)
gCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               <  2:176147/139‑1 (MQ=255)
gCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCCTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               >  2:1399286/1‑139 (MQ=255)
  tGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCTTCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGga                                                                                                                             <  1:969123/139‑1 (MQ=255)
    ttGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGgaga                                                                                                                           <  2:491194/139‑1 (MQ=255)
           tCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATacac                                                                                                                    >  1:1327985/1‑139 (MQ=255)
                    ccATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATc                                                                                                           >  1:567236/1‑139 (MQ=255)
                           gCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATcagccagc                                                                                                    >  2:1540086/1‑139 (MQ=255)
                               cccGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCTGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCATCGTACCGGAGACAt                                                                                                                        <  2:350613/115‑1 (MQ=255)
                               cccGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               <  1:440016/108‑1 (MQ=255)
                               cccGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               >  2:440016/1‑108 (MQ=255)
                               cccGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACAt                                                                                                                        >  1:350613/1‑115 (MQ=255)
                                 cgccgTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCg                                                                                              >  1:1473156/1‑139 (MQ=255)
                                               ccTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                <  1:710596/139‑1 (MQ=255)
                                               ccTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                <  2:567236/139‑1 (MQ=255)
                                                cTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGcc                                                                               >  1:1006439/1‑139 (MQ=255)
                                                  tGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAg                                                                                     <  2:1007273/131‑1 (MQ=255)
                                                  tGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAg                                                                                     >  1:1007273/1‑131 (MQ=255)
                                                    cGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCg                                                                           >  1:150372/1‑139 (MQ=255)
                                                      gATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGt                                                                         <  1:1399286/139‑1 (MQ=255)
                                                            ttCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCgcta                                                                   <  2:806338/139‑2 (MQ=255)
                                                                    cGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           <  2:904050/139‑1 (MQ=255)
                                                                    cGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           >  1:183025/1‑139 (MQ=255)
                                                                    cGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           >  1:1562336/1‑139 (MQ=255)
                                                                     gTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                <  1:1488474/117‑1 (MQ=255)
                                                                     gTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                >  2:1488474/1‑117 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGa                                                                                       >  1:121170/1‑108 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGa                                                                                       <  2:121170/108‑1 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                >  1:1482264/1‑115 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                <  2:1482264/115‑1 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAg                                                                             >  1:1197139/1‑118 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAg                                                                             <  2:1197139/118‑1 (MQ=255)
                                                                       ggCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCggg                                                        >  2:214674/1‑139 (MQ=255)
                                                                               tGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               <  2:692397/60‑1 (MQ=255)
                                                                               tGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCg                                                                                                                               >  1:692397/1‑60 (MQ=255)
                                                                               tGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTa                                                <  2:1473156/139‑1 (MQ=255)
                                                                                 cgcCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGTGAATAATTCCg                                              <  2:1006439/139‑1 (MQ=255)
                                                                                      ttGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTg                                                                                                 >  2:553074/1‑83 (MQ=255)
                                                                                      ttGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTg                                                                                                 <  1:553074/83‑1 (MQ=255)
                                                                                         ggTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTgg                                      <  2:1327985/139‑1 (MQ=255)
                                                                                              tCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCAt                                                                                                            <  2:507285/64‑1 (MQ=255)
                                                                                              tCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCAt                                                                                                            >  1:507285/1‑64 (MQ=255)
                                                                                              tCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                <  2:308744/92‑1 (MQ=255)
                                                                                              tCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                >  1:308744/1‑92 (MQ=255)
                                                                                              tCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCg                                 >  2:45301/1‑139 (MQ=255)
                                                                                                ggCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCg                                                                           >  1:503653/1‑95 (MQ=255)
                                                                                                ggCAAACAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCg                                                                           <  2:503653/95‑1 (MQ=255)
                                                                                                   aaaCAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCggg                                                        >  2:294344/1‑111 (MQ=255)
                                                                                                   aaaCAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCggg                                                        <  1:294344/111‑1 (MQ=255)
                                                                                                   aaaCAGCACGTTATTATGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTaaa                            <  2:150372/139‑1 (MQ=255)
                                                                                                                ttatGTTGATCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGCGat               <  2:1562336/139‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAg                                                                                     <  2:1653607/61‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAg                                                                                     >  1:1653607/1‑61 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGt                                                                         >  2:881251/1‑73 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGt                                                                         <  1:881251/73‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           >  1:229122/1‑87 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           <  2:907814/87‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           >  1:907814/1‑87 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGc                                                           <  2:229122/87‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACg                                              <  2:536824/100‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACg                                              >  1:536824/1‑100 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGc                  >  2:387029/1‑128 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGc                  <  1:387029/128‑1 (MQ=255)
                                                                                                                        aTCAAAAGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGCGATATGAGCaa       >  2:439840/1‑139 (MQ=255)
                                                                                                                          cAAAAGCCAACGTACCGGAGACATACACCACGCCATcagccagc                                                                                                    >  2:689117/1‑44 (MQ=255)
                                                                                                                          cAAAAGCCAACGTACCGGAGACATACACCACGCCATcagccagc                                                                                                    <  1:689117/44‑1 (MQ=255)
                                                                                                                             aaGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                >  1:795484/1‑61 (MQ=255)
                                                                                                                             aaGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGc                                                                                <  2:795484/61‑1 (MQ=255)
                                                                                                                             aaGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTaa                    <  2:1393529/121‑1 (MQ=255)
                                                                                                                             aaGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTaa                    >  1:1393529/1‑121 (MQ=255)
                                                                                                                             aaGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGAc  >  1:846338/1‑139 (MQ=255)
                                                                                                                             aaGCCAACGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGAc  <  1:1540086/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
GCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCATCGGCAAACAGCACGTTATTATGTTGATCAAAAGCCAGCGTACCGGAGACATACACCACGCCATCAGCCAGCGTGCCGGGAACGAAGGGGGCCAGCGGTGCGCTGCTGCCAGCGGGAATAATTACGGATTTTGGCATCGTTAAACTCCTTAAGCGATATGAGCAAAGGAC  >  NC_000913/1071929‑1072192

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