Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,822,898 G→C 9.3% A157A (GCG→GCC curA → NADPH‑dependent curcumin/dihydrocurcumin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,822,8980GC 9.3% 156.4 / 7.7 54A157A (GCG→GCCcurANADPH‑dependent curcumin/dihydrocurcumin reductase
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base C (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

CGATGGAAGCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCGGA  >  NC_002947/2822751‑2823044
                                                                                                                                                   |                                                                                                                                                  
cGATGGAAGCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTggcggcgg                                                                                                                                                   >  6:213895/1‑149 (MQ=255)
 gATGGAAGCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACACTGGTGGTggcggcggc                                                                                                                                                  >  8:263517/1‑149 (MQ=255)
      aaGCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCg                                                                                                                                             >  6:157235/1‑149 (MQ=255)
       aGCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCgg                                                                                                                                            <  2:68183/149‑1 (MQ=255)
        gCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCggg                                                                                                                                           >  8:213390/1‑149 (MQ=255)
                 ccGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgc                                                                                                                                  >  7:133232/1‑149 (MQ=255)
                  cgcgTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcg                                                                                                                                 <  7:105790/149‑1 (MQ=255)
                                      cgcATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGc                                                                                                             <  3:184615/149‑1 (MQ=255)
                                       gcATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTAACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCCCCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGc                                                                                                             >  4:263328/1‑148 (MQ=255)
                                         aTCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGc                                                                                                             >  5:152314/1‑146 (MQ=255)
                                             ttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCTCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCaa                                                                                                           <  5:244347/144‑1 (MQ=255)
                                             ttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCTCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCaa                                                                                                           >  6:244347/1‑144 (MQ=255)
                                             ttcttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATc                                                                                                       >  4:290468/1‑148 (MQ=255)
                                               cttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCaa                                                                                                     <  3:263328/148‑1 (MQ=255)
                                                ttATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAgg                                                                                                   <  3:21107/149‑1 (MQ=255)
                                                  aTGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcgc                                                                                                                                >  7:237395/1‑118 (MQ=255)
                                                  aTGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcgc                                                                                                                                <  8:237395/118‑1 (MQ=255)
                                                  aTGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAggg                                                                                                  <  8:166849/148‑1 (MQ=255)
                                                  aTGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAg                                                                                                    <  5:213895/146‑1 (MQ=255)
                                                        tGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACtgtg                                                                                                                         <  5:66047/119‑1 (MQ=255)
                                                        tGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACtgtg                                                                                                                         >  6:66047/1‑119 (MQ=255)
                                                           gCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAccacca                                                                                                                                          <  4:241296/99‑5 (MQ=255)
                                                           gCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAccacca                                                                                                                                          >  3:241296/1‑95 (MQ=255)
                                                              tCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCTACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCAtg                                                                                          <  1:55951/144‑1 (MQ=255)
                                                              tCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCAtg                                                                                          >  2:55951/1‑144 (MQ=255)
                                                                tttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCgg                                                                                   >  8:67868/1‑149 (MQ=255)
                                                                tttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAc                                                                                                                                               >  7:265450/1‑89 (MQ=255)
                                                                tttGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCCGCCGCCAc                                                                                                                                               <  8:265450/89‑1 (MQ=255)
                                                                  tGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCgg                                                                                   <  7:13066/147‑1 (MQ=255)
                                                                       tGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcgc                                                                                                                                >  5:272537/1‑97 (MQ=255)
                                                                       tGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGgcgc                                                                                                                                <  6:272537/97‑1 (MQ=255)
                                                                                tCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGACATGAGGGCGGTGTCGCCGGTGGcgcg                                                                   >  4:83708/1‑149 (MQ=255)
                                                                                           atGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCg                                                        >  1:162208/1‑149 (MQ=255)
                                                                                              ggCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGc                                                      <  7:213390/148‑1 (MQ=255)
                                                                                                 ttGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGc                                                  >  6:1199/1‑149 (MQ=255)
                                                                                                    cTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATAGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGt                                               >  7:204597/1‑149 (MQ=255)
                                                                                                       gATATTTGCCGGCCGCAGGCGGGTGAAACACTGGTGGTGGCGGCGGACACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCAAGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGAgggc                                              >  6:102383/1‑145 (MQ=255)
                                                                                                       gATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTc                                              <  3:290468/147‑1 (MQ=255)
                                                                                                        atatTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAg                                           <  6:152314/149‑1 (MQ=255)
                                                                                                         tatTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAgg                                          >  1:2241/1‑149 (MQ=255)
                                                                                                           ttGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGt                                         <  5:1199/148‑1 (MQ=255)
                                                                                                           ttGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTg                                        <  2:162208/149‑1 (MQ=255)
                                                                                                              gccggccgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTg                                     >  5:114706/1‑149 (MQ=255)
                                                                                                                    cgCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCaa                                                                                                     >  4:171984/1‑79 (MQ=255)
                                                                                                                       aGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGAt                            <  6:114706/149‑1 (MQ=255)
                                                                                                                       aGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGAt                            <  7:67868/149‑1 (MQ=255)
                                                                                                                           gggTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCt                        <  7:263517/149‑1 (MQ=255)
                                                                                                                             gTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGt                      <  4:60889/149‑1 (MQ=255)
                                                                                                                              tGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTTGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                     <  6:69499/149‑1 (MQ=255)
                                                                                                                              tGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                     <  5:157235/149‑1 (MQ=255)
                                                                                                                                 aaCCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGATGTGCTGGGTTTTGATGCCTGTCTa                   <  8:133232/148‑1 (MQ=255)
                                                                                                                                  aCCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAAt                                                            >  8:175208/1‑106 (MQ=255)
                                                                                                                                  aCCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAAt                                                            <  7:175208/106‑1 (MQ=255)
                                                                                                                                   cccTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTc                     >  8:184724/1‑144 (MQ=255)
                                                                                                                                        gtggtggCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCg           >  8:248162/1‑149 (MQ=255)
                                                                                                                                        ccggtggCCGCCGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCg           >  8:286581/3‑149 (MQ=255)
                                                                                                                                                ggcggcCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCg    >  2:94683/1‑148 (MQ=255)
                                                                                                                                                 gcggcCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCGGa  <  8:204597/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGATGGAAGCGATCTGACCGCGTTGGGTGAGTCGCCGGCGCATCCTTCTTATGCCTTGGGCATCTTTGGCATGCCGGGCTTCACCGCTTATATGGGCTTGCTCGATATTGGCCGGCCGCAGGCGGGTGAAACCCTGGTGGTGGCGGCGGCCACCGGGCCGGTCGGCGCGACTGTGGGGCAGATCGGCAAGATCAAGGGCTGCCATGTGGTCGGTGTCGCCGGTGGCGCGGAAAAATGCCGGCATGCGGTCGAGGTGCTGGGTTTTGATGCCTGTCTAGATCACCGTGCGCCGGA  >  NC_002947/2822751‑2823044

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: