Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,544,384 C→T pseudogene (381/381 nt) yddK ← pseudogene, leucine‑rich protein, putative glycoportein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,544,3840CT100.0% 222.8 / NA 70pseudogene (381/381 nt)yddKpseudogene, leucine‑rich protein, putative glycoportein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/35);  total (35/35)

GCTAACCTTCAAACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTTTTATA  >  NC_000913/1544250‑1544519
                                                                                                                                      |                                                                                                                                       
gCTAACCTTCAAACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTaaa                                                                                                                                     >  1:801102/1‑139 (MQ=255)
gCTAACCTTCAAACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTaaa                                                                                                                                     >  2:1000380/1‑139 (MQ=255)
           aaCGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGaa                                                                                                                          <  1:424544/139‑1 (MQ=255)
                ggTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTgga                                                                                                                     >  2:340671/1‑139 (MQ=255)
                    aaTCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGaa                                                                                                                 <  1:309732/139‑1 (MQ=255)
                        tttGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACg                                                                                                             >  1:207990/1‑139 (MQ=255)
                          tGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGTAGGAAAACGAt                                                                                                           >  1:397908/1‑139 (MQ=255)
                          tGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGAt                                                                                                           >  2:137866/1‑139 (MQ=255)
                                   gTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                  <  1:137866/139‑1 (MQ=255)
                                       tGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGt                                                                                              <  1:704870/139‑1 (MQ=255)
                                           tCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTAt                                                                                          >  1:113115/1‑139 (MQ=255)
                                           tCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTAt                                                                                          <  1:688710/139‑1 (MQ=255)
                                                    cAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                  <  2:257687/138‑1 (MQ=255)
                                                    cAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                  >  1:257687/1‑138 (MQ=255)
                                                         cGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTa                                                                            <  2:992917/139‑1 (MQ=255)
                                                          gTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTAt                                                                           >  1:122060/1‑139 (MQ=255)
                                                            gAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGt                                                                                      >  2:993320/1‑126 (MQ=255)
                                                            gAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATATGTCCAGAAGGTCTATCGGt                                                                                      <  1:993320/126‑1 (MQ=255)
                                                            gAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATAAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTg                                                                         >  1:656274/1‑139 (MQ=255)
                                                             aaTCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                                        <  1:753714/139‑1 (MQ=255)
                                                                         tACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                            <  2:626938/139‑1 (MQ=255)
                                                                                  aCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTTTCTTATTGCCTGATTGTAATAATCGttatt                                                   <  2:397908/139‑1 (MQ=255)
                                                                                  aCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCt                                                                      >  1:905185/1‑120 (MQ=255)
                                                                                  aCTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCt                                                                      <  2:905185/120‑1 (MQ=255)
                                                                                    tataTTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTAc                                                 <  2:207990/139‑1 (MQ=255)
                                                                                     atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                                    <  1:297944/135‑1 (MQ=255)
                                                                                     atatTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGttat                                                    >  2:297944/1‑135 (MQ=255)
                                                                                      tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                  <  1:46477/88‑1 (MQ=255)
                                                                                      tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGa                                                                                                  >  2:46477/1‑88 (MQ=255)
                                                                                      tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                                       <  2:272060/99‑1 (MQ=255)
                                                                                      tattatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                                       >  1:272060/1‑99 (MQ=255)
                                                                                       attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                              <  2:113115/139‑1 (MQ=255)
                                                                                       attatCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                              <  1:979576/139‑1 (MQ=255)
                                                                                           tCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagt                                          >  1:691750/1‑139 (MQ=255)
                                                                                            caacaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagtaagta                                         <  1:340671/139‑1 (MQ=255)
                                                                                              acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGAt                                                                   <  1:839159/111‑1 (MQ=255)
                                                                                              acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGAt                                                                   >  2:839159/1‑111 (MQ=255)
                                                                                              acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGta                                                               <  1:887813/115‑1 (MQ=255)
                                                                                              acaaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGta                                                               >  2:887813/1‑115 (MQ=255)
                                                                                                aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                  <  2:64863/94‑1 (MQ=255)
                                                                                                aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCct                                                                                  >  1:64863/1‑94 (MQ=255)
                                                                                                aaGACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGt                                     >  2:803602/1‑139 (MQ=255)
                                                                                                  gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTAt                                                                                          >  1:1007344/1‑84 (MQ=255)
                                                                                                  gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTa                                                                            <  1:601603/98‑1 (MQ=255)
                                                                                                  gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGttt                                   >  1:164278/1‑139 (MQ=255)
                                                                                                  gACCAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCCCTCTTa                                                                            >  2:601603/1‑98 (MQ=255)
                                                                                                  gACCAGACACTTATAGAATTAGGTATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTAt                                                                                          <  2:1007344/84‑1 (MQ=255)
                                                                                                     cAGACACTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTAtt                                <  1:803602/139‑1 (MQ=255)
                                                                                                        acacTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGa                                                                    >  2:947953/1‑100 (MQ=255)
                                                                                                        acacTTATAGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGa                                                                    <  1:947953/100‑1 (MQ=255)
                                                                                                          acTTATAGAATTAGTTATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATc                           <  2:122060/139‑1 (MQ=255)
                                                                                                             tataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTaa                        >  1:88246/1‑139 (MQ=255)
                                                                                                              ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                                        <  2:844871/90‑1 (MQ=255)
                                                                                                              ataGAATTAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGc                                                                        >  1:844871/1‑90 (MQ=255)
                                                                                                                    ttAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATg                 >  2:97970/1‑139 (MQ=255)
                                                                                                                    ttAGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATg                 >  1:506990/1‑139 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                                       <  2:544961/67‑1 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCgg                                                                                       >  1:544961/1‑67 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                            >  1:969478/1‑94 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGtaata                                                            <  2:969478/94‑1 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                              <  2:552119/108‑1 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACagt                                              >  1:552119/1‑108 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCAt                  <  2:175965/136‑1 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCAt                  >  1:175965/1‑136 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt               >  2:728882/1‑139 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt               >  2:420686/1‑139 (MQ=255)
                                                                                                                      aGTGATGATTTATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGtt               >  2:868494/1‑139 (MQ=255)
                                                                                                                              tttATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACtt       <  2:664909/139‑1 (MQ=255)
                                                                                                                              tttATTCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACtt       <  2:164278/139‑1 (MQ=255)
                                                                                                                                   tCATTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTTTtata  <  1:97970/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GCTAACCTTCAAACGGGGTTAATCTTTGAACGTATGTCATGTTTCAACCCTTCAGATCGTGAATCTAAAGGGTTACATATTAACTATATTATCAACAAGACCAGACACTTATAGAATTAGTGATGATTTATTCACTAAAAAGGCCATGAACTGGAGGAAAACGATCTGTCCAGAAGGTCTATCGGTTCCTCTCTTATTGCCTGATTGTAATAATCGTTATTACAGTAAGTATTGTTTATTTCATCTAAGGTCATGTTTTTAACTTTTATA  >  NC_000913/1544250‑1544519

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