Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,945,351 +C coding (1/600 nt)
coding (240/240 nt)
yifDA ←
yifDA ←
transcription regulator
transcription regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,945,3501.C99.2% 315.8 / ‑4.5 120E81G (GAG→GGG) yifDAtranscription regulator
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (47/72);  minor base A (0/1);  total (47/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 123 total aligned reads displayed.

ATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTT  >  NC_000913/3945277‑3945420
                                                                          |                                                                      
 ttAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACctcct                                                                     <  1:128611/76‑3 (MQ=255)
  tAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACctcctt                                                                    >  1:3330044/1‑73 (MQ=255)
     gCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACg                                                                  <  2:4756503/76‑1 (MQ=255)
     gCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACg                                                                  >  2:343863/1‑76 (MQ=255)
     gCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACg                                                                  <  2:4108237/76‑1 (MQ=255)
      cATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACgg                                                                 <  2:4240091/76‑1 (MQ=255)
      cATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACgg                                                                 <  1:3486927/76‑1 (MQ=255)
      cATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACgg                                                                 >  1:3090610/1‑76 (MQ=255)
       aTACATTCTCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACggg                                                                <  1:973640/76‑1 (MQ=255)
       aTACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACggg                                                                >  1:1945967/1‑76 (MQ=255)
        tACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggc                                                               <  1:3885798/76‑1 (MQ=255)
        aaCATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggc                                                               <  2:83120/75‑1 (MQ=255)
         aCATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggcg                                                              <  1:2746280/76‑1 (MQ=255)
           aTTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGggcggc                                                            >  2:2466916/1‑76 (MQ=255)
             tCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                          <  2:2439300/76‑1 (MQ=255)
             tCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                          >  1:3478471/1‑76 (MQ=255)
             tCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCt                                                          <  2:2462296/76‑1 (MQ=255)
              cccACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTg                                                         <  1:2405572/76‑1 (MQ=255)
                cacaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGcc                                                       <  1:3906549/76‑1 (MQ=255)
                cacaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGcc                                                       >  2:2277070/1‑76 (MQ=255)
                 acaACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCa                                                      <  1:3136293/76‑1 (MQ=255)
                  caACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAg                                                     >  1:2643651/1‑76 (MQ=255)
                    aCGAGGCGCTGGCCCCGATAGAAAACTCGTTATGTCGTGAGGTAGGCGCCACCTCCTTACGGGCGGCCTGCCAGGt                                                   <  2:1652705/76‑1 (MQ=255)
                    aCGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAggg                                                   >  2:3914075/1‑75 (MQ=255)
                     cGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTg                                                  >  2:3667979/1‑76 (MQ=255)
                     cGAGGCGCTCGCACCGATAGAAAACTCGTTATGGCGTGAGTGATGCGCCACCTACTTACGGGCGGCCTGCCAGGTg                                                  <  2:2317128/76‑1 (MQ=255)
                     cGAGGCGCTAGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTg                                                  <  2:2313224/76‑1 (MQ=255)
                      gAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCCGgcgc                                                 >  2:3281424/1‑76 (MQ=255)
                      gAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGc                                                 <  2:2602830/76‑1 (MQ=255)
                      gAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGc                                                 >  1:2965601/1‑76 (MQ=255)
                        ggCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTc                                               <  2:672135/76‑1 (MQ=255)
                         gcgcTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCa                                              >  2:2115510/1‑76 (MQ=255)
                         gcgcTGGCACCGATAGAAAACTAGTTATGTCGTGAGGTATGCGCAACCTCCTTACGGGCGACGTGCCAGGTGCtat                                              >  2:3073681/1‑74 (MQ=255)
                          cgcTGGCCCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCAt                                             <  2:1370191/76‑1 (MQ=255)
                          cgcTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCAt                                             <  1:1614144/76‑1 (MQ=255)
                            cTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGa                                           <  2:1792866/76‑1 (MQ=255)
                            cTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGa                                           <  1:2154090/76‑1 (MQ=255)
                              ggCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGc                                         <  1:3549097/76‑1 (MQ=255)
                               gCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCg                                        >  2:1130048/1‑76 (MQ=255)
                                cACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGt                                       >  1:4616040/1‑76 (MQ=255)
                                 aCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCCGGTGCTCATGAGCGtt                                      >  2:4449676/1‑76 (MQ=255)
                                 aCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGtt                                      >  1:159080/1‑76 (MQ=255)
                                 aCCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGtt                                      >  1:2053126/1‑76 (MQ=255)
                                  ccGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGttt                                     <  2:1390793/76‑1 (MQ=255)
                                    gATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCt                                   <  1:3215732/76‑1 (MQ=255)
                                    gATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCt                                   >  2:4794571/1‑76 (MQ=255)
                                    caTAGAAAACTCGTTATGTCGTGAGGTATGCGCCCCCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCt                                   <  1:4590102/75‑1 (MQ=255)
                                     aTAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGATCGTTTCTg                                  >  1:2730615/1‑76 (MQ=255)
                                      tAGAGAGCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGc                                 <  2:2928869/76‑1 (MQ=255)
                                      tAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTCCGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGc                                 <  1:3851575/76‑1 (MQ=255)
                                       aGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGCGCGGCCTGCCAGCTGCTCATGATCGTTTCtcct                                >  2:1144968/1‑73 (MQ=255)
                                         aaaaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATa                              <  1:690858/76‑1 (MQ=255)
                                         aaaaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATa                              <  1:2346609/76‑1 (MQ=255)
                                         aaaaCTAGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATa                              <  1:2365144/76‑1 (MQ=255)
                                          aaaCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATaa                             >  2:3569405/1‑76 (MQ=255)
                                            aCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAgg                           <  2:2533601/76‑1 (MQ=255)
                                            aCTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAgg                           <  2:4408974/76‑1 (MQ=255)
                                             cTCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAAGTGCTCATGAGCGTTTCTGCATAAGGc                          <  1:1371667/76‑1 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGCGCGTTTCTGCATACGGCa                         >  2:3148817/1‑76 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                         >  1:3974605/1‑76 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                         <  1:1924707/76‑1 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                         <  1:2059529/76‑1 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                         >  2:3122860/1‑76 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                         >  2:3462935/1‑76 (MQ=255)
                                              tCGTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCa                         <  1:621923/76‑1 (MQ=255)
                                                gtgaTGTCGTGAGGTATGCGCCCCCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagt                       <  2:181842/73‑1 (MQ=255)
                                                gTTATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagt                       <  1:4599041/76‑1 (MQ=255)
                                                  tATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtag                     >  1:1251143/1‑76 (MQ=255)
                                                  tATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtag                     <  2:3784595/76‑1 (MQ=255)
                                                  tATGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtag                     >  2:2432076/1‑76 (MQ=255)
                                                   aTGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCagtagt                    >  1:3312956/1‑76 (MQ=255)
                                                    tGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATACGGCAGTAGtt                   >  2:4599201/1‑76 (MQ=255)
                                                    tGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtt                   <  2:3312956/76‑1 (MQ=255)
                                                    tGTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtt                   <  2:3814861/76‑1 (MQ=255)
                                                     gTCGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGttt                  <  1:3169931/76‑1 (MQ=255)
                                                     gTCGGGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGttt                  <  2:268599/76‑1 (MQ=255)
                                                      tCGTGAGGTATGCCCCCCCTCCTGACGGGCGGCCTGCTAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGtttt                 <  1:3073681/76‑1 (MQ=255)
                                                       cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGttttt                <  2:2826120/76‑1 (MQ=255)
                                                       cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGttttt                >  2:709234/1‑76 (MQ=255)
                                                       cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGttttt                <  1:3462935/76‑1 (MQ=255)
                                                       cGTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCGCATGAGCGTTTCTGAATAAGGCAGGAGttttt                >  2:2838019/1‑76 (MQ=255)
                                                        gTGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTc               <  1:4683493/76‑1 (MQ=255)
                                                         tGAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCa              <  2:2741671/76‑1 (MQ=255)
                                                          gAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCAc             >  1:488411/1‑76 (MQ=255)
                                                          gAGGTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCAc             <  2:1626888/76‑1 (MQ=255)
                                                            ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCg           <  2:3877678/76‑1 (MQ=255)
                                                            ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCg           <  1:3293037/76‑1 (MQ=255)
                                                            ggTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCg           >  1:3357529/1‑76 (MQ=255)
                                                             gTATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCgg          >  1:3990937/1‑76 (MQ=255)
                                                              tATGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGc         <  1:2670243/76‑1 (MQ=255)
                                                               aTGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCa        <  2:2421778/76‑1 (MQ=255)
                                                                tGCGCCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAg       >  2:600804/1‑76 (MQ=255)
                                                                 gcgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGc      <  1:1106460/76‑1 (MQ=255)
                                                                 gcgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGc      <  1:2432226/76‑1 (MQ=255)
                                                                  cgcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCg     >  2:1337632/1‑76 (MQ=255)
                                                                   gcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCCTAAGGCAGTAGTTTTTCACCGGCAGCgg    >  2:470717/1‑76 (MQ=255)
                                                                   gcCACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCgg    >  1:1517204/1‑76 (MQ=255)
                                                                     cACCTCCTTACGGGCGTCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGtt  <  1:2970112/76‑1 (MQ=255)
                                                                     cACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGtt  >  2:4512262/1‑76 (MQ=255)
                                                                     cACCTCCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGtt  <  2:2520701/76‑1 (MQ=255)
                                                                          |                                                                      
ATTAAGCATACATTCCCACAACGAGGCGCTGGCACCGATAGAAAACTCGTTATGTCGTGAGGTATGCGCCACCTCTTACGGGCGGCCTGCCAGGTGCTCATGAGCGTTTCTGCATAAGGCAGTAGTTTTTCACCGGCAGCGGTT  >  NC_000913/3945277‑3945420

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