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% 691.6 190A809T (GCC→ACC) RCAP_rcc00783diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (91/97):  ref base (0/1):  total (91/99)

Only 100 of 203 total aligned reads displayed.

TGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACG  >  NC_014034/861544‑861738
                                                                                                  |                                                                                                
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCTCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                >  1:440779/1‑101 (MQ=255)
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                <  2:3434988/101‑1 (MQ=255)
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                >  1:168160/1‑101 (MQ=255)
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                <  1:2335160/101‑1 (MQ=255)
tGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCAcc                                                                                                >  1:1877290/1‑101 (MQ=255)
 gACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCg                                                                                               <  2:2642205/101‑1 (MQ=255)
 gACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCg                                                                                               <  2:2171036/101‑1 (MQ=255)
 gACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCg                                                                                               >  1:1948036/1‑101 (MQ=255)
   ccGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGc                                                                                             <  1:4128889/101‑1 (MQ=255)
    cGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCt                                                                                            >  1:2247754/1‑101 (MQ=255)
     gAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTa                                                                                           <  1:3302791/101‑1 (MQ=255)
     gAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTa                                                                                           <  1:2411777/101‑1 (MQ=255)
     gAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTa                                                                                           <  1:1986517/101‑1 (MQ=255)
       ggATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTAtt                                                                                         <  2:2927894/101‑1 (MQ=255)
          tCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCt                                                                                      <  1:2889776/101‑1 (MQ=255)
           cAGCTGGCCGACCGGGTGCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTg                                                                                     <  1:3923717/101‑1 (MQ=255)
              cTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCt                                                                                  >  1:4002714/1‑101 (MQ=255)
               tGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTg                                                                                 >  1:181789/1‑101 (MQ=255)
               tGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTg                                                                                 >  1:2262311/1‑101 (MQ=255)
                ggCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTgg                                                                                >  2:2792625/1‑101 (MQ=255)
                ggCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTgg                                                                                <  2:3746172/101‑1 (MQ=255)
                  ccGGCCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGcc                                                                              <  2:607739/101‑1 (MQ=255)
                  ccGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGcc                                                                              >  2:2344869/1‑101 (MQ=255)
                  ccGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGcc                                                                              >  2:1173449/1‑101 (MQ=255)
                   cGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCt                                                                             >  1:1580702/1‑101 (MQ=255)
                    gACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGGGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTa                                                                            >  2:2390316/1‑101 (MQ=255)
                     aCCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTAt                                                                           >  1:3735101/1‑101 (MQ=255)
                       cGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCt                                                                         <  1:2664153/101‑1 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                        <  1:1888660/101‑1 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                        >  1:1762114/1‑101 (MQ=255)
                        gggTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTg                                                                        <  2:3134989/101‑1 (MQ=255)
                             ccGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGcgcg                                                                   >  1:3880529/1‑101 (MQ=255)
                             ccGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGcgcg                                                                   <  2:2749101/101‑1 (MQ=255)
                              cGGGCGCTGTCGAGGCGCTTGGGCGGCTGCGCGCGCGCGGGCTGCGGCTGGCGACCCAGGACTTTGGCACCGGCTATTCCTGTCTCGGCTATCTGgcgccg                                                                  >  2:1588356/1‑99 (MQ=255)
                              cGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCgg                                                                  >  2:3284949/1‑101 (MQ=255)
                               gggCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGc                                                                 >  1:129911/1‑101 (MQ=255)
                                ggCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCt                                                                >  1:1546292/1‑101 (MQ=255)
                                ggCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCt                                                                <  2:411902/101‑1 (MQ=255)
                                      gTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGt                                                          <  1:1342226/101‑1 (MQ=255)
                                      gTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCCGt                                                          >  1:62718/1‑101 (MQ=255)
                                        cGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtct                                                        >  2:2571516/1‑101 (MQ=255)
                                        cGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtct                                                        >  2:553102/1‑101 (MQ=255)
                                        cGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtcn                                                        >  1:3101736/1‑100 (MQ=255)
                                         gAGGCGCTTGGGCGGTGGCGCGCGCGGGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGtctc                                                       <  2:1715304/101‑1 (MQ=255)
                                            gcgcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGcagatcggaag                                                    >  1:73633/1‑91 (MQ=37)
                                             cgcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATg                                                   <  1:1201314/101‑1 (MQ=255)
                                              gcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGCTGc                                                  >  1:139774/1‑101 (MQ=255)
                                              gcTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGc                                                  <  2:2505852/101‑1 (MQ=255)
                                                ttGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTg                                                >  2:1483642/1‑101 (MQ=255)
                                                ttGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCCATTCCTGTCTGGCCTATCTTGCGCGGCTGTCGGTGTCgggggtg                                                >  1:3417959/1‑95 (MQ=37)
                                                 tGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGa                                               <  1:67080/101‑1 (MQ=255)
                                                  gggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGaa                                              <  1:3586828/101‑1 (MQ=255)
                                                  gggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGaa                                              <  1:3640154/101‑1 (MQ=255)
                                                   ggCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTCAAg                                             >  2:1909289/1‑101 (MQ=255)
                                                    tcGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGa                                            <  1:3193108/100‑1 (MQ=255)
                                                     cGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGGCGGTCTCGATGCTGAAGaa                                           >  2:3192869/1‑100 (MQ=255)
                                                      ggTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTGGATGCTGAAGATc                                          >  2:2104056/1‑101 (MQ=255)
                                                        ttGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGa                                        >  1:795348/1‑101 (MQ=255)
                                                          gcgcgcgcgcGGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGCTGCTGAAGATCGAcg                                      >  2:146570/1‑101 (MQ=255)
                                                           cgcgcgcgcgGGCTGCGTCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGAcgc                                     <  2:3808231/101‑1 (MQ=255)
                                                              gcgcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCa                                  <  1:1975416/101‑1 (MQ=255)
                                                                gcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCTACGCGCAtt                                >  1:3381126/1‑101 (MQ=255)
                                                                gcgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAtt                                <  2:1031361/101‑1 (MQ=255)
                                                                 cgcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCAttt                               <  1:601516/101‑1 (MQ=255)
                                                                  gcgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTc                              <  2:2833683/101‑1 (MQ=255)
                                                                   cgGGCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCg                             >  1:1555404/1‑101 (MQ=255)
                                                                    gggCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGCTGCTGAAGATCGACGCGCATTTCGt                            >  1:3519436/1‑101 (MQ=255)
                                                                     ggCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTg                           <  1:1495782/101‑1 (MQ=255)
                                                                      gCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGGGCATTTCGTgc                          >  1:1959480/1‑101 (MQ=255)
                                                                      gCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgc                          >  2:4058817/1‑101 (MQ=255)
                                                                      gCTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgc                          >  1:3894558/1‑101 (MQ=255)
                                                                       cTGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGGGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgcg                         >  2:2370434/1‑101 (MQ=255)
                                                                        tGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgcgc                        >  2:3577996/1‑101 (MQ=255)
                                                                        tGCGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTgcgc                        >  2:853418/1‑101 (MQ=255)
                                                                          cGGCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGACGCTGAAGATCGACGCGCATTTCGTGCGCgg                      >  2:2849896/1‑101 (MQ=255)
                                                                            gCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCa                    >  1:1067714/1‑101 (MQ=255)
                                                                            gCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCa                    >  1:4032491/1‑101 (MQ=255)
                                                                            gCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCa                    >  1:45007/1‑101 (MQ=255)
                                                                            gCTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCa                    <  2:1273659/101‑1 (MQ=255)
                                                                            gCTGGCGAAAAACGATTACGGCACTCGCTATTCCTTTCTGGCCTATCTAGCGCTGCTGCCGGTATCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCa                    <  2:2090643/101‑1 (MQ=255)
                                                                             cTGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCAt                   >  1:680518/1‑101 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                  <  1:2692505/101‑1 (MQ=255)
                                                                              tGGCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                  <  2:3863229/101‑1 (MQ=255)
                                                                              gggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATc                  <  2:2992537/100‑1 (MQ=255)
                                                                               ggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCg                 >  2:2712611/1‑101 (MQ=255)
                                                                               ggCGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCg                 >  2:809403/1‑101 (MQ=255)
                                                                                 ggATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg               <  2:3009634/100‑1 (MQ=255)
                                                                                 cGATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAg               <  1:814812/101‑1 (MQ=255)
                                                                                  gATCGGCGATTTCGGCACCGGGTGTTCCTGNCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAgg              <  1:2104194/101‑1 (MQ=21)
                                                                                  gATCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAgg              <  2:1305066/101‑1 (MQ=255)
                                                                                   aTCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGc             >  2:3680039/1‑101 (MQ=255)
                                                                                   aTCGACGATTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGc             <  2:1513057/101‑1 (MQ=255)
                                                                                      gacgaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGGGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGa          <  2:1351816/101‑1 (MQ=255)
                                                                                         gaTTTCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCg       >  1:2807461/1‑101 (MQ=255)
                                                                                           tttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCa     >  2:432905/1‑101 (MQ=255)
                                                                                           tttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCa     <  2:2276612/101‑1 (MQ=255)
                                                                                            ttCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCaa    >  1:862992/1‑101 (MQ=255)
                                                                                             tCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAAc   <  1:30019/101‑1 (MQ=255)
                                                                                             tCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAAc   <  2:2734453/101‑1 (MQ=255)
                                                                                             tCGGCACCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAAc   >  1:1024988/1‑101 (MQ=255)
                                                                                                  |                                                                                                
TGACCGAGGATCAGCTGGCCGACCGGGTTCCGGGCGCGGTCGAGGCGCTTGGGCGGTTGCGCGCGCGCGGGCTGCGGCTGGCGATCGACGATTTCGGCGCCGGCTATTCCTGTCTGGCCTATCTGGCGCGGCTGCCGGTCTCGATGCTGAAGATCGACGCGCATTTCGTGCGCGGCATCGAGGCGGACCGCAACG  >  NC_014034/861544‑861738

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG