Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,900 G→C 8.9% P232A (CCA→GCA)  trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,9000GC 8.9% 167.5 / 6.4 56P232A (CCA→GCA) trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (3/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GCCGTTCAGGCTGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAC  >  NC_002947/92751‑93038
                                                                                                                                                     |                                                                                                                                          
gccgTTCAGGCTGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTgg                                                                                                                                            <  1:50237/150‑1 (MQ=255)
    atCAGGCTGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATTAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGa                                                                                                                                        >  2:193899/2‑150 (MQ=255)
           tGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                 >  2:855182/1‑149 (MQ=255)
            ggcggcATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctgg                                                                                                                                <  1:855182/150‑3 (MQ=255)
             gcggcATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGNAGCCGCNTCCGACGCAGGCNACCAGGCTGTCTGGCAGACGGccttcc                                                                                                                               <  1:809211/150‑1 (MQ=255)
                   tGCTTGTCGGTGTGTCAGCCGTGGCCACCGGCTTTAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAGTGCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTctct                                                                                                                         <  1:222846/150‑1 (MQ=255)
                        gTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCaa                                                                                                                    <  1:777695/150‑1 (MQ=255)
                            gtgtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGGGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTTTTGCAACTgg                                                                                                                >  1:620760/1‑150 (MQ=255)
                            gcgtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTgg                                                                                                                <  2:481741/148‑1 (MQ=255)
                             tggctACGCCGTGGCCACCGTCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTggg                                                                                                               <  2:620760/146‑1 (MQ=255)
                              gtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGc                                                                                                              <  2:752062/150‑1 (MQ=255)
                              gtgtACGCCGTGGCCACCGGCTTCAACGCCGATGAGCTGCACGCTCGGCTCTTCGAGGAACTAATGACAAAGGCCCATGGCATTGCAGCAACGGCGAAAACAAACAACTATCCGATCTGCCGTTCGACCTTCCTTCTCTTGCAACTGGGc                                                                                                              >  2:1243612/1‑150 (MQ=255)
                                 tACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTGCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCCGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACg                                                                                                           <  1:1125787/150‑1 (MQ=255)
                                  aCGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACgg                                                                                                          >  2:438789/1‑150 (MQ=255)
                                           ccACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTg                                                                                                                 <  1:425881/134‑1 (MQ=255)
                                           ccACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTg                                                                                                                 >  2:425881/1‑134 (MQ=255)
                                            cACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGc                                                                                                <  1:1243612/150‑1 (MQ=255)
                                                       aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGAc                                                                                     >  1:813565/1‑150 (MQ=255)
                                                       aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGGATCTTGGCCGTAGATCGAc                                                                                     >  2:582212/1‑150 (MQ=255)
                                                         gCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTg                                                                                   <  2:86979/150‑1 (MQ=255)
                                                            gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTctct                                                                                                                         <  2:1143607/109‑1 (MQ=255)
                                                            gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTctct                                                                                                                         >  1:1143607/1‑109 (MQ=255)
                                                            gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGaa                                                                                >  2:742665/1‑150 (MQ=255)
                                                              tgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGccttcct                                                                                                                              >  1:923050/1‑102 (MQ=255)
                                                              tgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGccttcct                                                                                                                              <  2:923050/102‑1 (MQ=255)
                                                                  cTGCGCGCTCTGCTCCGCGAGGAGCTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTcgcg                                                                          <  2:1036246/150‑1 (MQ=255)
                                                                   tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCCGGCTCTCTGGCAGACGGCCGTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCgcgc                                                                         <  2:247231/150‑1 (MQ=255)
                                                                   tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCgcgc                                                                         >  1:256852/1‑150 (MQ=255)
                                                                   tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                 <  1:57182/94‑2 (MQ=255)
                                                                   tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                 >  2:57182/1‑93 (MQ=255)
                                                                     cACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCCTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCAc                                                                       <  2:172995/150‑1 (MQ=255)
                                                                              ctctTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCggg                                                              >  2:850982/1‑150 (MQ=255)
                                                                               tctTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCgggg                                                             >  2:413020/1‑150 (MQ=255)
                                                                                      ggAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTctct                                                                                                    >  2:353009/1‑104 (MQ=255)
                                                                                      ggAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTctct                                                                                                    <  1:353009/104‑1 (MQ=255)
                                                                                       gcaCTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGAACACCATGGCGGGGTACGGgtg                                                     <  1:39133/148‑1 (MQ=255)
                                                                                          cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                  <  2:205836/150‑1 (MQ=255)
                                                                                          cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                  <  1:413020/150‑1 (MQ=255)
                                                                                           tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGgtgt                                                    >  2:1112363/1‑147 (MQ=255)
                                                                                           tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGgtgt                                                    <  1:1112363/147‑1 (MQ=255)
                                                                                                    aGGCCCATGGCAGTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTTTCTTGCAACCGGGCAAGGTTCTCTTTGCCGATGATCCACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACg                                        <  1:582212/150‑1 (MQ=255)
                                                                                                        ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                    >  1:1157989/1‑150 (MQ=255)
                                                                                                        ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                    <  2:256852/150‑1 (MQ=255)
                                                                                                        ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGAATGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                    >  2:1129128/1‑150 (MQ=255)
                                                                                                         cATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGcc                                   >  1:839401/1‑150 (MQ=255)
                                                                                                         cATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGcc                                   >  2:46516/1‑150 (MQ=255)
                                                                                                          aTGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCg                                  >  1:692033/1‑150 (MQ=255)
                                                                                                           tGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGa                                 >  1:746587/1‑150 (MQ=255)
                                                                                                               aTTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAg                             <  1:46516/150‑1 (MQ=255)
                                                                                                               aTTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAg                             >  1:805602/1‑150 (MQ=255)
                                                                                                                     gccgccgccGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGaa                       <  1:193899/150‑1 (MQ=255)
                                                                                                                          cgccgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGgtgt                  >  1:272994/1‑150 (MQ=255)
                                                                                                                              gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTg                                     >  2:1156516/1‑127 (MQ=255)
                                                                                                                              gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTg                                     <  1:1156516/127‑1 (MQ=255)
                                                                                                                                     gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg       <  1:438789/150‑1 (MQ=255)
                                                                                                                                     gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg       <  2:839401/150‑1 (MQ=255)
                                                                                                                                          catccTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAc  >  1:887627/5‑150 (MQ=255)
                                                                                                                                                     gCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                        >  1:948275/1‑85 (MQ=255)
                                                                                                                                                     gCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                        <  2:948275/85‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
GCCGTTCAGGCTGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAC  >  NC_002947/92751‑93038

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: