Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,680,491 G→A 7.1% L104L (CTG→TTG)  PP_1475 ← CAIB/BAIF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,680,4910GA 7.1% 163.7 / 2.5 56L104L (CTG→TTG) PP_1475CAIB/BAIF family protein
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base A (2/2);  total (29/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 = 6.73e-01

TAGCGCAGGCCGCCCATCGACTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAACCACCACAGCGAGGTGCC  >  NC_002947/1680343‑1680632
                                                                                                                                                    |                                                                                                                                             
tAGCGCAGGCCGCCCATCGACTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAg                                                                                                                                               >  4:261607/1‑149 (MQ=255)
           gccCATCGACTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGc                                                                                                                                    >  4:242818/1‑149 (MQ=255)
               aTCGACTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGaaa                                                                                                                                <  4:33511/149‑1 (MQ=255)
               aTCGACTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGa                                                                                                                                  >  7:77964/1‑147 (MQ=255)
                    cTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGtttt                                                                                                                           <  1:46501/149‑1 (MQ=255)
                    cTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGtttt                                                                                                                           <  3:183332/149‑1 (MQ=255)
                     tCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTc                                                                                                                          >  4:184694/1‑149 (MQ=255)
                       gCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGa                                                                                                                        >  8:97821/1‑149 (MQ=255)
                       gCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGa                                                                                                                        <  3:184694/149‑1 (MQ=255)
                       gCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGa                                                                                                                        <  7:97821/149‑1 (MQ=255)
                             ggCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGa                                                                                                                  <  3:261607/149‑1 (MQ=255)
                                aCCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGt                                                                                                               >  5:147401/1‑149 (MQ=255)
                                  cGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcg                                                                                                             >  6:74189/1‑149 (MQ=255)
                                  cGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcg                                                                                                             >  4:28597/1‑149 (MQ=255)
                                   gAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcgg                                                                                                            >  2:190195/1‑149 (MQ=255)
                                           gctggTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAAGATGtcggcctcggc                                                                                                     >  5:268617/1‑148 (MQ=255)
                                                  tttCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcg                                                                                                             >  3:172097/1‑133 (MQ=255)
                                                  tttCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcg                                                                                                             <  4:172097/133‑1 (MQ=255)
                                                  tttCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGgc                                                                                             >  4:61228/1‑149 (MQ=255)
                                                  tttCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCACGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCcagcag                                                                                               >  4:31984/1‑147 (MQ=255)
                                                    tCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGgcg                                                                                            <  3:61228/148‑1 (MQ=255)
                                                       tCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcg                                                                                                             >  5:199860/1‑128 (MQ=255)
                                                       tCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGtcg                                                                                                             <  6:199860/128‑1 (MQ=255)
                                                       tCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTg                                                                                        <  5:74189/149‑1 (MQ=255)
                                                        cGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCcagca                                                                                                >  6:165932/1‑140 (MQ=255)
                                                        cGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCcagca                                                                                                <  5:165932/140‑1 (MQ=255)
                                                         ggCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGa                                                                                       <  3:28597/148‑1 (MQ=255)
                                                          gCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGTCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAAGATGTCGGCCTCGGCCAGCAGGCGCTTGAgg                                                                                     <  6:268617/149‑1 (MQ=255)
                                                           cccGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCcagcag                                                                                               >  3:98395/1‑138 (MQ=255)
                                                           cccGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCcagcag                                                                                               <  4:98395/138‑1 (MQ=255)
                                                               gTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATc                                                                                <  3:242818/149‑1 (MQ=255)
                                                                    gCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCgg                                                                            <  8:157684/148‑1 (MQ=255)
                                                                    gCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGc                                                                           >  1:196099/1‑149 (MQ=255)
                                                                          gCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCgg                                                                     <  6:147401/149‑1 (MQ=255)
                                                                            ccGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCgggg                                                                   <  2:196099/149‑1 (MQ=255)
                                                                              gACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTg                                                                 >  4:150808/1‑149 (MQ=255)
                                                                                 aGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCtt                                                              <  1:190195/149‑1 (MQ=255)
                                                                                         cATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCttgaggtt                                                        >  5:202000/1‑147 (MQ=255)
                                                                                                 ggCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGc                                              <  3:150808/149‑1 (MQ=255)
                                                                                                        gTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAACACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCttgaggttga                                                      <  5:155565/134‑1 (MQ=255)
                                                                                                        gTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAACACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCttgaggttga                                                      >  6:155565/1‑134 (MQ=255)
                                                                                                        gTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAACACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCttgaggttga                                                      <  5:159278/134‑1 (MQ=255)
                                                                                                        gTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAACACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCttgaggttga                                                      >  6:159278/1‑134 (MQ=255)
                                                                                                            aGCGCATGCAGCACGTCCCCCCCCAACCCAAGTTTTTCCAGCTCGCCCGGGCGAAAGTTTTCGATCTGGATGTCGGCCTCGCCCAGCAGGCGGGTGAGGATATCGCGGCCCTCGTGGTGCTTGAGGGTGAGCGTGAGCGACTGCttgtt                                   >  4:92538/1‑149 (MQ=255)
                                                                                                                   gcagcaCCTCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGcc                            >  7:85267/1‑149 (MQ=255)
                                                                                                                      gcaCATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTg                          >  3:36876/1‑148 (MQ=255)
                                                                                                                      gcaCATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTg                          >  7:164481/1‑148 (MQ=255)
                                                                                                                               aCCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACTAac                   <  8:119451/146‑1 (MQ=255)
                                                                                                                               aCCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAac                   >  7:119451/1‑146 (MQ=255)
                                                                                                                                 cccAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAaccacct              >  3:8903/1‑148 (MQ=255)
                                                                                                                                   cAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGa                                                                                       <  6:264308/74‑1 (MQ=255)
                                                                                                                                        cAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAaccac                <  3:248130/140‑1 (MQ=255)
                                                                                                                                        cAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAaccac                >  4:248130/1‑140 (MQ=255)
                                                                                                                                        cAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAACCACCACAGCGAg       >  6:80273/1‑149 (MQ=255)
                                                                                                                                        cAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAACCACCACAGCGAg       >  3:178488/1‑149 (MQ=255)
                                                                                                                                              tttCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAACCACCACAGCGAGGTGcc  <  6:202000/148‑1 (MQ=255)
                                                                                                                                              tttCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAACCACCACAGCGAGGTGcc  <  8:85267/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TAGCGCAGGCCGCCCATCGACTCGCCTACGGCACCGAAGCCTGGCTGGTCTTTCATCGGCCCGGTCTGGCCAAAGCCCGACAGGCGCACCATCACCAGGCGCGGGTTCAGCGCATGCAGCACATCCCACCCCAAACCAAGTTTTTCCAGCACGCCCGGGCGAAAGTTTTCGATCAGGATGTCGGCCTCGGCCAGCAGGCGCTTGAGGATATCGCGGCCCTCGGGGTGCTTGAGGTTGAGGGTGAGCGACTGCTTGTTGCGGGCCTGCACGAACCACCACAGCGAGGTGCC  >  NC_002947/1680343‑1680632

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: