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% 1056.1 290A809T (GCC→ACC) RCAP_rcc00783diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (133/153):  ref base (0/0):  total (133/157)

Only 100 of 299 total aligned reads displayed.

TGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACGC  >  NC_014034/861544‑861739
                                                                                                  |                                                                                                 
tggccGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                 >  2:1556801/4‑101 (MQ=255)
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                 <  2:3447936/101‑1 (MQ=255)
 gACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCg                                                                                                <  1:1155762/101‑1 (MQ=255)
    cGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCt                                                                                             >  1:1422275/1‑101 (MQ=255)
     gAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTa                                                                                            >  1:404161/1‑101 (MQ=255)
      aGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTAt                                                                                           <  2:3064087/101‑1 (MQ=255)
       ggATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTAtt                                                                                          <  1:1625772/101‑1 (MQ=255)
        gATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTc                                                                                         >  2:887390/1‑101 (MQ=255)
        gATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTc                                                                                         <  1:942897/101‑1 (MQ=255)
         aTCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTcc                                                                                        >  2:1656483/1‑101 (MQ=255)
           cAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTg                                                                                      >  1:611681/1‑101 (MQ=255)
            aGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGt                                                                                     <  1:2230012/101‑1 (MQ=255)
            aGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGt                                                                                     >  1:4290454/1‑101 (MQ=255)
              cTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCt                                                                                   <  1:2342151/101‑1 (MQ=255)
              cTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCt                                                                                   >  1:5992476/1‑101 (MQ=255)
                 gCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGc                                                                                <  2:1276849/101‑1 (MQ=255)
                    gagcGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCGGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                             <  2:3332998/98‑1 (MQ=255)
                    gACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                             >  1:1250112/1‑101 (MQ=255)
                    gACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                             >  2:240478/1‑101 (MQ=255)
                    gACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                             >  1:6890617/1‑101 (MQ=255)
                    gACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                             >  1:168513/1‑101 (MQ=255)
                     gccGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTAt                                                                            <  2:1848525/100‑1 (MQ=255)
                      ccGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATc                                                                           <  2:5357834/101‑1 (MQ=255)
                       cGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCt                                                                          <  2:298158/101‑1 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                         <  2:756651/101‑1 (MQ=255)
                          gTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGgc                                                                       <  2:3262469/101‑1 (MQ=255)
                               gggCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGc                                                                  >  1:67403/1‑101 (MQ=255)
                                  cgcgGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGc                                                               <  2:6045351/101‑1 (MQ=255)
                                  cgcgGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGc                                                               <  2:4837731/101‑1 (MQ=255)
                                  cgcgGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGc                                                               <  2:5141722/101‑1 (MQ=255)
                                  cgcgGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGc                                                               <  2:6852625/101‑1 (MQ=255)
                                   gcgGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGcc                                                              <  1:3465270/101‑1 (MQ=255)
                                     ggTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCgg                                                            >  2:3865293/1‑101 (MQ=255)
                                     ggTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCgg                                                            <  1:5052079/101‑1 (MQ=255)
                                     ggTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCgg                                                            >  1:1111152/1‑101 (MQ=255)
                                        cGATGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtct                                                         >  2:5937866/1‑101 (MQ=255)
                                        cGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtct                                                         <  1:385925/101‑1 (MQ=255)
                                         gAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtctc                                                        >  1:1989032/1‑101 (MQ=255)
                                          aGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCg                                                       >  1:5146550/1‑101 (MQ=255)
                                          aGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCg                                                       >  2:1973785/1‑101 (MQ=255)
                                          aGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCg                                                       >  1:1509043/1‑101 (MQ=255)
                                          aGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCg                                                       <  2:4498868/101‑1 (MQ=255)
                                             cgcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATg                                                    <  2:508604/101‑1 (MQ=255)
                                               cTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCCCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCt                                                  <  1:695753/101‑1 (MQ=255)
                                               cTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCt                                                  <  2:7039029/101‑1 (MQ=255)
                                               cTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCt                                                  >  1:3769698/1‑101 (MQ=255)
                                                ttGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGGGCGGCTGCCGGTCTCGATGccg                                                 >  1:663867/1‑99 (MQ=255)
                                                ttGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTg                                                 <  2:4397825/101‑1 (MQ=255)
                                                ttGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTg                                                 <  2:4675230/101‑1 (MQ=255)
                                                 tGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGa                                                <  1:590073/101‑1 (MQ=255)
                                                  gggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGaa                                               <  1:7084047/101‑1 (MQ=255)
                                                  gggCGGTGGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGaa                                               <  1:2930789/101‑1 (MQ=255)
                                                   ggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAg                                              >  1:4564779/1‑101 (MQ=255)
                                                    gCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGa                                             <  1:3044668/101‑1 (MQ=255)
                                                     gggTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGAt                                            <  2:1928181/100‑1 (MQ=255)
                                                     cGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGAt                                            <  1:4698574/101‑1 (MQ=255)
                                                     cGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGAt                                            <  2:4275598/101‑1 (MQ=255)
                                                       gTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCg                                          <  1:4467715/101‑1 (MQ=255)
                                                         tGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAc                                        >  1:4899652/1‑101 (MQ=255)
                                                          gcgcgcgcgcGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAcg                                       >  1:2088460/1‑101 (MQ=255)
                                                          gcgcgcgcgcGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAcg                                       >  1:3360752/1‑101 (MQ=255)
                                                              gcgcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATNCTGAAGATCGACGCGCa                                   >  2:2001304/1‑101 (MQ=255)
                                                                gcgcgCGCTGGGGGTGGCGATCGCTGA.TTCGGCCCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAttt                                <  2:2272011/101‑1 (MQ=25)
                                                                 cgctGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAttt                                <  1:5759303/101‑1 (MQ=255)
                                                                 cgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAttt                                <  2:5445557/101‑1 (MQ=255)
                                                                 cgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAttt                                <  2:6675606/101‑1 (MQ=255)
                                                                  gcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTc                               <  1:1003773/101‑1 (MQ=255)
                                                                     ggCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTg                            >  1:7048773/1‑101 (MQ=255)
                                                                       cTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgcg                          <  1:3271541/101‑1 (MQ=255)
                                                                        tGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgcgc                         >  1:393887/1‑101 (MQ=255)
                                                                        tGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgcgc                         >  2:6818406/1‑101 (MQ=255)
                                                                        tGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGGCTCGATGCTGAAGATCGACGCGCATTTCGTgcgc                         >  2:2100018/1‑101 (MQ=255)
                                                                             cTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCAt                    <  1:4768053/101‑1 (MQ=255)
                                                                             cTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCAt                    <  1:647495/101‑1 (MQ=255)
                                                                             cTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCAt                    <  1:2649303/101‑1 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCCCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                   <  1:4835956/101‑1 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                   <  1:1310248/101‑1 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                   >  2:5613373/1‑101 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACAGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCCTGCGCGGCATc                   >  2:7081369/1‑101 (MQ=255)
                                                                               ggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCg                  >  1:1911003/1‑101 (MQ=255)
                                                                               ggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCg                  >  1:2894039/1‑101 (MQ=255)
                                                                                gCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGa                 <  1:5510734/101‑1 (MQ=255)
                                                                                 cgttgtccgaTTTCGGCCCCTGTTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg                <  1:5702348/94‑1 (MQ=25)
                                                                                 cGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg                <  2:2567278/101‑1 (MQ=255)
                                                                                 cGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg                <  2:1788156/101‑1 (MQ=255)
                                                                                 cGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg                >  2:6427259/1‑101 (MQ=255)
                                                                                  gATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAgg               >  2:3732502/1‑101 (MQ=255)
                                                                                  gATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGAGCATTTCGTGCGCGGCATCGAgg               >  2:4044056/1‑101 (MQ=255)
                                                                                   aTCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAggg              >  2:4002248/1‑100 (MQ=255)
                                                                                    tCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCg             <  2:1180868/101‑1 (MQ=255)
                                                                                      gccgaTTCCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGa           <  2:393883/99‑1 (MQ=255)
                                                                                      gacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCggc           >  2:1349035/1‑100 (MQ=255)
                                                                                      gacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGa           >  1:1382153/1‑101 (MQ=255)
                                                                                      gacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGa           >  1:5116039/1‑101 (MQ=255)
                                                                                      gacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGa           >  1:618807/1‑101 (MQ=255)
                                                                                       acgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGAc          >  1:65602/1‑101 (MQ=255)
                                                                                          aTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGc       >  1:3634037/1‑101 (MQ=255)
                                                                                            ttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCaa     >  1:244075/1‑101 (MQ=255)
                                                                                            ttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCaa     <  1:846661/101‑1 (MQ=255)
                                                                                              cGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACg   <  2:2893930/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACGC  >  NC_014034/861544‑861739

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG