Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,016,722 T→G 9.7% intergenic (‑38/+189) PP_5728 ← / ← PP_5267 rhs repeat‑associated core domain‑containing protein/cytochrome c5

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,016,7220TG 9.7% 154.6 / 4.6 51intergenic (‑38/+189)PP_5728/PP_5267rhs repeat‑associated core domain‑containing protein/cytochrome c5
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base G (3/2);  total (26/25)
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

ATAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAA  >  NC_002947/6016574‑6016869
                                                                                                                                                    |                                                                                                                                                   
aTAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGttctt                                                                                                                                                     >  3:165601/1‑149 (MQ=255)
  aGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGttctt                                                                                                                                                     >  7:27092/1‑147 (MQ=255)
  aGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGACAGTcttcttaa                                                                                                                                                   >  7:188138/1‑147 (MQ=255)
         aGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTATCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCt                                                                                                                                            <  5:3637/149‑1 (MQ=255)
          ggCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCt                                                                                                                                            >  2:21305/1‑148 (MQ=255)
               aTAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGTAGCGTCAGTGTTCTTCAGAACCCTGaaaaa                                                                                                                                      <  1:24059/149‑1 (MQ=255)
                   gCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATTGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGa                                                                                                                                  <  2:52689/149‑1 (MQ=255)
                    ccAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGCGTTCTTCAGAACCCTGAAAAATAGaa                                                                                                                                 <  5:243952/149‑1 (MQ=255)
                         tGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGt                                                                                                                            <  8:27092/149‑1 (MQ=255)
                          gCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGt                                                                                                                            >  4:228229/1‑148 (MQ=255)
                            cGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGacgctt                                                                                                                                     >  6:189471/1‑132 (MQ=255)
                            cGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGacgctt                                                                                                                                     <  5:189471/137‑6 (MQ=255)
                                           tggATCAGCTTTTGCGTCTGTGTTCGTCCGGCTGCATGTTTCTCGGGACAGTTATCGCCCAGTTCCATTGTTATCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACACTGTAAAATAGATCTTTTAAAAATGGTATTTCTCAt                                                                                                          <  8:188138/148‑1 (MQ=255)
                                             gATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATgg                                                                                                        <  5:75034/149‑1 (MQ=255)
                                               tCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGGCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGaa                                                                                                      <  3:45221/149‑1 (MQ=255)
                                                  gcgtttGCTGCTGTGTTCGACCTTCTTCATGTTGCTCGGGACAGTTATCGCACAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATATAACTTGTAAAAATGCTAGTTCTCATGGAATCa                                                                                                   <  4:162413/146‑1 (MQ=255)
                                                          ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                           >  7:225188/1‑149 (MQ=255)
                                                          ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                           <  3:228229/149‑1 (MQ=255)
                                                          ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                           >  4:196728/1‑149 (MQ=255)
                                                           gCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGac                                                                                                                                         <  3:196975/102‑2 (MQ=255)
                                                           gCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGac                                                                                                                                         >  4:196975/1‑101 (MQ=255)
                                                                 ttCGACCGACTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTg                                                                                    >  4:179154/1‑149 (MQ=255)
                                                                         gCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCt                                                                            >  1:204545/1‑149 (MQ=255)
                                                                          cTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTg                                                                           <  5:246161/149‑1 (MQ=255)
                                                                            gCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACGCTGAAAAATAGaa                                                                                                                                 <  6:24918/93‑1 (MQ=255)
                                                                            gCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGaa                                                                                                                                 >  5:24918/1‑93 (MQ=255)
                                                                            gCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTg                                                                           >  8:134534/1‑147 (MQ=255)
                                                                             cATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGga                                                                         <  6:246058/148‑1 (MQ=255)
                                                                              aTGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGgag                                                                        >  4:42715/1‑148 (MQ=255)
                                                                                gTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTCGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACa                                                                     >  6:227720/1‑149 (MQ=255)
                                                                                 ttGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACa                                                                     <  4:165601/148‑1 (MQ=255)
                                                                                 ttGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTCGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACAt                                                                    <  5:227720/149‑1 (MQ=255)
                                                                                     tCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                >  4:109892/1‑149 (MQ=255)
                                                                                       gggACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                <  3:237308/147‑1 (MQ=255)
                                                                                         gACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGa                                                             >  6:58728/1‑148 (MQ=255)
                                                                                              ttATCGCCCAGATACATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCAATGGACATgg                                                        >  2:71000/1‑148 (MQ=255)
                                                                                                aTCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGctc                                                     >  5:172747/1‑149 (MQ=255)
                                                                                                   gCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTa                                                   <  8:225188/148‑1 (MQ=255)
                                                                                                          tCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGggttg                                           <  3:42715/149‑1 (MQ=255)
                                                                                                           ccATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGggttgg                                          >  4:68892/1‑149 (MQ=255)
                                                                                                                      caccacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCt                               >  3:146018/1‑149 (MQ=255)
                                                                                                                       accacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTg                              >  2:185440/1‑149 (MQ=255)
                                                                                                                       accacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAACGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTg                              >  7:28203/1‑149 (MQ=255)
                                                                                                                         cacTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGa                            <  3:196728/149‑1 (MQ=255)
                                                                                                                            ttCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCt                         >  2:250523/1‑149 (MQ=255)
                                                                                                                            ttCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCt                         >  3:111696/1‑149 (MQ=255)
                                                                                                                              cAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                >  5:53863/1‑108 (MQ=255)
                                                                                                                              cAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCAt                                                                <  6:53863/108‑1 (MQ=255)
                                                                                                                                    gAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCggg                 >  3:5333/1‑149 (MQ=255)
                                                                                                                                      aGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGca               <  3:68892/149‑1 (MQ=255)
                                                                                                                                         ttCAGGGTTCTGAAGAACACTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCAATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGcacca            >  6:161463/2‑149 (MQ=255)
                                                                                                                                          tCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGcaccac           <  4:111696/149‑1 (MQ=255)
                                                                                                                                            aGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGca               <  4:201114/143‑1 (MQ=255)
                                                                                                                                            aGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGca               >  3:201114/1‑143 (MQ=255)
                                                                                                                                                  cttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCtgtg                                 <  7:242612/119‑1 (MQ=255)
                                                                                                                                                  cttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCtgtg                                 >  8:242612/1‑119 (MQ=255)
                                                                                                                                                   ttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACaaa  <  2:204545/149‑1 (MQ=255)
                                                                                                                                                   ttcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACa    >  4:214611/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATAGCCATAAGGCGAATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAA  >  NC_002947/6016574‑6016869

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: