Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,208,517 C→T P129L (CCT→CTT) 
*201* (TAG→TAA
tfaP →
tfaE ←
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,5170CT100.0% 209.9 / NA 70P129L (CCT→CTT) 
*201* (TAG→TAA
tfaP
tfaE
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/32);  total (38/32)

TAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCT  >  NC_000913/1208383‑1208647
                                                                                                                                      |                                                                                                                                  
tAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACg                                                                                                                                >  2:310529/1‑139 (MQ=255)
        aGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGc                                                                                                                        >  2:1338467/1‑139 (MQ=255)
          tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAg                                                                                                                      >  1:797978/1‑139 (MQ=255)
          tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAg                                                                                                                      <  1:581232/139‑1 (MQ=255)
                ggTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttaggc                                                                                                                >  2:944056/1‑135 (MQ=255)
                    aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGa                                                                                                           >  2:782034/1‑139 (MQ=255)
                    aGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGa                                                                                                           >  2:1247669/1‑139 (MQ=255)
                      gCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATg                                                                                                         >  2:775108/1‑139 (MQ=255)
                                gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGCATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               <  2:1052700/139‑1 (MQ=255)
                                gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               <  2:797978/139‑1 (MQ=255)
                                gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               <  1:1029823/139‑1 (MQ=255)
                                gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCACATGACATCCGgcg                                                                                               <  1:1081797/139‑1 (MQ=255)
                                    aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTCTCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                           >  1:1144100/1‑139 (MQ=255)
                                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCgg                                                                                                                         <  2:560261/100‑1 (MQ=255)
                                              gaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCgg                                                                                                                         >  1:560261/1‑100 (MQ=255)
                                                    cTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGcc                                                                           >  2:492471/1‑139 (MQ=255)
                                                             tataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcgcg                                                                  >  2:433467/1‑139 (MQ=255)
                                                              ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                 <  1:1338467/139‑1 (MQ=255)
                                                                 aTTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                              >  2:576431/1‑139 (MQ=255)
                                                                   ttGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                            <  1:782034/139‑1 (MQ=255)
                                                                    tgtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                            <  2:666542/138‑1 (MQ=255)
                                                                    tgtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                            >  1:666542/1‑138 (MQ=255)
                                                                     gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGAt                                                                                                          <  1:464705/91‑1 (MQ=255)
                                                                     gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGAt                                                                                                          >  2:464705/1‑91 (MQ=255)
                                                                     gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTa                                                          >  1:1184226/1‑139 (MQ=255)
                                                                       ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                        >  1:506403/1‑139 (MQ=255)
                                                                        gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                       >  1:771894/1‑139 (MQ=255)
                                                                            gATTATCTGGAAGCACTGGAACGGGTTGCTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGTCATCCGGCGCGGTGCTGGTATCAGTTGCCGGCACCGTGTCAGTGTAATCCGGc                                                   >  1:959746/1‑139 (MQ=14)
                                                                                      aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTcc                                                                                                                           <  1:178232/58‑1 (MQ=255)
                                                                                      aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTcc                                                                                                                           >  2:178232/1‑58 (MQ=255)
                                                                                      aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgc                                                                                              <  1:596524/86‑1 (MQ=255)
                                                                                      aaGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgc                                                                                              >  2:596524/1‑86 (MQ=255)
                                                                                           cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               >  2:1424174/1‑80 (MQ=17)
                                                                                           cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               <  1:1424174/80‑1 (MQ=17)
                                                                                           cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTg                                    >  2:1120033/1‑139 (MQ=255)
                                                                                            tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                   <  1:310529/139‑1 (MQ=255)
                                                                                             ggAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                                     >  1:1546418/1‑57 (MQ=255)
                                                                                             ggAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                                     <  2:1546418/57‑1 (MQ=255)
                                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgc                                                                                              >  2:556989/1‑72 (MQ=17)
                                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcgc                                                                                              <  1:556989/72‑1 (MQ=17)
                                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               >  1:1454587/1‑71 (MQ=17)
                                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               <  2:1454587/71‑1 (MQ=17)
                                                                                                    gTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                           >  2:521919/1‑139 (MQ=255)
                                                                                                      tGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                         <  1:775108/139‑1 (MQ=255)
                                                                                                       gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                        <  1:944056/139‑1 (MQ=255)
                                                                                                       gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                        <  1:1247669/139‑1 (MQ=255)
                                                                                                       gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                        <  1:492471/139‑1 (MQ=255)
                                                                                                        aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCAc                                                                      >  2:1611578/1‑92 (MQ=255)
                                                                                                        aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCAc                                                                      <  1:1611578/92‑1 (MQ=255)
                                                                                                        aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcct                       >  1:212230/1‑139 (MQ=255)
                                                                                                            aTCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                   <  1:1120033/139‑1 (MQ=255)
                                                                                                             tCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTc                  <  2:1314899/139‑1 (MQ=255)
                                                                                                              ccggTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                           <  2:1257628/126‑1 (MQ=255)
                                                                                                              ccAGTGCGCCAGATATTGAATTGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                        >  1:1137487/1‑100 (MQ=11)
                                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTact                                                        <  2:1137487/100‑3 (MQ=11)
                                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                           >  1:1257628/1‑129 (MQ=255)
                                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                   <  2:631593/137‑1 (MQ=255)
                                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGt                   >  1:631593/1‑137 (MQ=255)
                                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACAACCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                            >  1:549125/1‑96 (MQ=11)
                                                                                                              ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACAACCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATg                                                            <  2:549125/96‑1 (MQ=11)
                                                                                                                  tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                                     <  2:1413165/36‑1 (MQ=14)
                                                                                                                  tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                                     >  1:1413165/1‑36 (MQ=14)
                                                                                                                  tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               <  2:1068235/57‑1 (MQ=11)
                                                                                                                  tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                               >  1:1068235/1‑57 (MQ=11)
                                                                                                                  tGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt             >  2:1458507/1‑139 (MQ=255)
                                                                                                                   gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTa            >  1:1515427/1‑139 (MQ=255)
                                                                                                                      ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGtt              <  1:1225570/134‑1 (MQ=255)
                                                                                                                      ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGtt              >  2:1225570/1‑134 (MQ=255)
                                                                                                                          atatTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCgg     >  1:821589/1‑139 (MQ=255)
                                                                                                                             ttGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCt  <  2:328470/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
TAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCT  >  NC_000913/1208383‑1208647

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