Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,728,522 T→G 100% R635R (CGA→CGC RCAP_rcc02539 ← diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,728,5220T→G100.0% 760.5 215R635R (CGA→CGCRCAP_rcc02539diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (87/128):  ref base (0/0):  total (87/128)

Only 100 of 243 total aligned reads displayed.

GGGATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCTCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCGCGCGCT  >  NC_014034/2728425‑2728619
                                                                                                 |                                                                                                 
gggATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCgg                                                                                                <  2:233177/101‑1 (MQ=255)
gggATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCgg                                                                                                >  2:1541287/1‑101 (MQ=255)
 ggATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGGCGCggg                                                                                               >  1:4780985/1‑100 (MQ=255)
  gATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCgggg                                                                                              >  1:138746/1‑99 (MQ=255)
   aTTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGc                                                                                             <  1:897504/101‑1 (MQ=255)
    tttGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCa                                                                                            <  1:1030311/101‑1 (MQ=255)
     ttGACGGCGATGCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAg                                                                                           <  2:1188930/101‑1 (MQ=255)
     ttGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAg                                                                                           >  1:2372235/1‑101 (MQ=255)
       gACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGaa                                                                                         >  2:1806061/1‑101 (MQ=255)
       gACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGaa                                                                                         <  1:2939252/101‑1 (MQ=255)
        ccGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGaaa                                                                                        <  1:4985142/100‑1 (MQ=255)
        aCGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGaaa                                                                                        >  1:2717519/1‑101 (MQ=255)
          gCCGATCCGGGGCGAACGGGAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGAGGGAGGCGCGGAGCACAAAAc                                                                                      >  2:468997/1‑101 (MQ=255)
            cGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAAcgc                                                                                    >  2:358779/1‑101 (MQ=255)
             gATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAAcgcg                                                                                   <  2:3795907/101‑1 (MQ=255)
              aTCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCACCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAAAgcgc                                                                                  >  1:4050805/1‑101 (MQ=255)
                ccGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCaa                                                                                >  2:1074041/1‑101 (MQ=255)
                ccGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCaa                                                                                <  2:2726738/101‑1 (MQ=255)
                 cGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCaaa                                                                               >  2:3143993/1‑101 (MQ=255)
                 cGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCaaa                                                                               <  1:1100346/101‑1 (MQ=255)
                  gAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCaaaa                                                                              >  1:3610337/1‑101 (MQ=255)
                  gAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCcaaa                                                                              >  2:1870386/1‑101 (MQ=255)
                  gAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCaaaa                                                                              <  1:4623430/101‑1 (MQ=255)
                   aGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAAt                                                                             <  2:1637047/101‑1 (MQ=255)
                   aGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAAt                                                                             <  2:178561/101‑1 (MQ=255)
                   aGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAAt                                                                             <  2:3472864/101‑1 (MQ=255)
                   aGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAAt                                                                             <  2:162676/101‑1 (MQ=255)
                   aGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAAt                                                                             <  2:5028096/101‑1 (MQ=255)
                    ggCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATc                                                                            <  1:4264479/101‑1 (MQ=255)
                    ggCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATc                                                                            <  1:2930127/101‑1 (MQ=255)
                      cGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcg                                                                          >  2:753274/1‑101 (MQ=255)
                      cGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcg                                                                          <  2:2363998/101‑1 (MQ=255)
                      cGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcg                                                                          <  2:3005256/101‑1 (MQ=255)
                       gAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCcgc                                                                         >  1:4371741/1‑101 (MQ=255)
                       gAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcgc                                                                         <  2:4490770/101‑1 (MQ=255)
                       gAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcgc                                                                         >  1:5020483/1‑101 (MQ=255)
                       gAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGGCGCGGGGCAGAAAACGCGCAAAAACcgc                                                                         >  1:2277644/1‑101 (MQ=255)
                           ggCAAGGCGCGAAAGCGCCACAGTCGATTTCTGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGca                                                                     >  1:4488689/1‑101 (MQ=255)
                            gCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcat                                                                    >  2:2028441/1‑101 (MQ=255)
                             cAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatc                                                                   >  1:1988963/1‑101 (MQ=255)
                               aGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGcatcca                                                                 >  1:1557092/1‑101 (MQ=255)
                               aGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatcca                                                                 <  1:3165751/101‑1 (MQ=255)
                                ggCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGcatccat                                                                >  1:220307/1‑101 (MQ=255)
                                ggCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatccat                                                                <  1:2519324/101‑1 (MQ=255)
                                ggCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatccat                                                                <  2:2826229/101‑1 (MQ=255)
                                 gcgcGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatccatc                                                               <  1:3865296/101‑1 (MQ=255)
                                  cgcgAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCg                                                              <  1:4202921/101‑1 (MQ=255)
                                      aaaGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGCATCCCTCGGGCt                                                          >  1:3888906/1‑101 (MQ=255)
                                       aaGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTc                                                         >  2:1503004/1‑101 (MQ=255)
                                        aGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCg                                                        <  1:2183698/101‑1 (MQ=255)
                                          cgcCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGAc                                                      <  1:1638974/101‑1 (MQ=255)
                                              acaGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACcaca                                                  <  1:1867891/101‑1 (MQ=255)
                                               caGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAg                                                 <  1:3366933/101‑1 (MQ=255)
                                               caGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAg                                                 <  1:4732174/101‑1 (MQ=255)
                                                aGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACCGa                                                >  2:4395239/1‑101 (MQ=255)
                                                 gTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGAt                                               <  1:68246/101‑1 (MQ=255)
                                                  tCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATa                                              <  1:1001286/101‑1 (MQ=255)
                                                  tCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATa                                              <  1:3307084/101‑1 (MQ=255)
                                                   cGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAg                                             <  1:1593540/101‑1 (MQ=255)
                                                    gATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGGATCCATCGGGCTCGACCACAGATAGc                                            >  2:2462014/1‑101 (MQ=255)
                                                    gATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGc                                            <  1:2779986/101‑1 (MQ=255)
                                                     aTTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGcc                                           >  2:2591441/1‑101 (MQ=255)
                                                        tCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTg                                        >  1:3111967/1‑101 (MQ=255)
                                                         ccGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGc                                       >  2:277141/1‑101 (MQ=255)
                                                            gTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCg                                    <  2:3955783/101‑1 (MQ=255)
                                                              ccTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGAc                                  >  2:3378982/1‑101 (MQ=255)
                                                              ccTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGAc                                  <  1:1915332/101‑1 (MQ=255)
                                                               cTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCCCAGATAGCCCTGCCCGACa                                 >  2:4037223/1‑101 (MQ=255)
                                                                tGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGAcct                                >  2:4949914/1‑99 (MQ=255)
                                                                 gTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATc                               <  1:2057138/101‑1 (MQ=255)
                                                                 gTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATc                               <  1:2633917/101‑1 (MQ=255)
                                                                 gTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGATCGACCACAGCTAACCCTGCCCGACATc                               >  2:2480118/1‑101 (MQ=255)
                                                                  tCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCg                              <  1:3795626/101‑1 (MQ=255)
                                                                  tCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCg                              <  1:169827/101‑1 (MQ=255)
                                                                  tCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCg                              <  1:2823734/101‑1 (MQ=255)
                                                                  tCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATAGGGCTCGACCACAGATAGGCCTGCACGACATCg                              >  1:4359217/1‑101 (MQ=255)
                                                                    ccccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGCCATCGAc                            >  1:2839423/1‑101 (MQ=255)
                                                                    ccccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGAc                            <  1:3392450/101‑1 (MQ=255)
                                                                    ccccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGAc                            <  1:1466558/101‑1 (MQ=255)
                                                                     cccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACg                           <  1:57958/101‑1 (MQ=255)
                                                                     cccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACg                           >  1:3995072/1‑101 (MQ=255)
                                                                      ccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGc                          <  1:2482886/101‑1 (MQ=255)
                                                                        gCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGccc                        >  1:4138983/1‑101 (MQ=255)
                                                                           cTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGa                     >  2:1388754/1‑101 (MQ=255)
                                                                            tCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGAt                    >  2:2637964/1‑101 (MQ=255)
                                                                            tCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGAt                    <  2:881734/101‑1 (MQ=255)
                                                                              aaCCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCg                  >  1:4440598/1‑101 (MQ=255)
                                                                                ccGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGATATCGACGCCCAGATCGAc                <  2:625588/101‑1 (MQ=255)
                                                                                ccGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGAc                <  2:3305106/101‑1 (MQ=255)
                                                                                ccGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGAc                <  2:218168/101‑1 (MQ=255)
                                                                                 cGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGAcc               <  2:4137940/101‑1 (MQ=255)
                                                                                   cTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAg             <  2:1756148/101‑1 (MQ=255)
                                                                                    tGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGc            <  2:1209866/101‑1 (MQ=255)
                                                                                      cggcggGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCAt          >  2:1289744/1‑101 (MQ=255)
                                                                                       ggcgggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATc         <  1:829948/101‑1 (MQ=255)
                                                                                         cgggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATcgc       >  2:2851426/1‑101 (MQ=255)
                                                                                         cgggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATcgc       >  1:751106/1‑101 (MQ=255)
                                                                                          gggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCATAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATcgcg      <  2:4914505/101‑1 (MQ=255)
                                                                                             cgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCgcgcgc   >  2:3438833/1‑101 (MQ=255)
                                                                                             cgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCgcgcgc   <  2:335138/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
GGGATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCTCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCGCGCGCT  >  NC_014034/2728425‑2728619

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG