Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,077,580 C→T 6.7% G435S (GGC→AGC)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,077,5800CT 6.7% 183.3 / 2.0 59G435S (GGC→AGC) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGGGTAGATTAGATTCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGA  >  NC_000913/1077434‑1077723
                                                                                                                                                  |                                                                                                                                               
cTGCGGGTAGATTAGATTCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCtt                                                                                                                                               <  2:384485/149‑1 (MQ=255)
    gggTAGATTAGATTCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAACCCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCAc                                                                                                                                           >  4:263616/1‑149 (MQ=255)
       tagattagatGCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGGCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGGCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAg                                                                                                                                        <  4:61881/149‑1 (MQ=255)
         gattagattCGGCAGCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCAGCTGCTCACGCTTAATCTCACTATCCCAGTACGCGCCTTTCACCAGGc                                                                                                                                      <  1:177151/149‑1 (MQ=255)
            tagattCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGaa                                                                                                                                   >  7:91389/1‑149 (MQ=255)
              gattCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCAc                                                                                                                                           <  1:252642/139‑1 (MQ=255)
              gattCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCAc                                                                                                                                           >  2:252642/1‑139 (MQ=255)
                   ggCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAg                                                                                                                            >  2:74494/1‑149 (MQ=255)
                            agcagcTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCCCCAGGCGAATCATCAgacggcgacg                                                                                                                   >  5:158079/1‑149 (MQ=255)
                                  tttttCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCAATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCg                                                                                                             >  6:305852/1‑149 (MQ=255)
                                   ttttCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGc                                                                                                                                      <  4:28219/123‑1 (MQ=255)
                                   ttttCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGc                                                                                                                                      >  3:28219/1‑123 (MQ=255)
                                    tttCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCgg                                                                                                            <  1:74494/148‑1 (MQ=255)
                                       cGCACAGGCGAGATCAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCCGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTg                                                                                                         <  8:18943/148‑1 (MQ=255)
                                        gCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCGTCAAGGCCGTCCATCTGCGCAGGCGTAATTTCACTATCCCAGTACGCGCCTGTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGc                                                                                                       >  4:4643/1‑149 (MQ=255)
                                          acaGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGc                                                                                                       >  5:188962/1‑147 (MQ=255)
                                          acaGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCga                                                                                                     <  4:89448/149‑1 (MQ=255)
                                             ggCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGat                                                                                                  <  8:228651/149‑1 (MQ=255)
                                                   aTAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCGGACGGCGACGGCTGCGGGTGGCGAGatcaat                                                                                              >  3:239833/1‑147 (MQ=255)
                                                   aTAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTg                                                                                                         <  1:249404/136‑1 (MQ=255)
                                                   aTAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTg                                                                                                         >  2:249404/1‑136 (MQ=255)
                                                     aaGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCg                                                                                                             >  3:94021/1‑130 (MQ=255)
                                                     aaGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCg                                                                                                             <  4:94021/130‑1 (MQ=255)
                                                       gAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAGTCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGt                                                                                         >  6:206825/1‑148 (MQ=255)
                                                           cGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAatcat                                                                                                                               <  6:311043/106‑1 (MQ=255)
                                                           cGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAatcat                                                                                                                               >  5:311043/1‑106 (MQ=255)
                                                             tCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGAt                                                                                  <  8:91389/149‑1 (MQ=255)
                                                             tCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGAt                                                                                  >  1:152524/1‑149 (MQ=255)
                                                                gTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCAc                                                                               >  3:363395/1‑149 (MQ=255)
                                                                   taCACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTCCGCGCTTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCaa                                                                            >  4:23920/1‑149 (MQ=255)
                                                                   taCACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCTTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCaa                                                                            >  4:29034/1‑149 (MQ=255)
                                                                    acacCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAAc                                                                           >  7:367924/1‑149 (MQ=255)
                                                                     cacCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACg                                                                          <  3:171124/149‑1 (MQ=255)
                                                                              ggTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGtt                                                                 >  8:191221/1‑149 (MQ=255)
                                                                              ggTATAAACCGGATAACCGTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGtt                                                                 >  6:169843/1‑149 (MQ=255)
                                                                                    aaCCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGa                                                            >  3:54675/1‑148 (MQ=255)
                                                                                     aCCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCGGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATa                                                          <  4:239833/149‑1 (MQ=255)
                                                                                       cGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCCCAAACGGGCAGCGTTTTTGAAAAg                                                        >  4:3643/1‑149 (MQ=255)
                                                                                             aCCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGaa                                                  <  3:263616/149‑1 (MQ=255)
                                                                                              ccTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAAt                                                 <  6:188962/149‑1 (MQ=255)
                                                                                               cTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGTGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCGCCAACGGGCAGCGTTTTTGATAAGCCTGAAt                                                 >  1:50397/1‑148 (MQ=255)
                                                                                                ttCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATaa                                               >  5:184459/1‑149 (MQ=255)
                                                                                                 tCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCCGTACGCGCCTTTCACCCGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACTGGCAGCGTTTTTGATAAGCCTGAaaaa                                               >  6:126995/1‑148 (MQ=255)
                                                                                                      gCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAAc                                         <  6:158079/149‑1 (MQ=255)
                                                                                                       ccGGCCAGCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAAcc                                        <  8:277726/149‑1 (MQ=255)
                                                                                                         gTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGa                                      >  8:4941/1‑149 (MQ=255)
                                                                                                          tCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGAt                                     <  6:184459/149‑1 (MQ=255)
                                                                                                          gccATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGAt                                     <  8:297461/148‑1 (MQ=255)
                                                                                                            cATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGCCGGCGACGGCTGCGGGTGGCGAGATCAATCTGGGAATCGATCACCCCCGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGc                                   >  4:228410/1‑149 (MQ=255)
                                                                                                            cATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATg                                    <  8:158060/148‑1 (MQ=255)
                                                                                                            cATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGc                                   <  3:3643/149‑1 (MQ=255)
                                                                                                                   gcACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCa                            <  7:191221/149‑1 (MQ=255)
                                                                                                                     aCGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGc                          <  4:54675/149‑1 (MQ=255)
                                                                                                                        cTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGTGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCGCCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTg                       <  2:50397/149‑1 (MQ=255)
                                                                                                                         ttAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGc                      >  2:249415/1‑149 (MQ=255)
                                                                                                                                     tCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGc            <  5:206825/147‑1 (MQ=255)
                                                                                                                                       ccAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCg         >  8:221621/1‑148 (MQ=255)
                                                                                                                                             cgcgCTTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACaga  <  3:23920/149‑1 (MQ=255)
                                                                                                                                             cgcgCTTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACaga  <  3:29034/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CTGCGGGTAGATTAGATTCGGCACCGCCAGCAGCTTTTTCGCACAGGCGAGATAAGAAACGTCGGTATACACCTTGCGGGTATAAACCGGATAACCTTCAAGGCCGTCCATCTGCGCACGCTTAATTTCACTATCCCAGTACGCGCCTTTCACCAGGCGAATCATCAGACGGCGACGGCTGCGGGTGGCGAGATCAATCAGGTAATCGATCACCAACGGGCAGCGTTTTTGATAAGCCTGAATAACAAAACCGATGCCGTTCCAGCCTGCCAGTTCCGGCTCGAAACAGA  >  NC_000913/1077434‑1077723

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