Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 861,642 G→A 100% A809T (GCC→ACC)  RCAP_rcc00783 → diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034861,6420G→A100.0% 649.7 183A809T (GCC→ACC) RCAP_rcc00783diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (89/92):  ref base (0/0):  total (89/94)

Only 100 of 198 total aligned reads displayed.

TGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACGC  >  NC_014034/861544‑861739
                                                                                                  |                                                                                                 
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                 >  1:3143416/1‑101 (MQ=255)
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                 >  2:574022/1‑101 (MQ=255)
   ccGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGc                                                                                              <  2:790045/101‑1 (MQ=255)
   ccGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGc                                                                                              <  2:4078626/101‑1 (MQ=255)
    cGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCt                                                                                             >  2:1725609/1‑101 (MQ=255)
      aGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTAt                                                                                           <  1:861216/101‑1 (MQ=255)
      aGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTAt                                                                                           >  1:3010331/1‑101 (MQ=255)
       ggATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTAtt                                                                                          >  1:2694462/1‑101 (MQ=255)
        gATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTc                                                                                         >  2:3272992/1‑101 (MQ=255)
         aTCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCGGGCGATCGACGATTTCGGCACCGGCTATTcc                                                                                        <  1:2550059/101‑1 (MQ=255)
          tCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCt                                                                                       <  2:2081723/101‑1 (MQ=255)
            aGCTGGCCGACCGGGTGCCGGGCGCGGTCGAGGCGCTTGGGCGGGTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGt                                                                                     <  2:2342044/101‑1 (MQ=255)
              cTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCt                                                                                   >  2:2793233/1‑101 (MQ=255)
               tGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTg                                                                                  <  2:3670150/101‑1 (MQ=255)
               tGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTg                                                                                  >  2:3042894/1‑101 (MQ=255)
                   cGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCt                                                                              <  1:3276017/101‑1 (MQ=255)
                   cGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCt                                                                              >  2:1782261/1‑101 (MQ=255)
                    gcccGGGTTCCGGGCGCGGGCGAGGCGCTTGGGGGGGTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                             <  2:2746346/99‑1 (MQ=38)
                    gACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCtt                                                                             >  1:3044201/1‑100 (MQ=255)
                      ccGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATc                                                                           <  2:2225931/101‑1 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGTGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                         >  2:3995295/1‑101 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                         >  2:120992/1‑101 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                         >  1:2401747/1‑101 (MQ=255)
                        gggTGCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                         <  1:1004602/101‑1 (MQ=255)
                         ggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTgg                                                                        >  1:1124685/1‑101 (MQ=255)
                         ggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATCTCGGCACCGGCTATTCCTGTCTGGCCTATCTgg                                                                        >  1:3862541/1‑101 (MQ=255)
                              cGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCgg                                                                   >  2:146297/1‑101 (MQ=255)
                                ggCGGGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCt                                                                 <  1:4172703/101‑1 (MQ=255)
                                ggCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCt                                                                 >  1:3427352/1‑101 (MQ=255)
                                ggCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCt                                                                 <  2:1613774/101‑1 (MQ=255)
                                 gcgcGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTg                                                                <  1:574090/101‑1 (MQ=255)
                                   gcgGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGcc                                                              <  1:2470478/101‑1 (MQ=255)
                                     ggTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCgg                                                            >  1:23784/1‑101 (MQ=255)
                                      gTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTggccgt                                                           >  1:4021447/1‑96 (MQ=255)
                                      gTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGt                                                           >  1:3594019/1‑101 (MQ=255)
                                      gTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGt                                                           >  1:3029562/1‑101 (MQ=255)
                                        cGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtct                                                         >  1:1690254/1‑101 (MQ=255)
                                         gAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTTTTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtctc                                                        <  2:3141789/101‑1 (MQ=255)
                                         gAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtctc                                                        <  2:2944665/101‑1 (MQ=255)
                                         gAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtctc                                                        <  2:3534144/101‑1 (MQ=255)
                                         gAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtctc                                                        <  1:3114729/101‑1 (MQ=255)
                                            gcgcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGAt                                                     <  2:2196848/101‑1 (MQ=255)
                                            gcgcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGAt                                                     <  2:697464/101‑1 (MQ=255)
                                                  gggCGGTTGCGCGCGCGCGGGCTGCGGCTGGGGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGaa                                               >  2:191432/1‑101 (MQ=255)
                                                   ggCGGTTGCGCGCGCGCGGGGTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAg                                              <  2:899323/101‑1 (MQ=255)
                                                   ggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAg                                              <  1:1221806/101‑1 (MQ=255)
                                                   ggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAg                                              <  2:2445536/101‑1 (MQ=255)
                                                   ggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCGGGCGCGGCTGCCGGTCTCGATGCTGAAg                                              <  2:2561669/101‑1 (MQ=255)
                                                    gCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGATGa                                             >  2:3748152/1‑101 (MQ=255)
                                                    gCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGa                                             <  2:931435/101‑1 (MQ=255)
                                                     cGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGAt                                            <  1:1445584/101‑1 (MQ=255)
                                                     cGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGAt                                            >  2:1634919/1‑101 (MQ=255)
                                                      ggTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATc                                           >  1:3634536/1‑101 (MQ=255)
                                                        ttGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGa                                         <  2:1347826/101‑1 (MQ=255)
                                                        ttGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGa                                         <  2:2578364/101‑1 (MQ=255)
                                                         tggcgctgggcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAc                                        <  1:3332116/92‑1 (MQ=25)
                                                         tGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAc                                        <  1:2248245/101‑1 (MQ=255)
                                                         tGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAc                                        >  1:3217795/1‑101 (MQ=255)
                                                         tGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAc                                        <  2:220249/101‑1 (MQ=255)
                                                             cgcgcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACgcgc                                    >  1:2870623/1‑101 (MQ=255)
                                                               cgcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAt                                  <  1:2640547/101‑1 (MQ=255)
                                                               cgcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAt                                  <  1:447355/101‑1 (MQ=255)
                                                                gcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAtt                                 >  2:2790533/1‑101 (MQ=255)
                                                                  gcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTc                               <  2:3986354/101‑1 (MQ=255)
                                                                  gcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTc                               <  2:319079/101‑1 (MQ=255)
                                                                  gcgGGCTGCGGCGGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTc                               <  2:1058977/101‑1 (MQ=255)
                                                                   cgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCg                              >  1:2735204/1‑101 (MQ=255)
                                                                    gggCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGt                             <  1:6716/101‑1 (MQ=255)
                                                                    atcttGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGt                             <  1:4291532/97‑1 (MQ=255)
                                                                      gttccGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgc                           <  2:180764/97‑1 (MQ=255)
                                                                      gCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgc                           <  1:2178563/101‑1 (MQ=255)
                                                                      gCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGGCTCGATGCTGAAGATCGACGCGCATTTCGTgg                           >  2:3375270/1‑100 (MQ=255)
                                                                        tGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGtagat                         >  2:4291263/1‑97 (MQ=255)
                                                                           ggCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGc                      >  2:4017387/1‑101 (MQ=255)
                                                                           ggCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGc                      <  1:922887/101‑1 (MQ=255)
                                                                           ggCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGc                      >  2:3944054/1‑101 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                   >  2:4421988/1‑101 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                   >  1:878700/1‑101 (MQ=255)
                                                                               ggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTTAAGATCGACGCGCATTTCGTGCGCGGCATCg                  <  1:3981858/101‑1 (MQ=255)
                                                                               ggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCg                  <  1:1210434/101‑1 (MQ=255)
                                                                                gctaTCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGa                 <  2:3330237/98‑1 (MQ=255)
                                                                                 cGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg                >  1:1946423/1‑101 (MQ=255)
                                                                                 cGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg                <  1:2790739/101‑1 (MQ=255)
                                                                                  gATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAgg               <  1:3745212/101‑1 (MQ=255)
                                                                                  gATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAgg               >  1:2081115/1‑101 (MQ=255)
                                                                                  gATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAgg               >  2:3692869/1‑101 (MQ=255)
                                                                                   aTCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGc              >  1:2576311/1‑101 (MQ=255)
                                                                                   aTCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATAGACGCGCATTTCGTGCGCGGCATCGAGGc              >  1:422126/1‑101 (MQ=255)
                                                                                     cgacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATATACGCGCATTTCTTGCGCGGCATCGAGGCgg            >  2:2689348/1‑101 (MQ=25)
                                                                                      gacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGa           >  1:1801369/1‑101 (MQ=255)
                                                                                       acgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGAc          <  1:3911933/101‑1 (MQ=255)
                                                                                        cgaTTTCGGCCCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGAcc         <  1:2034469/101‑1 (MQ=255)
                                                                                        cgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGAcc         >  2:2140144/1‑101 (MQ=255)
                                                                                         gaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCg        >  1:3504748/1‑101 (MQ=255)
                                                                                           tttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCa      <  1:4045711/101‑1 (MQ=255)
                                                                                           tttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCCGCATCGAGGCGGACCGCa      <  2:4205961/101‑1 (MQ=255)
                                                                                            ttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCaa     <  1:2196877/101‑1 (MQ=255)
                                                                                            ttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCaa     <  2:1965649/101‑1 (MQ=255)
                                                                                             tCGGCCCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAAc    <  2:878615/101‑1 (MQ=255)
                                                                                             tCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAAc    <  2:1658556/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACGC  >  NC_014034/861544‑861739

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG