Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 410,479 A→G S112S (TCA→TCG mak → manno(fructo)kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913410,4790AG100.0% 294.4 / NA 80S112S (TCA→TCGmakmanno(fructo)kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/40);  total (40/40)

CAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCG  >  NC_000913/410382‑410575
                                                                                                 |                                                                                                
cAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGa                                                                                                <  2:1857563/100‑1 (MQ=255)
cAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGa                                                                                                <  1:1777556/100‑1 (MQ=255)
  aCCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAg                                                                                              <  1:268971/100‑1 (MQ=255)
  aCCTGGCTCAACGGTCAGCCAGGCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGACAAATGACGCTAACTGTCTGGCGGTTTCGGAAg                                                                                              <  2:2783155/100‑1 (MQ=255)
   ccTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGc                                                                                             >  2:505011/1‑100 (MQ=255)
    cTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCa                                                                                            <  2:2027071/100‑1 (MQ=255)
      ggCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGt                                                                                          <  1:1915275/100‑1 (MQ=255)
       gCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGTCAAATGACGCTAACTGTCTGGCGGCTTCGGAAGCAGTa                                                                                         <  2:1840031/100‑1 (MQ=255)
        ntCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAg                                                                                        >  2:2219392/2‑100 (MQ=255)
         tCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGa                                                                                       <  1:513643/100‑1 (MQ=255)
          cAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGAt                                                                                      <  2:879017/100‑1 (MQ=255)
            aCGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATgg                                                                                    <  1:1765899/100‑1 (MQ=255)
            aCGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATgg                                                                                    <  1:2127734/100‑1 (MQ=255)
             cGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGgc                                                                                   <  2:2253078/100‑1 (MQ=255)
                tCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGcgcg                                                                                <  2:137087/100‑1 (MQ=255)
                tCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGcgcg                                                                                <  2:2855526/100‑1 (MQ=255)
                tCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGcgcg                                                                                >  2:2944926/1‑100 (MQ=255)
                 cAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCgg                                                                               >  1:1488717/1‑100 (MQ=255)
                  aGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGc                                                                              <  1:517329/100‑1 (MQ=255)
                  aGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGc                                                                              >  1:2309719/1‑100 (MQ=255)
                   gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGc                                                                              >  2:2313170/1‑99 (MQ=255)
                   gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCa                                                                             >  2:421862/1‑100 (MQ=255)
                   gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCa                                                                             >  2:1166154/1‑100 (MQ=255)
                   gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCa                                                                             >  1:2845174/1‑100 (MQ=255)
                   gCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCa                                                                             >  2:846505/1‑100 (MQ=255)
                    ccATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAg                                                                            >  1:569356/1‑100 (MQ=255)
                     cATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGc                                                                           >  2:467730/1‑100 (MQ=255)
                      aTTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCg                                                                          >  2:1315949/1‑100 (MQ=255)
                         cGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGa                                                                       >  2:1896426/1‑100 (MQ=255)
                          gATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAg                                                                      >  1:1403049/1‑100 (MQ=255)
                          gATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAg                                                                      <  1:846508/100‑1 (MQ=255)
                           aTAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgc                                                                     <  1:753890/100‑1 (MQ=255)
                            tAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGGGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcg                                                                    >  2:2010092/1‑100 (MQ=255)
                            tAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcg                                                                    >  1:1199432/1‑100 (MQ=255)
                            ctaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcg                                                                    <  1:1888132/98‑1 (MQ=255)
                             aaaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCgccatcta                                                                   >  2:2568563/1‑93 (MQ=255)
                              aaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAgcgag                                                                  >  2:1888121/1‑98 (MQ=255)
                              aaGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCa                                                                  <  1:171374/100‑1 (MQ=255)
                                  cTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACg                                                              >  1:2745021/1‑100 (MQ=255)
                                    tAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGt                                                            >  1:696809/1‑100 (MQ=255)
                                     aaGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTa                                                           <  1:152875/100‑1 (MQ=255)
                                     aaGCGCGAGGTTGAAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTa                                                           <  1:1181648/100‑1 (MQ=255)
                                          cgAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGc                                                      >  1:960602/1‑100 (MQ=255)
                                             ggTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGt                                                   <  2:175221/100‑1 (MQ=255)
                                              gTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTg                                                  >  1:319721/1‑100 (MQ=255)
                                               ttGCAGCGGGAAGTGCGGCTGGCAAATGGCGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGa                                                 <  1:1881890/100‑1 (MQ=255)
                                               ttGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGa                                                 <  2:1199426/100‑1 (MQ=255)
                                                 gCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTGTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGAt                                                <  1:415396/99‑1 (MQ=255)
                                                  cAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTa                                              >  2:2254939/1‑100 (MQ=255)
                                                    gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                            <  1:870945/100‑1 (MQ=255)
                                                    gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                            >  1:472746/1‑100 (MQ=255)
                                                    gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                            <  2:2237928/100‑1 (MQ=255)
                                                    gCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATc                                            >  1:1648911/1‑100 (MQ=255)
                                                          aGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACg                                      <  2:2525483/100‑1 (MQ=255)
                                                           gTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACgg                                     <  1:2069577/100‑1 (MQ=255)
                                                           gTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACgg                                     <  1:2326972/100‑1 (MQ=255)
                                                              cGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGAt                                  <  1:2219404/100‑1 (MQ=255)
                                                                  tGGCAAATGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATgcgg                              <  1:1896437/100‑1 (MQ=255)
                                                                        aTGACGCTAACTGTCTTGCGGTTTCGGAAGCAGTGGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGGTGTGATTATCGGCACGGGATgcggcg                            >  2:1151479/1‑96 (MQ=255)
                                                                         tGACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGTGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGc                       >  2:2724858/1‑100 (MQ=255)
                                                                          gACGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCg                      <  1:2407339/100‑1 (MQ=255)
                                                                          gACGCGAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCg                      <  1:2313184/100‑1 (MQ=255)
                                                                            cGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTg                    <  1:2402939/100‑1 (MQ=255)
                                                                            cGCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTg                    >  1:2136695/1‑100 (MQ=255)
                                                                             gCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTgg                   >  1:506621/1‑100 (MQ=255)
                                                                             gCTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTgg                   >  1:2818789/1‑100 (MQ=255)
                                                                              cTAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGc                  <  2:2507722/100‑1 (MQ=255)
                                                                               tAACTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCa                 >  1:2799090/1‑100 (MQ=255)
                                                                                  cggTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTc              >  1:2007800/3‑100 (MQ=255)
                                                                                  cTGTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTc              >  1:1445672/1‑100 (MQ=255)
                                                                                    gTCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCaa            >  1:1656262/1‑100 (MQ=255)
                                                                                     tCTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAAt           <  1:2441349/100‑1 (MQ=255)
                                                                                      cTGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGGAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATg          <  2:2099451/100‑1 (MQ=255)
                                                                                       tGGCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgg         >  1:390526/1‑100 (MQ=255)
                                                                                        ggCGGTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGGGTGGCATTCAATggg        >  1:2946377/1‑100 (MQ=255)
                                                                                           ggTTTCGCAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATg          <  2:2439842/95‑1 (MQ=255)
                                                                                            gTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcggg    <  2:2695072/100‑1 (MQ=255)
                                                                                            gTTTCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcggg    >  1:2128891/1‑100 (MQ=255)
                                                                                             tttCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATgggcgggc   >  1:2535951/1‑100 (MQ=255)
                                                                                              ttCGGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGggcgggcg  >  1:13086/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAACCTGGCTCAACGGTCAGCCATTCGATAAAGACTTAAGCGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCG  >  NC_000913/410382‑410575

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