Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,227,983 C→T 6.8% E92K (GAA→AAA)  PP_5624 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,227,9830CT 6.8% 171.3 / 2.0 59E92K (GAA→AAA) PP_5624hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base T (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

ACAGCATGACGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAGGTTTTG  >  NC_002947/4227837‑4228125
                                                                                                                                                  |                                                                                                                                              
aCAGCATGACGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCtg                                                                                                                                              >  5:331987/1‑149 (MQ=255)
  aGCATGACGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCtgtg                                                                                                                                            >  3:133407/1‑149 (MQ=255)
  aGCATGACCAGCTTCCTGATCAGACCGAGTAGTAGCGAGTATCTACCCCGAAAGTTGGATCGCCATTATCCTCGTGCTCCAACTCCACATACCTCTAAGTCATCAGCTATAGGATATGCCTGACATCAGAACGATCATCACGATCtgtg                                                                                                                                            >  4:11317/1‑149 (MQ=255)
    cATGACGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGt                                                                                                                                          <  4:87731/149‑1 (MQ=255)
       gACGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCg                                                                                                                                        >  6:288234/1‑148 (MQ=255)
        aCGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCCGACATCAGATCGATCATCACGTTCTGTGGTCGaa                                                                                                                                      >  2:55069/1‑149 (MQ=255)
         cGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGaa                                                                                                                                      >  1:147088/1‑148 (MQ=255)
         cGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGaa                                                                                                                                      <  7:290164/148‑1 (MQ=255)
         cGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAAt                                                                                                                                     <  4:133407/149‑1 (MQ=255)
          gAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAAt                                                                                                                                     >  6:30219/1‑148 (MQ=255)
                ccTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTc                                                                                                                                <  7:78136/147‑1 (MQ=255)
                ccTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTc                                                                                                                                <  6:331987/147‑1 (MQ=255)
                ccTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCtt                                                                                                                              >  4:5334/1‑149 (MQ=255)
                  tGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGAAATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTgg                                                                                                                            >  2:77705/1‑149 (MQ=255)
                   gATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCTTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCtgt                                                                                                                                             <  6:243722/131‑1 (MQ=255)
                   gATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATAATCACGTTCtgt                                                                                                                                             >  5:243722/1‑131 (MQ=255)
                           cGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTCACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCt                                                                                                                   >  8:150437/1‑149 (MQ=255)
                             agtagtAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGCTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTgc                                                                                                                 <  6:56350/149‑1 (MQ=255)
                              gtagtaGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTgca                                                                                                                <  2:147088/149‑1 (MQ=255)
                                    gCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGAt                                                                                                          >  8:254695/1‑149 (MQ=255)
                                                cccGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTgg                                                                                                                            <  1:2963/119‑1 (MQ=255)
                                                cccGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTgg                                                                                                                            >  2:2963/1‑119 (MQ=255)
                                                  cGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTcc                                                                                            <  3:5334/149‑1 (MQ=255)
                                                               gCCATTATCGTCGTGCTCCAACTCCCCATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATAACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGcagctcagct                                                                               >  3:323619/1‑149 (MQ=255)
                                                                ccATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGcagctcagct                                                                               >  2:279549/1‑148 (MQ=255)
                                                                         tcgtGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCAtt                                                                     >  7:29866/1‑149 (MQ=255)
                                                                              cTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCa                                                                 <  5:155507/148‑1 (MQ=255)
                                                                               tCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTc                                                                                                                                >  2:68771/1‑84 (MQ=255)
                                                                               tCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCa                                                                 >  7:257932/1‑147 (MQ=255)
                                                                                 cAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCaga                                                               <  7:254695/147‑1 (MQ=255)
                                                                                    cTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTCACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCctc                                                          <  7:150437/149‑1 (MQ=255)
                                                                                     tCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCctct                                                         >  6:270257/1‑149 (MQ=255)
                                                                                       cacaTTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCctct                                                         <  8:75956/147‑1 (MQ=255)
                                                                                       cacaTTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCtctc                                                        <  3:11317/148‑1 (MQ=255)
                                                                                          aTTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTccc                                                      >  3:268264/1‑147 (MQ=255)
                                                                                           ttCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCAcgc                                                    >  7:313581/1‑148 (MQ=255)
                                                                                               tCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGAt                                                 >  5:324030/1‑147 (MQ=255)
                                                                                               tCTTAGTGATCAGCTATAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATaa                                               <  5:30219/149‑1 (MQ=255)
                                                                                                   aGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGaa                                            <  8:257932/148‑1 (MQ=255)
                                                                                                   aGTGATCAGCTCTAGGATATGGCCGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGaa                                            <  1:55069/148‑1 (MQ=255)
                                                                                                      gATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTACAAGCAGCTCAGATAAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCg                                        >  7:226306/1‑149 (MQ=255)
                                                                                                        tCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTTTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATaa                                               >  3:210297/1‑140 (MQ=255)
                                                                                                        tCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTTTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATaa                                               <  4:210297/140‑1 (MQ=255)
                                                                                                        tCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGa                                      <  5:270257/149‑1 (MQ=255)
                                                                                                               tAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGcagct                                                                                    <  5:262347/96‑1 (MQ=255)
                                                                                                               tAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGcagct                                                                                    >  6:262347/1‑96 (MQ=255)
                                                                                                               tAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCtc                                 <  6:324030/147‑1 (MQ=255)
                                                                                                                      tGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCt                        >  1:138710/1‑149 (MQ=255)
                                                                                                                             cATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGAcc                 <  4:323619/149‑1 (MQ=255)
                                                                                                                                cAGATCGATCATCACGTTTTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATaa                                               >  8:255114/1‑116 (MQ=255)
                                                                                                                                cAGATCGATCATCACGTTTTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATaa                                               <  7:255114/116‑1 (MQ=255)
                                                                                                                                 aGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCt             >  8:191099/1‑149 (MQ=255)
                                                                                                                                   atcgatcATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTg            >  2:130423/1‑148 (MQ=255)
                                                                                                                                     cgatcATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTgg           >  4:218910/1‑147 (MQ=255)
                                                                                                                                      gatcATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAg        <  8:29866/149‑1 (MQ=255)
                                                                                                                                      aactatcaCGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAg        <  1:218281/145‑1 (MQ=255)
                                                                                                                                       atcatcACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAgg       >  3:50579/1‑149 (MQ=255)
                                                                                                                                           tcaCGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAg        >  2:218281/1‑144 (MQ=255)
                                                                                                                                             aCGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAGGTTTTg  <  1:77705/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ACAGCATGACGAGCTTCCTGATCAGTCCGAGTAGTAGCGAGTATCTAGCCCGAAAGTTGGATCGCCATTATCGTCGTGCTCCAACTCCACATTCGTCTTAGTGATCAGCTCTAGGATATGGCTGACATCAGATCGATCATCACGTTCTGTGGTCGAATCAGTCTTGGTCAGAAGCTGCAGCAGATTGTTTTTTGACTCCAAGCAGCTCAGCTCAGACTCATTCGCAGAGCCTCTCCCGCGATAAGAACTCGGATCCTCATCATCGCTTTCGACCTTCTGGAAGGTTTTG  >  NC_002947/4227837‑4228125

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