Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,930,568 +G intergenic (+66/‑373) kup → / → rbsD low affinity potassium transport system/RbsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,930,5671.G98.6% 389.5 / ‑0.8 142intergenic (+65/‑374)kup/rbsDlow affinity potassium transport system/RbsD
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (66/74);  total (67/75)
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.76e-01

Only 100 of 152 total aligned reads displayed.

ATGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATG  >  NC_000913/3930493‑3930637
                                                                            |                                                                      
atGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTg                                                                       <  2:239942/76‑1 (MQ=255)
atGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTg                                                                       <  1:4164901/76‑1 (MQ=255)
 tGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTgg                                                                      <  2:1787936/76‑2 (MQ=255)
 tGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTgg                                                                      <  2:3215730/76‑2 (MQ=255)
  gTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTggg                                                                     >  1:832091/1‑74 (MQ=255)
  gTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTggc                                                                     <  2:2677662/76‑3 (MQ=255)
   tCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTggcg                                                                    <  1:992584/76‑4 (MQ=255)
   tCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTggcg                                                                    <  2:892719/76‑4 (MQ=255)
      aTAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCgaga                                                                  <  2:4729171/76‑1 (MQ=255)
        agcgtccgtcGGGTACAGATTGATCAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaaa                                                                <  1:43371/74‑1 (MQ=255)
        aaCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaaa                                                                >  1:378452/1‑76 (MQ=255)
        aaCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaaa                                                                <  1:153522/76‑1 (MQ=255)
        aaCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaaa                                                                >  1:5637372/1‑76 (MQ=255)
        aaCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaaa                                                                >  2:3633828/1‑76 (MQ=255)
         aCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAAc                                                               <  1:3599909/76‑1 (MQ=255)
         aCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAAc                                                               <  2:2705550/76‑1 (MQ=255)
         aCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAAc                                                               <  1:837727/76‑1 (MQ=255)
          cgtccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACg                                                              <  2:1985008/76‑1 (MQ=255)
           gtccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGc                                                             >  1:3659408/1‑76 (MQ=255)
           gtccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGc                                                             >  1:1010674/1‑76 (MQ=255)
           gtccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGc                                                             >  2:4693132/1‑76 (MQ=255)
            tccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGCGAAACGcc                                                            >  2:3017550/1‑76 (MQ=255)
               gtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaaa                                                         <  2:1522944/76‑1 (MQ=255)
                tcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAAc                                                        <  1:894855/76‑1 (MQ=255)
                 cGGGTACAGATTGACCAGCTGTCGCGCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCCAACg                                                       <  2:2187533/76‑1 (MQ=255)
                 cGGGTACAGATTGAACAACTGTCGCCTACTTTCTGGCGCGTGGTGGCAAGTTATTGTTGGCGAGAAACGCCAAACg                                                       <  2:2132986/76‑1 (MQ=255)
                  gggTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGt                                                      >  2:2544276/1‑76 (MQ=255)
                   ggTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTa                                                     <  1:5539509/76‑1 (MQ=255)
                      aCAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaa                                                  >  1:4012752/1‑76 (MQ=255)
                       cAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaag                                                 >  2:3865906/1‑76 (MQ=255)
                       cAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaag                                                 >  1:4803740/1‑76 (MQ=255)
                        aGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaaga                                                <  1:742302/76‑1 (MQ=255)
                         gATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACg                                                       <  1:5784802/68‑1 (MQ=255)
                         gATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAgaagaa                                               <  1:3986541/76‑1 (MQ=255)
                         gATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAgaagaa                                               <  2:142408/76‑1 (MQ=255)
                          aTTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGaagaag                                              <  1:1360068/76‑1 (MQ=255)
                          aTTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGaagaag                                              >  1:5242855/1‑76 (MQ=255)
                           ttGCACATCTGTCGCCCACTTTCTGGCGCGTGGTGGCACGTTATGGTTGGCGAGAACCGCCAAACGTAGAAGACGt                                             <  2:4246382/76‑1 (MQ=255)
                           ttGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGt                                             >  2:2788304/1‑76 (MQ=255)
                            tGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCGAAACGTAGAAGAAGtt                                            <  1:3448816/76‑1 (MQ=255)
                            tGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGtt                                            <  1:1740568/76‑1 (MQ=255)
                            tGAACAACTGTCGCCCACGTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCGAAACGTAGAAGAAGtg                                            <  1:3481375/76‑2 (MQ=255)
                             gAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGttt                                           >  1:4833862/1‑76 (MQ=255)
                             gAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGttt                                           >  2:4078573/1‑76 (MQ=255)
                               acaacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTc                                         >  2:5617385/1‑76 (MQ=255)
                                caacTGTCGCGCACTTACTGGCGCGTGGTGGCAAGTTATGTTCGGCTCGAAACGCCAAACTTAGACTAAGTTTTcc                                        <  1:2558662/76‑1 (MQ=255)
                                caacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTcc                                        <  2:2498387/76‑1 (MQ=255)
                                caacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTcc                                        >  2:5006602/1‑76 (MQ=255)
                                caacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTcc                                        >  1:1540886/1‑76 (MQ=255)
                                 aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCa                                       <  2:2274500/76‑1 (MQ=255)
                                  acTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCAc                                      >  1:2106333/1‑76 (MQ=255)
                                  acTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCAc                                      >  2:3098319/1‑76 (MQ=255)
                                  acTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCAc                                      <  2:1860812/76‑1 (MQ=255)
                                    tGTCGCCCACGTTCTGGCTCGGGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTGCCACCg                                    <  2:1723221/76‑1 (MQ=255)
                                     gTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCCCCGc                                   <  1:77293/76‑1 (MQ=255)
                                      tCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGCACACCCCCACAAGAAGAAGAAGATTTCCACCGcc                                  >  2:1587366/1‑75 (MQ=255)
                                      tCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCt                                  <  2:3052516/76‑1 (MQ=255)
                                      tCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCt                                  <  2:5783467/76‑1 (MQ=255)
                                      tCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCt                                  >  2:4054173/1‑76 (MQ=255)
                                         cccACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCg                               >  1:5270563/1‑76 (MQ=255)
                                         cccACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCg                               <  1:2249176/76‑1 (MQ=255)
                                          ccACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCCCCGCTGCgg                              <  1:4398460/76‑1 (MQ=255)
                                           cACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGt                             >  1:2358793/1‑76 (MQ=255)
                                           cACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGt                             >  1:4570275/1‑76 (MQ=255)
                                           cACTGTCTGGCGCGTGGTGGCAAGTTATGGCTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCGGCGGt                             >  1:4927599/1‑76 (MQ=255)
                                            aCTTTCTGGCGCGTGGTGGCACGTTATGGTTGGCGAGAAACGCCAAACATAGAAGAAGTTTTCCACCGCTGCGGTc                            >  2:5227836/1‑76 (MQ=255)
                                            aCTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTc                            >  2:659451/1‑76 (MQ=255)
                                            aCTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTc                            >  1:4201087/1‑76 (MQ=255)
                                            aCTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTc                            <  1:5533171/76‑1 (MQ=255)
                                              tttCTGGCGCGTGGTGGCCAGTTATGGTTGGCGAGGCCCGCTGATCGTAGAAGAAGTTTTCCTCGGCTGCGGTCtt                          <  2:5257011/76‑2 (MQ=255)
                                              tttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTg                          <  2:4892785/76‑1 (MQ=255)
                                              tttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTg                          >  1:5366283/1‑76 (MQ=255)
                                              tttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTg                          <  2:1675859/76‑1 (MQ=255)
                                               ttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTgg                         >  1:2927705/1‑76 (MQ=255)
                                                 cgggCGGGGGGTGGCAGGTTATGGTGGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGaa                       <  2:5709534/74‑1 (MQ=255)
                                                 cTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGaa                       <  2:5576600/76‑1 (MQ=255)
                                                  tGGCGCGTGGTGGCAAGTTATGGTTGGCTAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAg                      >  1:2943052/1‑76 (MQ=255)
                                                    gcgcGTGGTGGCAGGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGa                    <  1:2177454/76‑1 (MQ=255)
                                                    gcgcGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGa                    <  2:2439815/76‑1 (MQ=255)
                                                    gcgcGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGa                    >  2:3865678/1‑76 (MQ=255)
                                                     cgcgTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGAt                   >  1:585337/1‑76 (MQ=255)
                                                      gcgTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGAtt                  >  1:2734062/1‑76 (MQ=255)
                                                      gcgTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGAtt                  >  1:5351490/1‑76 (MQ=255)
                                                       cgTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTa                 <  2:1710859/76‑1 (MQ=255)
                                                        gtggtgGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTaa                >  1:4991335/1‑76 (MQ=255)
                                                         tggtggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAg               >  2:5651003/1‑76 (MQ=255)
                                                          ggtggCAAGTTATGGTTGGCGAGACACGCCCAACGTAGACGAAGTTTTCCACCGCTTCGTTCTGGCATGATTacgc              <  2:43249/76‑4 (MQ=255)
                                                           gtggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGtt             >  1:1244872/1‑76 (MQ=255)
                                                             ggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGc           <  2:1267558/76‑1 (MQ=255)
                                                             ggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGc           >  1:5691952/1‑76 (MQ=255)
                                                                aaGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCgg        >  2:2353719/1‑76 (MQ=255)
                                                                aaGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCgg        >  2:5605763/1‑76 (MQ=255)
                                                                aaGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCgg        >  1:1804975/1‑76 (MQ=255)
                                                                aaGTTATGGTTGGCGAGAAACACCAAACGTAGAAGAAGTTTTCCACCGCTGCGCTCTGGAAGGATCAAGTTGCCgg        >  2:4646047/1‑76 (MQ=255)
                                                                 aGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGga       <  2:2815547/76‑1 (MQ=255)
                                                                   ttATGGTTGGCGAGAAACCCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgatg     >  1:3349652/1‑76 (MQ=255)
                                                                    tATGGTTGGCGAGAAACGCCCAACGTAGAGGAAGCTCTCCACCGCTGCGGTGCGGCAGGCATAAGCTGGCGgatga    >  2:607472/1‑76 (MQ=255)
                                                                    tATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgatga    >  1:1100513/1‑76 (MQ=255)
                                                                      tGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCCCCGCTGCGGTCTGGAAGGATTAAGTTGCCGGatgatg  <  2:1370013/76‑1 (MQ=255)
                                                                      tGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGGatgatg  <  2:4991334/76‑1 (MQ=255)
                                                                            |                                                                      
ATGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATG  >  NC_000913/3930493‑3930637

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