Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,558,932 2 bp→GC coding (153‑154/831 nt) glpG ← GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,558,9320CG100.0% 268.5 / NA 103V52L (GTG→CTG) glpGGlpG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (57/46);  total (57/46)
*NC_0009133,558,9330GC100.0% 277.5 / NA 103D51E (GAC→GAGglpGGlpG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (57/46);  total (57/46)

Only 100 of 104 total aligned reads displayed.

GTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCACGTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGG  >  NC_000913/3558865‑3559008
                                                                   ||                                                                           
gTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGccc                                                                      >  1:2342002/1‑76 (MQ=255)
gTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGccc                                                                      >  1:3609422/1‑76 (MQ=255)
gTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGccc                                                                      <  1:1954097/76‑1 (MQ=255)
   tGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGtc                                                                   >  2:3402388/1‑75 (MQ=255)
   tGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTa                                                                   >  2:732137/1‑76 (MQ=255)
   tGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTa                                                                   <  1:4453091/76‑1 (MQ=255)
   tGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTa                                                                   >  2:3900536/1‑76 (MQ=255)
    ggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAc                                                                  >  1:2712522/1‑76 (MQ=255)
    ggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAc                                                                  <  1:4439929/76‑1 (MQ=255)
     gcctgcctGCCAGCTCGCCGCCAGATACCGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcg                                                                 >  2:3616903/1‑76 (MQ=255)
     gcctgcctGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcg                                                                 <  1:2694810/76‑1 (MQ=255)
      cctgcctgCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcgc                                                                <  1:1608283/76‑1 (MQ=255)
      cctgcctgCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAcgc                                                                >  1:2699133/1‑76 (MQ=255)
        tgcctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACgcgc                                                              >  1:4513822/1‑76 (MQ=255)
        tgcctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACgcgc                                                              >  1:468722/1‑76 (MQ=255)
        tgcctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGCACgcgc                                                              >  2:2436674/1‑76 (MQ=255)
         ccctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCt                                                             <  1:3417235/75‑1 (MQ=255)
          cctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                            <  1:3173361/76‑1 (MQ=255)
           ctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCg                                                           >  1:4767151/1‑76 (MQ=255)
           ctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCCCg                                                           >  2:2658558/1‑76 (MQ=255)
            tgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGCACGCGCTCgc                                                          >  2:1629023/1‑75 (MQ=255)
            cgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCgg                                                          <  1:1014976/75‑1 (MQ=255)
             gccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGc                                                         <  1:2561992/76‑1 (MQ=255)
              ccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGcc                                                        <  1:3558412/76‑1 (MQ=255)
              ccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGcc                                                        >  2:3619147/1‑76 (MQ=255)
               cAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                       <  1:4693728/76‑1 (MQ=255)
               cAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                       <  1:2629551/76‑1 (MQ=255)
               cAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                       >  2:4677667/1‑76 (MQ=255)
               cAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGCACGCGCTCGGCCt                                                       >  2:1346556/1‑76 (MQ=255)
                  cTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTggg                                                    <  1:3286220/76‑1 (MQ=255)
                  cTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTggg                                                    >  2:2865078/1‑76 (MQ=255)
                   tCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCCCTCGGCCTGGGa                                                   >  2:283708/1‑76 (MQ=255)
                   tCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCCCCGCCCGTACGCGCTCCGCCTggcc                                                   >  2:1099072/1‑74 (MQ=255)
                   tCGCCGCCAGATAACGCGGATCTGCCGGGTCTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTgggc                                                   >  2:2262553/1‑75 (MQ=255)
                   ccgccgcCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCGCGTACGCGCTCGGCCTGGGa                                                   <  2:4719285/75‑1 (MQ=255)
                    cgccgcCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGAc                                                  <  1:3795805/76‑1 (MQ=255)
                     gccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACt                                                 <  2:1223258/76‑1 (MQ=255)
                     cccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACt                                                 <  2:2173834/75‑1 (MQ=255)
                      ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTc                                                <  2:3017998/76‑1 (MQ=255)
                      ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTc                                                <  2:3752698/76‑1 (MQ=255)
                      ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTc                                                >  1:4227103/1‑76 (MQ=255)
                       cgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCa                                               >  2:1559590/1‑76 (MQ=255)
                       cgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCa                                               >  2:3526615/1‑76 (MQ=255)
                        gccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCac                                              >  2:2582195/1‑75 (MQ=255)
                        gccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCAt                                              >  2:3892722/1‑76 (MQ=255)
                         ccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATc                                             >  1:1254416/1‑76 (MQ=255)
                         ccAGAGAGCGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCTCCCGTACGCGCTCGGCCTGGGACTCATc                                             <  2:4610392/76‑1 (MQ=255)
                           aGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCg                                           <  2:3925871/76‑1 (MQ=255)
                           aGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCg                                           <  1:4080594/76‑1 (MQ=255)
                             aTAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgcc                                         >  1:1566485/1‑76 (MQ=255)
                              tAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgcca                                        <  1:4707494/76‑1 (MQ=255)
                               aaCGCGGATCTGCCGGGTTTTCGAGAAACCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccag                                       <  1:582116/76‑1 (MQ=255)
                               aaCGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccag                                       >  1:4140489/1‑76 (MQ=255)
                                aCGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagc                                      <  1:2716169/76‑1 (MQ=255)
                                aCGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagc                                      >  1:985937/1‑76 (MQ=255)
                                 cgcgGATCTGCCGGGTTTTCGAGAAGACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcc                                     <  1:455033/76‑1 (MQ=255)
                                 cgcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcc                                     <  2:4109422/76‑1 (MQ=255)
                                 cgcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcc                                     <  1:1609525/76‑1 (MQ=255)
                                  gcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCTGCCTGGGACTCCTCCgcccgcca                                    >  2:2011139/1‑76 (MQ=255)
                                   cgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGccagccag                                   >  1:2355348/1‑76 (MQ=255)
                                    gggtCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGa                                  <  1:607913/73‑1 (MQ=255)
                                    ggATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGa                                  <  2:4320097/76‑1 (MQ=255)
                                      aTCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGCACTCATCCGCCAGCCAGACa                                >  1:1301201/1‑76 (MQ=255)
                                       tCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCCGCCAGACac                               >  2:2096816/1‑75 (MQ=255)
                                       tCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACAt                               <  2:2355347/76‑1 (MQ=255)
                                          gCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGc                            <  2:3814102/76‑1 (MQ=255)
                                           ccGGGGTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCt                           >  1:2993803/1‑76 (MQ=255)
                                            cGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCtt                          >  2:4384940/1‑76 (MQ=255)
                                            cGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCtt                          >  1:3783816/1‑76 (MQ=255)
                                             gggTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCttt                         >  2:526076/1‑76 (MQ=255)
                                             gggTGTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCttt                         <  2:2044255/76‑1 (MQ=255)
                                             cggTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCttt                         <  1:363466/75‑1 (MQ=255)
                                              ggTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTg                        >  1:12694/1‑76 (MQ=255)
                                              ggTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTg                        >  2:1274829/1‑76 (MQ=255)
                                               gTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTgg                       >  1:4382607/1‑76 (MQ=255)
                                                ttttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGt                      <  1:4648125/76‑1 (MQ=255)
                                                ttttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGt                      >  1:1862033/1‑76 (MQ=255)
                                                ttttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGt                      <  2:4351307/76‑1 (MQ=255)
                                                 tttCGAGAAACCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGtt                     <  1:2865078/76‑1 (MQ=255)
                                                 tttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGtt                     >  1:1480554/1‑76 (MQ=255)
                                                 tttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGtt                     <  1:2968527/76‑1 (MQ=255)
                                                  tgcGAGCAAACGCGCCAGCTCCGCCCGTACGCTCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTa                    <  2:2231022/74‑1 (MQ=255)
                                                   tCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAt                   <  1:3900536/76‑1 (MQ=255)
                                                    cGAGAAAACGCGCCAGCTCCGCCCGTCCGCGCCCGCCCTGGCACTCACCCCCCACCCCCACACCGCTTTGGTCAtt                  >  2:357097/1‑75 (MQ=255)
                                                    cGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtg                  >  2:2872056/1‑76 (MQ=255)
                                                    cGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtg                  <  1:4079354/76‑1 (MQ=255)
                                                      agaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTTGTTAtgtt                >  1:2731628/1‑76 (MQ=255)
                                                      agaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgtt                <  1:3663768/76‑1 (MQ=255)
                                                      agaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgtt                <  2:4355276/76‑1 (MQ=255)
                                                      agaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgtt                >  1:2140781/1‑76 (MQ=255)
                                                       gaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCCGCCAGACATCCCTTTGGTTAtgttg               >  2:6532/1‑76 (MQ=255)
                                                        aaaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttgt              >  2:2576068/1‑76 (MQ=255)
                                                        aaaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGGTAtgttgt              >  1:2685305/1‑76 (MQ=255)
                                                         aaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATgttgtt             >  2:472184/1‑76 (MQ=255)
                                                         aaaCGCGCCAGCCCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATgttgtt             <  2:2173491/76‑1 (MQ=255)
                                                          aaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGttgttg            >  1:2525096/1‑76 (MQ=255)
                                                          aaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGttgttg            >  2:3459944/1‑76 (MQ=255)
                                                          aaCGCGCCAGCTCCCCCAGTACGCGCTCGGCCTGGGCCCCATCCTCCAGCCTGACCTAGATTTGGCTATgttgttc            >  1:3168154/1‑75 (MQ=255)
                                                           aCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGa           >  2:633919/1‑76 (MQ=255)
                                                                    gTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAgg  <  2:3932203/76‑1 (MQ=255)
                                                                   ||                                                                           
GTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCACGTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAATCGTGAGG  >  NC_000913/3558865‑3559008

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