Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,143,641 G→A H243Y (CAT→TAT)  rne ← endoribonuclease; RNA‑binding protein;RNA degradosome binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,143,6410GA98.8% 277.4 / ‑4.3 81H243Y (CAT→TAT) rneendoribonuclease; RNA‑binding protein;RNA degradosome binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (42/38);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCT  >  NC_000913/1143544‑1143737
                                                                                                 |                                                                                                
aTCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACt                                                                                                <  1:229004/100‑1 (MQ=255)
       aCTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTg                                                                                          >  2:1133493/1‑99 (MQ=255)
         tctaTCTGGTAGCGGCTGACCAGCGGGAGCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGcc                                                                                        <  1:224550/96‑1 (MQ=255)
          cGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGCGCGACCTAATGCAGCGATATACTGACGTGCCAg                                                                                      <  1:1093890/100‑1 (MQ=255)
            aTCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGtt                                                                                    >  1:1481429/1‑100 (MQ=255)
              cTGGTAGTGGCTTAACAGCGGGATCTCGCCTGGGTAAAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGACAGtt                                                                                    <  1:520130/98‑1 (MQ=255)
              cTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTc                                                                                   >  1:1141494/1‑99 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTGGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                                <  2:1190266/100‑1 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGa                                                                                 >  1:1572822/1‑99 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                                <  1:1202728/100‑1 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                                <  1:1087259/100‑1 (MQ=255)
                ggTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                                <  2:1410944/100‑1 (MQ=255)
                 gTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGc                                                                               <  2:1543874/100‑1 (MQ=255)
                 gTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGc                                                                               <  2:479963/100‑1 (MQ=255)
                 gTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGc                                                                               >  2:1794923/1‑100 (MQ=255)
                  tAGTGGCTGAACAGCGGGATCTCGTCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAg                                                                                <  2:392321/98‑1 (MQ=255)
                   aGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCAc                                                                             >  1:337059/1‑100 (MQ=255)
                      ggCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACttt                                                                          <  1:1221333/100‑1 (MQ=255)
                        cTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTACTGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCg                                                                        >  1:1166375/1‑100 (MQ=255)
                          gAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCggg                                                                      <  2:714856/100‑1 (MQ=255)
                           aaCAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                     >  2:1420597/1‑100 (MQ=255)
                            aCAGCGGGAATTCGCCGGTGTACAGTTTGTTTTTGCTGCTGAAAGCCGGGCGACGTAATACAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                     <  1:478895/99‑1 (MQ=255)
                             cAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGt                                                                     <  1:1603308/98‑1 (MQ=255)
                              aGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTAt                                                                  <  1:1115619/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                 <  2:960567/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                 <  1:979531/100‑1 (MQ=255)
                               gCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                 >  2:1579800/1‑100 (MQ=255)
                               gCGGGATCTCGCCGGCGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATc                                                                 <  2:1793430/100‑1 (MQ=255)
                                cGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCg                                                                <  2:1077575/100‑1 (MQ=255)
                                  gtaTCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATAAACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCAAt                                                              >  2:1684361/3‑100 (MQ=255)
                                   gATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATa                                                             <  2:1141505/100‑1 (MQ=255)
                                   gATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATa                                                             <  1:31415/100‑1 (MQ=255)
                                    aTCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATaa                                                            <  1:1265619/100‑1 (MQ=255)
                                     tctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAg                                                           <  1:1667577/100‑1 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCATTTCGAGCACTTTCGGGTTATCGATaa                                                            >  1:57050/1‑98 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAgg                                                          >  2:119475/1‑100 (MQ=255)
                                      ctcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAgg                                                          >  2:705374/1‑100 (MQ=255)
                                       tcGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGa                                                         >  1:94230/1‑100 (MQ=255)
                                           cGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTc                                                     <  1:1796428/100‑1 (MQ=255)
                                            ggTGTACAGCTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCg                                                    >  1:1701597/1‑100 (MQ=255)
                                               gtACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCg                                                 <  2:67488/100‑1 (MQ=255)
                                                tACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGa                                                <  2:1392456/100‑1 (MQ=255)
                                                   aGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATg                                              <  1:1077563/99‑1 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                            >  1:53406/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                            >  1:1292464/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTc                                            >  2:1127745/1‑100 (MQ=255)
                                                    gTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCAGGTTATCGATAAGGATTTCGCCGATGTc                                            >  2:92658/1‑100 (MQ=255)
                                                       tGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCAGATGTCCTg                                         >  2:170972/1‑100 (MQ=255)
                                                        gATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGa                                        <  1:284957/100‑1 (MQ=255)
                                                         aTTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGAc                                       >  1:1115810/1‑100 (MQ=255)
                                                          ttttGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGAc                                       >  1:845462/1‑99 (MQ=255)
                                                              gctgctGAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaag                                  <  2:352771/100‑1 (MQ=255)
                                                               ctgctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGATATCGATAAGGATTTCGCCGATGTCCTGACgtaagt                                 >  2:97547/1‑100 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagtaa                               >  1:107564/1‑100 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagtaa                               >  1:582049/1‑100 (MQ=255)
                                                                 gctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACgtaagta                                >  1:1433406/1‑99 (MQ=255)
                                                                  ctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAAt                              >  1:876861/1‑100 (MQ=255)
                                                                  ctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAAt                              >  1:1290879/1‑100 (MQ=255)
                                                                  ctgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCAAATGTCCTGAAgtaagaaa                               >  2:666344/1‑99 (MQ=255)
                                                                   tgAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTCCGCCGATGTCCTGACGTAAGTAATc                             >  1:927438/1‑100 (MQ=255)
                                                                    gAAATCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATcg                            >  1:249605/1‑100 (MQ=255)
                                                                      naTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGGACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATcgcg                          >  1:1007560/2‑100 (MQ=255)
                                                                       aTCCGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCgg                         >  2:593247/1‑100 (MQ=255)
                                                                          cGGGCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAAt                      >  1:596938/1‑100 (MQ=255)
                                                                            ggCGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAAt                      >  1:1169529/1‑98 (MQ=255)
                                                                              cGACCTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTACTGACGTAAATAATCGCGGaa                       >  2:116627/1‑95 (MQ=255)
                                                                                 ccTAATGCAGCGATATACTGGCGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGaa               >  1:1107189/1‑100 (MQ=255)
                                                                                  cTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAAc              <  1:1420586/100‑1 (MQ=255)
                                                                                  cTAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAAc              <  2:61959/100‑1 (MQ=255)
                                                                                   tAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACg             <  1:1746992/100‑1 (MQ=255)
                                                                                   tAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACg             <  2:1290890/100‑1 (MQ=255)
                                                                                   tAATGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACg             >  2:1136627/1‑100 (MQ=255)
                                                                                    aaTGCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGa            >  1:1561748/1‑100 (MQ=255)
                                                                                      tGCAGCGATATATTGACGTACCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGCTGTCCTGACGTAAGTAATCGCGGACTGCGCGAACGAtt          >  1:695232/1‑100 (MQ=255)
                                                                                       gCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa         >  2:245310/1‑100 (MQ=255)
                                                                                       gCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa         >  1:27532/1‑100 (MQ=255)
                                                                                       gCAGCGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTa         <  1:119475/100‑1 (MQ=255)
                                                                                           cGATATACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGtt     >  1:392042/1‑100 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc   <  1:1579789/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc   <  2:1701610/100‑1 (MQ=255)
                                                                                             atatACTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGc   <  2:1445446/100‑1 (MQ=255)
                                                                                              tataCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGct  <  2:1292475/100‑1 (MQ=255)
                                                                                                 |                                                                                                
ATCTGTGACTCGATCTGGTAGTGGCTGAACAGCGGGATCTCGCCGGTGTACAGTTTGATTTTGCTGCTGAAATCCGGGCGACCTAATGCAGCGATATGCTGACGTGCCAGTTCGAGCACTTTCGGGTTATCGATAAGGATTTCGCCGATGTCCTGACGTAAGTAATCGCGGAATGCGCGAACGATTACGTTGCT  >  NC_000913/1143544‑1143737

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