Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,115 G→A *255* (TAG→TAA yabI → DedA family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,1150GA100.0% 226.3 / NA 69*255* (TAG→TAAyabIDedA family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/32);  total (37/32)

TCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAA  >  NC_000913/71987‑72246
                                                                                                                                |                                                                                                                                   
tCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTc                                                                                                                            <  1:1342893/138‑1 (MQ=255)
tCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTc                                                                                                                            >  2:1342893/1‑138 (MQ=255)
tCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCt                                                                                                                           >  2:171717/1‑139 (MQ=255)
tCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCt                                                                                                                           >  1:955167/1‑139 (MQ=255)
    gcgGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGt                                                                                                                             <  2:374654/133‑1 (MQ=255)
    gcgGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGt                                                                                                                             >  1:374654/1‑133 (MQ=255)
    gcgGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTtgtg                                                                                                                       >  2:168405/1‑139 (MQ=255)
      ggTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCg                                                                                                                     >  2:1581998/1‑139 (MQ=255)
          gTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGtt                                                                                                                 >  1:1415542/1‑139 (MQ=255)
                  tgGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATg                                                                                                         >  2:26062/1‑139 (MQ=255)
                  tgGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATg                                                                                                         >  1:895970/1‑139 (MQ=255)
                         cGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACg                                                                                                    <  1:1059856/137‑1 (MQ=255)
                         cGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACg                                                                                                    >  2:1059856/1‑137 (MQ=255)
                         cGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCt                                                                                                  >  1:485433/1‑139 (MQ=255)
                                 aTTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGt                                                                                          <  1:168405/139‑1 (MQ=255)
                                        ccATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                   <  2:718585/139‑1 (MQ=255)
                                                      gcgcTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGc                                                                     <  2:1415542/139‑1 (MQ=255)
                                                       cgcTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCa                                                                    <  2:955167/139‑1 (MQ=255)
                                                       cgcTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCa                                                                    <  1:500741/139‑1 (MQ=255)
                                                        gcTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCAc                                                                   <  2:198789/139‑1 (MQ=255)
                                                         cTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACt                                                                  >  1:1131269/1‑139 (MQ=255)
                                                         cTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACt                                                                  >  1:928859/1‑139 (MQ=255)
                                                            gtggtgTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGc                                                                                               <  1:1239694/107‑1 (MQ=255)
                                                            gtggtgTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGc                                                                                               >  2:1239694/1‑107 (MQ=255)
                                                                 gtTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGc                                                                                               >  1:408946/1‑102 (MQ=255)
                                                                 gtTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGc                                                                                               <  2:408946/102‑1 (MQ=255)
                                                                 gtTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAg                                                          >  1:926724/1‑139 (MQ=255)
                                                                   tAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGt                                                        <  1:26062/139‑1 (MQ=255)
                                                                         gCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATa                                                  >  1:572836/1‑139 (MQ=255)
                                                                          ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACg                                                                                                    >  1:1445208/1‑88 (MQ=255)
                                                                          ccACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACg                                                                                                    <  2:1445208/88‑1 (MQ=255)
                                                                           cACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAg                                                <  2:895970/139‑1 (MQ=255)
                                                                            aCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAgg                                               >  2:912536/1‑139 (MQ=255)
                                                                                 ttGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTc                                          >  1:1031229/1‑139 (MQ=255)
                                                                                          gtgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTtgtg                                                                                                                       <  1:343168/53‑1 (MQ=255)
                                                                                          gtgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTtgtg                                                                                                                       >  2:343168/1‑53 (MQ=255)
                                                                                          gtgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCcgcg                                                                                           <  1:802991/81‑1 (MQ=255)
                                                                                          gtgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCcgcg                                                                                           >  2:802991/1‑81 (MQ=255)
                                                                                           tgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCcgc                                                                                            >  2:1607402/1‑79 (MQ=255)
                                                                                           tgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCcgc                                                                                            <  1:1607402/79‑1 (MQ=255)
                                                                                           tgtATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATc                                >  1:593281/1‑139 (MQ=255)
                                                                                              atatCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                             <  2:1131269/139‑1 (MQ=255)
                                                                                              atatCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                             <  1:437042/139‑1 (MQ=255)
                                                                                               tatCGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                            >  1:294725/1‑139 (MQ=255)
                                                                                                  cGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGt                         <  2:572836/139‑1 (MQ=255)
                                                                                                  cGATATTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGt                         <  1:912536/139‑1 (MQ=255)
                                                                                                     tatTTTGCGTAAAGTGGTTGGGGTTTAAGCGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCa                      >  2:45969/1‑139 (MQ=255)
                                                                                                     tatTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAgg                                                         <  2:16839/104‑1 (MQ=255)
                                                                                                     tatTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAgg                                                         >  1:16839/1‑104 (MQ=255)
                                                                                                     tatTTTGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCa                      >  2:1514461/1‑139 (MQ=255)
                                                                                                       ttttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGc                                                                                                        <  1:149699/55‑1 (MQ=255)
                                                                                                       ttttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGc                                                                                                        >  2:149699/1‑55 (MQ=255)
                                                                                                       ttttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                  <  2:209536/77‑1 (MQ=255)
                                                                                                       ttttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                  >  1:209536/1‑77 (MQ=255)
                                                                                                       ttttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATa                    >  1:985448/1‑139 (MQ=255)
                                                                                                         ttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaaaa                                                                        >  1:653162/1‑85 (MQ=255)
                                                                                                         ttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaaaa                                                                        <  2:653162/85‑1 (MQ=255)
                                                                                                         ttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaa                                                                          <  1:904648/83‑1 (MQ=255)
                                                                                                         ttGCGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACaa                                                                          >  2:904648/1‑83 (MQ=255)
                                                                                                            cGTAAAGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAAc               >  1:549642/1‑139 (MQ=255)
                                                                                                               aaaGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCg            <  1:171717/139‑1 (MQ=255)
                                                                                                                 aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACa                                                                           >  1:322283/1‑74 (MQ=255)
                                                                                                                 aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACa                                                                           <  2:322283/74‑1 (MQ=255)
                                                                                                                 aGTGGTTGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTa          <  2:1031229/139‑1 (MQ=255)
                                                                                                                      ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                     >  2:602038/1‑59 (MQ=255)
                                                                                                                      ttGGGGTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                     <  1:602038/59‑1 (MQ=255)
                                                                                                                         gggTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATATTAACCCGTAATCCCCaa  <  2:1589072/139‑1 (MQ=255)
                                                                                                                           gTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa                           <  2:640164/112‑1 (MQ=255)
                                                                                                                           gTTTAAGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACa                           >  1:640164/1‑112 (MQ=255)
                                                                                                                                |                                                                                                                                   
TCCCGCGGTCGTTTGTTGTGGCTGACGCCGTTGATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAA  >  NC_000913/71987‑72246

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: