Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,077,514 Δ1 bp intergenic (‑165/+149) cmtB ← / ← tktA PTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,077,5130C.100.0% 496.3 / NA 108intergenic (‑164/+150)cmtB/tktAPTS system mannitol‑specific enzyme IIA component/transketolase 1 isozyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (68/40);  total (68/40)

Only 100 of 110 total aligned reads displayed.

GCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCC  >  NC_000913/3077421‑3077599
                                                                                            |                                                                                      
  aaTTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATAcgc                                                                                      >  2:2563212/1‑91 (MQ=255)
       aGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATc                                                                                <  2:1369520/93‑1 (MQ=255)
        gCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCg                                                                               <  1:692708/93‑1 (MQ=255)
         cAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGt                                                                              <  2:138960/93‑1 (MQ=255)
           gggAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGa                                                                            >  2:4396572/1‑93 (MQ=255)
           gggAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGa                                                                            >  2:345901/1‑93 (MQ=255)
            ggAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAt                                                                           <  1:4324765/93‑1 (MQ=255)
            ggAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGAt                                                                           <  1:2706328/93‑1 (MQ=255)
               aaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCa                                                                        <  2:4429868/93‑1 (MQ=255)
               aaCCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCa                                                                        <  1:1341753/93‑1 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                                      >  1:914451/1‑93 (MQ=255)
                 ccGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCata                                                                      <  1:2483414/93‑1 (MQ=255)
                   ggTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAtt                                                                    >  1:1079953/1‑93 (MQ=255)
                   ggTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATAtt                                                                    <  2:638059/93‑1 (MQ=255)
                    gtgtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTa                                                                   <  2:4553890/93‑1 (MQ=255)
                      gtTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACa                                                                 <  2:3479098/93‑1 (MQ=255)
                       ttttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACAt                                                                >  1:288746/1‑93 (MQ=255)
                        tttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATa                                                               <  2:4406360/93‑1 (MQ=255)
                         ttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATaa                                                              <  2:2226790/93‑1 (MQ=255)
                         ttGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATaa                                                              <  1:61685/93‑1 (MQ=255)
                           gAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAtt                                                            <  2:1174932/93‑1 (MQ=255)
                           gAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAAtt                                                            <  1:3097899/93‑1 (MQ=255)
                             aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                                          <  1:2527473/93‑1 (MQ=255)
                             aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                                          <  2:2582334/93‑1 (MQ=255)
                             aaaCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGa                                                          <  1:371687/93‑1 (MQ=255)
                               aCAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGata                                                        <  1:1019009/93‑1 (MQ=255)
                                 aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                                      <  2:2444921/93‑1 (MQ=255)
                                 aGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatt                                                      >  2:3368824/1‑93 (MQ=255)
                                   tGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtattta                                                    <  1:1433188/93‑1 (MQ=255)
                                    gAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttat                                                   >  2:4298318/1‑93 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                  <  1:1629087/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                  <  2:1290932/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                  <  1:3528356/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                  <  1:2380115/93‑1 (MQ=255)
                                     aaTTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGAtatttatt                                                  >  2:4202949/1‑93 (MQ=255)
                                         attaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCat                                              >  2:6882/1‑93 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                             >  2:1846466/1‑93 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                             >  1:4312612/1‑93 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                             >  1:3285158/1‑93 (MQ=255)
                                          ttaATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCata                                             >  2:2352953/1‑93 (MQ=255)
                                           taATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCatat                                            <  2:4486264/93‑1 (MQ=255)
                                            aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTTCATAATTGATATTTATTGCAtata                                           >  1:4602935/1‑93 (MQ=255)
                                            aaTTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCAtata                                           <  1:2695853/93‑1 (MQ=255)
                                              tttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATaaa                                         <  2:1232297/93‑1 (MQ=255)
                                               ttCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAat                                        >  2:3427754/1‑93 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                       >  1:4185055/1‑93 (MQ=255)
                                                tCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAata                                       <  1:4518009/93‑1 (MQ=255)
                                                  ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAtt                                     >  2:4098593/1‑93 (MQ=255)
                                                  ttATATAACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATAtt                                     >  1:13456/1‑93 (MQ=255)
                                                    atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                                   >  1:4435528/1‑93 (MQ=255)
                                                    atataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtg                                   <  1:2581665/93‑1 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                  >  2:3502972/1‑93 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                  >  2:3502971/1‑93 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                  >  1:4271622/1‑93 (MQ=255)
                                                     tataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgt                                  >  2:2140382/1‑93 (MQ=255)
                                                      ataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgtg                                 >  1:1251266/1‑93 (MQ=255)
                                                      ataACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTtgtg                                 <  2:1437402/93‑1 (MQ=255)
                                                       taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                                <  1:1093305/93‑1 (MQ=255)
                                                       taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                                >  2:3538668/1‑93 (MQ=255)
                                                       taACATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGa                                >  2:923385/1‑93 (MQ=255)
                                                        aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                               >  2:1908800/1‑93 (MQ=255)
                                                        aaCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGAt                               >  2:2500500/1‑93 (MQ=255)
                                                         aCATTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATc                              >  2:2526206/1‑93 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                            <  1:1601542/93‑1 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                            >  2:4145326/1‑93 (MQ=255)
                                                           aTTGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTa                            <  2:4009223/93‑1 (MQ=255)
                                                            ttGAGTTATAGATATAACATCAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                           >  2:4581624/1‑93 (MQ=255)
                                                            ttGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                           <  2:964968/93‑1 (MQ=255)
                                                            ttGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                           >  1:238988/1‑93 (MQ=255)
                                                            ttGAGTTATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTac                           >  2:916644/1‑93 (MQ=255)
                                                                 ttATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAacac                      >  2:2110558/1‑93 (MQ=255)
                                                                  tataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACg                     >  1:977766/1‑93 (MQ=255)
                                                                  tataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACg                     >  2:3512123/1‑93 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                   <  1:2050231/93‑1 (MQ=255)
                                                                    taGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGcc                   >  1:3144421/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                 >  2:3598731/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                 >  2:1636342/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                 >  2:87464/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                 >  1:4196777/1‑93 (MQ=255)
                                                                      gATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCtt                 >  1:4221813/1‑93 (MQ=255)
                                                                       atatAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTa                <  1:3708615/93‑1 (MQ=255)
                                                                       atatAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTa                >  2:2982290/1‑93 (MQ=255)
                                                                          taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCt             >  1:3072981/1‑93 (MQ=255)
                                                                          taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCt             >  1:2610855/1‑93 (MQ=255)
                                                                          taACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCt             >  2:1065300/1‑93 (MQ=255)
                                                                             caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTAtt          >  2:2923590/1‑93 (MQ=255)
                                                                              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg         >  2:2077901/1‑93 (MQ=255)
                                                                              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg         <  2:2897996/93‑1 (MQ=255)
                                                                              aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTg         >  1:87709/1‑93 (MQ=255)
                                                                               acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGc        >  2:59781/1‑93 (MQ=255)
                                                                                caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCt       >  1:651475/1‑93 (MQ=255)
                                                                                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt      <  2:532989/93‑1 (MQ=255)
                                                                                 aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCtt      >  1:2472156/1‑93 (MQ=255)
                                                                                  aaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTg     >  2:1977082/1‑93 (MQ=255)
                                                                                   aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGt    >  2:4227640/1‑93 (MQ=255)
                                                                                   aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGt    >  2:350451/1‑93 (MQ=255)
                                                                                    tGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTc   >  1:3962838/1‑93 (MQ=255)
                                                                                    tGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTc   >  1:3621703/1‑93 (MQ=255)
                                                                                     gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTcc  >  1:3820928/1‑93 (MQ=255)
                                                                                     gTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTcc  >  1:255784/1‑93 (MQ=255)
                                                                                            |                                                                                      
GCAATTGAGCAGGGAAACCGGTGTTTTGAAAACAGTGAATTATTAATTTCTTATATAACATTGAGTTATAGATATAACAACAAATGTCAATACCGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCC  >  NC_000913/3077421‑3077599

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