Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 61,392 A→G Y625H (TAT→CAT)  rapA ← RNA polymerase‑binding ATPase and RNAP recycling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091361,3920AG98.3% 214.5 / ‑3.2 58coding (1873/2907 nt)rapARNA polymerase‑binding ATPase and RNAP recycling factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/27);  minor base . (1/0);  total (31/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 = 8.40e-01

AGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACAT  >  NC_000913/61257‑61527
                                                                                                                                       |                                                                                                                                       
aGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAg                                                                                                                                     <  1:455017/140‑1 (MQ=255)
    tttGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCa                                                                                                                                 <  2:1011863/140‑1 (MQ=255)
                   aaaGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAg                                                                                                                  <  2:876709/140‑1 (MQ=255)
                       ccTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGt                                                                                                              >  1:649530/1‑140 (MQ=255)
                       ccTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGt                                                                                                              >  1:921/1‑140 (MQ=255)
                            ggTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTtct                                                                                                         <  1:149105/140‑1 (MQ=255)
                            ggTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTtct                                                                                                         <  2:296092/140‑1 (MQ=255)
                                 gATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGa                                                                                                    >  1:309559/1‑140 (MQ=255)
                                  gtCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGAt                                                                                                   <  2:994552/139‑1 (MQ=255)
                                  aTCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGAt                                                                                                   >  2:47121/1‑140 (MQ=255)
                                   tCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGTCCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATa                                                                                                  <  2:787163/140‑1 (MQ=255)
                                          aaGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGcaca                                                                                           <  1:29679/140‑1 (MQ=255)
                                            gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATg                                                                                         <  2:474960/140‑1 (MQ=255)
                                            gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATg                                                                                         <  2:223134/140‑1 (MQ=255)
                                            gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATg                                                                                         <  2:748603/140‑1 (MQ=255)
                                            gCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATg                                                                                         <  1:465162/140‑1 (MQ=255)
                                              cAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGa                                                                                       <  2:242438/140‑1 (MQ=255)
                                                gATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATc                                                                                     <  2:8809/140‑1 (MQ=255)
                                                        aTCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATc                                                                             <  1:1002093/140‑1 (MQ=255)
                                                             aTCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTgcg                                                                        >  2:221243/1‑140 (MQ=255)
                                                             aTCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTgcg                                                                        >  1:991391/1‑140 (MQ=255)
                                                               cGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGcc                                                                      <  1:27495/140‑1 (MQ=255)
                                                               cGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGcc                                                                      >  1:578750/1‑140 (MQ=255)
                                                               cGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGcc                                                                      >  1:3585/1‑140 (MQ=255)
                                                                      taCGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGa                                                               <  2:581576/140‑1 (MQ=255)
                                                                        cGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATa                                                             <  1:52680/140‑1 (MQ=255)
                                                                           tATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGa                                                          >  2:212301/1‑140 (MQ=255)
                                                                           tATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGa                                                          >  2:820775/1‑140 (MQ=255)
                                                                            aTCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGAt                                                         >  1:494799/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  2:508912/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  2:656236/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  1:280429/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  2:803902/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  2:874082/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  2:963843/1‑140 (MQ=255)
                                                                              cGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATcc                                                       >  1:827043/1‑140 (MQ=255)
                                                                                 aaaTAGTGCGTCCAGCCTGGCGCCTCTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCGCCAGCACCGACTCTGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAga                                                    >  2:826051/1‑140 (MQ=255)
                                                                                   aTAGTGCGTCCGGCCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACGGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAgacg                                                  >  2:224757/1‑140 (MQ=255)
                                                                                       tGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCa                                              >  1:964889/1‑140 (MQ=255)
                                                                                       tGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCa                                              >  2:691176/1‑140 (MQ=255)
                                                                                       tGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCa                                              >  1:147656/1‑140 (MQ=255)
                                                                                       tGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCa                                              <  2:309559/140‑1 (MQ=255)
                                                                                       tGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCa                                              >  1:584861/1‑140 (MQ=255)
                                                                                       tGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCaa                                             >  1:10113/1‑140 (MQ=255)
                                                                                         cGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCaaa                                            >  1:987222/1‑139 (MQ=255)
                                                                                               gTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACgctg                                      <  1:248347/140‑1 (MQ=255)
                                                                                               gTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACgctg                                      <  1:221243/140‑1 (MQ=255)
                                                                                               gTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACgctg                                      >  2:1017690/1‑140 (MQ=255)
                                                                                                    gCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCa                                 <  1:833393/140‑1 (MQ=255)
                                                                                                    gCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCa                                 <  1:473491/140‑1 (MQ=255)
                                                                                                             cTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCg                        >  1:379541/1‑140 (MQ=255)
                                                                                                             cTCAAATGCATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATACCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCg                        <  2:760412/140‑1 (MQ=255)
                                                                                                                     cATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAAt                >  2:383244/1‑140 (MQ=255)
                                                                                                                     cATCCAGACCTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAAt                >  1:943277/1‑140 (MQ=255)
                                                                                                                            aCCTTCGTGATGCCAGCGCCCCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGt         <  1:224757/140‑1 (MQ=255)
                                                                                                                             ccTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTc        <  2:964889/140‑1 (MQ=255)
                                                                                                                              cTTCGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCa       <  1:963843/140‑1 (MQ=255)
                                                                                                                                 cGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAAc    >  1:229793/1‑140 (MQ=255)
                                                                                                                                 cGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAAc    >  2:701679/1‑140 (MQ=255)
                                                                                                                                 cGTGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAAc    >  2:687308/1‑140 (MQ=255)
                                                                                                                                   tGATGCCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACAt  <  1:691176/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AGTTTTTGATCAGATCGTCAAAGCCTTCGGTTTGATCCGGTGAAGCCAGATAGTTAATCAGATCGTTGTATACGCTATCGTAAATAGTGCGTCCGGTCGGGCAGGTGTGCTCAAATGCATCCAGACCTTCGTGATACCAGCGCACCAGCACCGACTGAGCGGTTTTCTCCAGATAAGGCACATGGATCTGAATATCGTGCGCCTGGCCGATACGATCCAGACGACCAATACGCTGCTCCAGTAGATCCGGGTTGAATGGCAGGTCAAACAT  >  NC_000913/61257‑61527

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: