Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,013,717 +G coding (3/999 nt) ysgA ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,013,7171.G99.0% 268.1 / ‑3.5 102coding (3/999 nt)ysgAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (48/53);  total (48/54)
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.72e-01

Only 100 of 107 total aligned reads displayed.

CCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTT  >  NC_000913/4013643‑4013785
                                                                            |                                                                     
ccTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcg                                                                       >  1:3046789/1‑75 (MQ=255)
 cTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcga                                                                      <  1:2159362/76‑3 (MQ=255)
  ttGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcgat                                                                     <  1:1489247/76‑4 (MQ=255)
  ttGTGAAGGAATGGAGGTAAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGcgat                                                                     <  1:2781545/76‑4 (MQ=255)
      gAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcag                                                                 >  1:3989575/1‑76 (MQ=255)
       gggggATGGAGGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                                <  2:1353147/74‑1 (MQ=255)
       aaGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                                >  2:3632667/1‑76 (MQ=255)
       aaGGAATAGAGGTAAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagc                                                                <  1:2764776/76‑1 (MQ=255)
        aGGGATAGAGGGGAACCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagca                                                               <  1:483823/76‑1 (MQ=255)
         ggaaTAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                                              <  1:2552810/73‑1 (MQ=255)
          ggATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAgg                                                             >  1:443969/1‑76 (MQ=255)
            aTAGAGGTGAAGCCCGCCCCAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAgggg                                                           >  2:4033366/1‑74 (MQ=255)
            aTAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                           >  2:4001659/1‑76 (MQ=255)
            aTAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTg                                                           >  2:2633133/1‑76 (MQ=255)
             tAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCCGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGc                                                          >  2:2654500/1‑76 (MQ=255)
               gagGTGAAGCCCGCCACAGTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                        >  2:2588340/1‑76 (MQ=255)
               gagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCCGGGGCaa                                                        >  2:1042748/1‑76 (MQ=255)
               gagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                        >  1:4725738/1‑76 (MQ=255)
               gagGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaa                                                        >  1:4137142/1‑76 (MQ=255)
                agGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaaa                                                       >  2:2969184/1‑76 (MQ=255)
                agGGGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCaaa                                                       <  2:3755929/76‑1 (MQ=255)
                 ggTGACCCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAAt                                                      <  2:1409377/76‑1 (MQ=255)
                   tGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                                    >  1:2990366/1‑76 (MQ=255)
                    gAAGCCCGCCACAATCGCGTCGTCCCGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                                                   <  1:2398624/76‑1 (MQ=255)
                      aGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGa                                                 <  2:2717372/76‑1 (MQ=255)
                        cccGCCACAATCGCGTCGTCCGGGGTGTGAACGGTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                               <  2:4860245/76‑1 (MQ=255)
                        cccGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAtt                                               <  1:2877790/76‑1 (MQ=255)
                          cGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgt                                             >  2:1347370/1‑76 (MQ=255)
                             cacaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                          <  2:3311165/76‑1 (MQ=255)
                             cacaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgt                                          >  1:385997/1‑76 (MQ=255)
                              acaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCCGCCGGCGCAAATCCCGACtgctgtg                                         >  2:485844/1‑76 (MQ=255)
                              acaATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGGGCAAATCCAGATtgttgtg                                         >  2:461207/1‑76 (MQ=255)
                               caATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                                        >  1:4372997/1‑76 (MQ=255)
                               caATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                                        >  1:4471823/1‑76 (MQ=255)
                                aaTCGCGTCGTCCGGGGTCTGAACGCTGGTCGAATCTATATGCGATGCAGCAGGTGCAACTCCAGATtgttgtgtt                                       >  1:1217580/1‑76 (MQ=255)
                                aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                       <  1:2776708/76‑1 (MQ=255)
                                aaTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                                       >  1:1172904/1‑76 (MQ=255)
                                 aTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgttg                                      <  1:4513800/76‑1 (MQ=255)
                                   cgcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                    <  1:4157516/76‑1 (MQ=255)
                                   cgcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                    <  1:2196233/76‑1 (MQ=255)
                                   cgcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                    <  1:3163537/76‑1 (MQ=255)
                                   cgcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGGGCAAATCCAGATTGATGTgttgtt                                    >  1:818008/1‑76 (MQ=255)
                                   cgcgTCGGCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                                    <  2:281673/76‑1 (MQ=255)
                                   cgagtcgccCGGGGGATGAACGATGGTCGAAGCGAGAGGCGCTGCGGCAGGTGCAAATCCAGATTGTTGTgctgtt                                    <  2:4776286/73‑1 (MQ=255)
                                    gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                   <  2:2486374/76‑1 (MQ=255)
                                    gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                   >  2:2325259/1‑76 (MQ=255)
                                    gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                   >  1:1813672/1‑76 (MQ=255)
                                    gcgTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg                                   >  1:3153995/1‑76 (MQ=255)
                                      gtcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGcc                                 <  1:1319226/76‑1 (MQ=255)
                                       tcttcCGGGGTCTGAACGCTGGTCGAAGCGAGAGTCGTTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa                                <  1:1021794/73‑1 (MQ=255)
                                       tcgtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa                                >  2:3696958/1‑76 (MQ=255)
                                        cggccGGGGTCGGAACGCTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCACATCCAGATCGTTGTGTTCTTGCCAt                               <  2:1000853/73‑1 (MQ=255)
                                         gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                              >  1:4392668/1‑76 (MQ=255)
                                         gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                              >  1:2196168/1‑76 (MQ=255)
                                         gtcCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                              >  2:3469051/1‑76 (MQ=255)
                                         gtcCGGGGTCTGAACAATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATg                              >  1:51506/1‑76 (MQ=255)
                                           ccGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                            <  2:2913643/76‑1 (MQ=255)
                                           ccGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                            <  2:2647635/76‑1 (MQ=255)
                                           ccGGGGTCTGAACGATGGGCGATGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt                            <  2:4248262/76‑1 (MQ=255)
                                            cGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTa                           >  2:1558012/1‑76 (MQ=255)
                                             ggggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt                          <  2:3424341/76‑1 (MQ=255)
                                              gggTCTGAACGATGGTCGAAGCGAGAGGAGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAtt                         <  1:4348186/76‑1 (MQ=255)
                                              gggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAtt                         >  1:3940212/1‑76 (MQ=255)
                                               ggTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc                        <  2:773534/76‑1 (MQ=255)
                                                gTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtct                       <  2:3005163/76‑1 (MQ=255)
                                                 tCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc                      <  2:1625659/76‑1 (MQ=255)
                                                  cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                     <  2:733975/76‑1 (MQ=255)
                                                  cTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTcc                     <  2:3485090/76‑1 (MQ=255)
                                                   tGAACGATGGTCGCAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCg                    <  2:2589966/76‑1 (MQ=255)
                                                    gAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGt                   >  1:2092266/1‑76 (MQ=255)
                                                     aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                  <  1:4330995/76‑1 (MQ=255)
                                                     aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                  >  2:3708723/1‑76 (MQ=255)
                                                      aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAc                <  1:2391114/76‑1 (MQ=255)
                                                       tgATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc               >  2:4110009/2‑76 (MQ=255)
                                                        gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc              >  2:1031699/1‑75 (MQ=255)
                                                        gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc              >  2:3458173/1‑75 (MQ=255)
                                                        gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc              <  1:2034820/76‑2 (MQ=255)
                                                         aTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccct             <  1:4398673/76‑3 (MQ=255)
                                                          tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt            >  2:4705074/1‑73 (MQ=255)
                                                          tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt            <  1:2194445/76‑4 (MQ=255)
                                                          tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt            <  2:3858800/76‑4 (MQ=255)
                                                          tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt            >  1:2878326/1‑73 (MQ=255)
                                                             tCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata          <  1:22832/76‑1 (MQ=255)
                                                             tCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata          >  1:3135495/1‑76 (MQ=255)
                                                              cGAAGGGAGAGGCGAGGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa         <  2:261203/76‑1 (MQ=255)
                                                              cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa         <  1:3820872/76‑1 (MQ=255)
                                                              cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa         >  1:2893933/1‑76 (MQ=255)
                                                              cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa         >  2:4453002/1‑76 (MQ=255)
                                                                aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa       >  1:3947962/1‑76 (MQ=255)
                                                                 aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa      <  2:4046991/76‑1 (MQ=255)
                                                                 aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa      <  2:3178773/76‑1 (MQ=255)
                                                                  gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt     <  1:1717032/76‑1 (MQ=255)
                                                                  gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt     <  1:1348947/76‑1 (MQ=255)
                                                                  gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt     <  2:2246510/76‑1 (MQ=255)
                                                                  gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt     <  1:3562185/76‑1 (MQ=255)
                                                                  gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt     >  2:1917826/1‑76 (MQ=255)
                                                                   cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg    >  2:1853417/1‑76 (MQ=255)
                                                                     agagGCGATGCAGCAGGTGCAACTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGtt  <  2:4747963/76‑1 (MQ=255)
                                                                     agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                  <  1:3961743/60‑1 (MQ=255)
                                                                     agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGtt  <  2:2092265/76‑1 (MQ=255)
                                                                            |                                                                     
CCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTT  >  NC_000913/4013643‑4013785

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