Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,227 T→C R737R (CGT→CGC aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,2270TC100.0% 241.5 / NA 67R737R (CGT→CGCaceEpyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/38);  total (29/38)

AACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCT  >  NC_000913/125096‑125349
                                                                                                                                   |                                                                                                                          
aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                     >  2:545008/1‑139 (MQ=255)
aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                     >  2:183140/1‑139 (MQ=255)
      aaCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               <  1:1159155/139‑1 (MQ=255)
       aCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                              <  1:545008/139‑1 (MQ=255)
       aCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                              <  1:713958/139‑1 (MQ=255)
       aCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                              <  1:117417/139‑1 (MQ=255)
                 gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               <  2:208825/128‑1 (MQ=255)
                 gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               >  1:208825/1‑128 (MQ=255)
                 gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                    >  1:1211716/1‑139 (MQ=255)
                 gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                    <  1:183140/139‑1 (MQ=255)
                        aTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGa                                                                                             <  1:806654/139‑1 (MQ=255)
                                     cTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGAtt                                                                                       <  1:353535/132‑1 (MQ=255)
                                     cTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGAtt                                                                                       >  2:353535/1‑132 (MQ=255)
                                        aGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                     <  1:1259657/131‑1 (MQ=255)
                                        aGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                     >  2:1259657/1‑131 (MQ=255)
                                          gAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                           <  2:1211716/139‑1 (MQ=255)
                                             ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               >  2:1125975/1‑100 (MQ=255)
                                             ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               <  1:1125975/100‑1 (MQ=255)
                                             ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTg                                                                        >  2:132524/1‑139 (MQ=255)
                                              gTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGa                                                                       <  2:1403095/139‑1 (MQ=255)
                                                  ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                    <  1:167185/90‑1 (MQ=255)
                                                  ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                    >  2:167185/1‑90 (MQ=255)
                                                  ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGtt                                                                   <  2:1122966/139‑1 (MQ=255)
                                                     tAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtat                                                                <  2:914196/139‑1 (MQ=255)
                                                     tAAAGGTATCTACAAACTCGAAACTATTGAAGGTAACAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                         >  1:658489/1‑130 (MQ=255)
                                                     tAAAGGTATCTACAAACTCGAAACTATTGAAGGTAACAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                         <  2:658489/130‑1 (MQ=255)
                                                       aaGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAg                                                              >  2:1063950/1‑139 (MQ=255)
                                                        aGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                             >  1:344714/1‑139 (MQ=255)
                                                        aGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                             <  1:1113330/139‑1 (MQ=255)
                                                           tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                          <  2:249957/139‑1 (MQ=255)
                                                           tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                          <  1:132524/139‑1 (MQ=255)
                                                           tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                          <  2:344714/139‑1 (MQ=255)
                                                           tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                          >  2:1196260/1‑139 (MQ=255)
                                                                   aaCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCtt                                                  >  2:290382/1‑139 (MQ=255)
                                                                    aCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTc                                                 >  1:652641/1‑139 (MQ=255)
                                                                       cGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                           >  1:1214253/1‑94 (MQ=255)
                                                                       cGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGaaa                                                                                           <  2:1214253/94‑1 (MQ=255)
                                                                         aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCa                                                >  2:306020/1‑135 (MQ=255)
                                                                         aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCa                                                <  1:306020/135‑1 (MQ=255)
                                                                         aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGa                                            >  1:1233302/1‑139 (MQ=255)
                                                                         aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGa                                            >  1:888359/1‑139 (MQ=255)
                                                                               ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               <  2:379241/66‑1 (MQ=255)
                                                                               ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               >  1:379241/1‑66 (MQ=255)
                                                                                    ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                 <  2:1233302/139‑1 (MQ=255)
                                                                                    ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                 <  1:621606/139‑1 (MQ=255)
                                                                                    ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                 <  2:888359/139‑1 (MQ=255)
                                                                                       aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               >  1:1440775/1‑58 (MQ=255)
                                                                                       aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                               <  2:1440775/58‑1 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGc                                                                                  >  1:1204886/1‑85 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGc                                                                                  <  2:1204886/85‑1 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                          >  1:73560/1‑109 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                          <  2:73560/109‑1 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCAcc                                              <  2:688722/121‑1 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                            <  2:1258983/139‑1 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                            >  1:1263988/1‑139 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                            <  1:290382/139‑1 (MQ=255)
                                                                                         cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                            >  1:371894/1‑139 (MQ=255)
                                                                                         cAAAGGTAAAATTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCAcc                                              >  1:688722/1‑121 (MQ=255)
                                                                                                 aaGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                              <  2:1363806/49‑1 (MQ=255)
                                                                                                 aaGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                              >  1:1363806/1‑49 (MQ=255)
                                                                                                       aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCTAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACg              <  1:469452/139‑1 (MQ=255)
                                                                                                       aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACg              <  2:1311995/139‑1 (MQ=255)
                                                                                                                   ccGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                    <  2:244093/41‑1 (MQ=255)
                                                                                                                   ccGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                    >  1:244093/1‑41 (MQ=255)
                                                                                                                   ccGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGAAACAtgct  <  1:1416560/139‑1 (MQ=255)
                                                                                                                      gTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                         <  2:1046072/65‑1 (MQ=255)
                                                                                                                      gTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                         >  1:1046072/1‑65 (MQ=255)
                                                                                                                                   |                                                                                                                          
AACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCT  >  NC_000913/125096‑125349

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: