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% 613.4 178R635R (CGA→CGCRCAP_rcc02539diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  new base (71/107):  ref base (0/0):  total (71/107)

Only 100 of 200 total aligned reads displayed.

CGGGATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCTCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCGCGCGCT  >  NC_014034/2728424‑2728619
                                                                                                  |                                                                                                 
   gATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAg                                                                                              <  1:145361/101‑1 (MQ=255)
   gATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAg                                                                                              >  2:2458401/1‑101 (MQ=255)
    aTTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGc                                                                                             <  1:1143078/101‑1 (MQ=255)
    aTTTGACGCCGATCCGAGGCGAACGGCAAGGCGAGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGCGc                                                                                             >  2:3503813/1‑101 (MQ=255)
     tttGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCa                                                                                            <  1:2941760/101‑1 (MQ=255)
      ttGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAg                                                                                           <  2:160364/101‑1 (MQ=255)
       tGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGa                                                                                          >  2:2122815/1‑101 (MQ=255)
       tGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGa                                                                                          <  1:751266/101‑1 (MQ=255)
       tGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGa                                                                                          >  2:1839593/1‑101 (MQ=255)
       tGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCATCAACCGCTGCGGCGGGCGTCGCGGAGCAGa                                                                                          <  1:2974341/101‑1 (MQ=255)
           gCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGGCGCGGGGCAGAAAAc                                                                                      >  2:1276448/1‑101 (MQ=255)
            ccGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAAcg                                                                                     <  1:4096574/101‑1 (MQ=255)
              gATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAAagcg                                                                                   >  1:3565584/1‑101 (MQ=255)
               aTCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACgcgc                                                                                  <  2:1164266/101‑1 (MQ=255)
                tCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCATAAAAcgcgga                                                                                 >  2:2040604/1‑99 (MQ=255)
                 ccGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCaa                                                                                <  2:839441/101‑1 (MQ=255)
                  cGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCCGAAAACGCGCaaa                                                                               >  1:215352/1‑101 (MQ=255)
                    aGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCGGAAAACGCGCAAAAt                                                                             >  1:2115523/1‑101 (MQ=255)
                     ggCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGACGTCGCGGAGCAGAAAACGCGCAAAATc                                                                            <  1:260784/101‑1 (MQ=255)
                       cGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCcg                                                                          >  2:3885237/1‑101 (MQ=255)
                       cGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcg                                                                          <  2:3295457/101‑1 (MQ=255)
                       cGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcg                                                                          <  2:4015024/101‑1 (MQ=255)
                         aaCGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCcgcg                                                                        <  2:2294668/101‑1 (MQ=255)
                          aCGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCgg                                                                       <  2:1053980/101‑1 (MQ=255)
                            ggCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGca                                                                     >  1:4063950/1‑101 (MQ=255)
                             gCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcat                                                                    <  1:24559/101‑1 (MQ=255)
                               aaGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCCCGGcatcc                                                                  >  2:1629720/1‑101 (MQ=255)
                                aGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatcca                                                                 >  2:2634245/1‑101 (MQ=255)
                                  gcgcGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGcatccatc                                                               <  1:2848191/101‑1 (MQ=255)
                                   cgcgAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCg                                                              <  2:3594617/101‑1 (MQ=255)
                                   cgcgAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCg                                                              <  2:913300/101‑1 (MQ=255)
                                       aaaGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGCATCCCTCGGGCt                                                          >  1:3759873/1‑101 (MQ=255)
                                       aaaGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCt                                                          >  2:3363872/1‑101 (MQ=255)
                                       aaaGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCt                                                          <  1:433733/101‑1 (MQ=255)
                                        aaGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTc                                                         <  1:3721066/101‑1 (MQ=255)
                                         aGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCCGGCTCg                                                        >  1:1683538/1‑101 (MQ=255)
                                          gcgcCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGa                                                       <  2:1818646/101‑1 (MQ=255)
                                           cgcCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATACGCGGCATCCATCGGGCTCgcc                                                      >  2:814977/1‑99 (MQ=255)
                                           cgcCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGAAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGAc                                                      <  2:145343/101‑1 (MQ=255)
                                               acaGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACcaca                                                  <  1:4185705/101‑1 (MQ=255)
                                                caGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAg                                                 >  2:48110/1‑101 (MQ=255)
                                                 aGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGa                                                <  2:2738574/101‑1 (MQ=255)
                                                   tCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACATATa                                              <  2:2566320/101‑1 (MQ=255)
                                                   tCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATa                                              <  2:430184/101‑1 (MQ=255)
                                                    cGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAAAGCGCCAAATCCGCGGCATCCCTCGGGCTCGACCACAGAAAg                                             >  2:2363169/1‑101 (MQ=255)
                                                    cGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAg                                             >  1:2790818/1‑101 (MQ=255)
                                                    cGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAg                                             >  1:3137668/1‑101 (MQ=255)
                                                    cGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAg                                             >  1:396815/1‑101 (MQ=255)
                                                    cGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAg                                             <  1:3810698/101‑1 (MQ=255)
                                                       tttCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGccc                                          >  1:681276/1‑101 (MQ=255)
                                                       tttCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGccc                                          <  2:2734995/101‑1 (MQ=255)
                                                        ttCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCt                                         >  2:3069180/1‑101 (MQ=255)
                                                        ttCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCt                                         <  1:3356778/101‑1 (MQ=255)
                                                         tCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTg                                        >  1:2366841/1‑101 (MQ=255)
                                                          ccGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGc                                       >  1:3648858/1‑101 (MQ=255)
                                                           cGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCCCGACCGCAGATAGCCCTGcc                                      >  2:2019188/1‑101 (MQ=255)
                                                           cGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGcc                                      <  1:1043063/101‑1 (MQ=255)
                                                           cGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGcc                                      <  2:1809127/101‑1 (MQ=255)
                                                            ggTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGccc                                     <  1:3364117/101‑1 (MQ=255)
                                                             gTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCg                                    <  2:2407744/101‑1 (MQ=255)
                                                             gTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCg                                    <  1:927484/101‑1 (MQ=255)
                                                              tCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGa                                   >  1:1221357/1‑101 (MQ=255)
                                                               ccTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCCAc                                  >  1:179082/1‑101 (MQ=255)
                                                               ccTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGAc                                  <  1:1628263/101‑1 (MQ=255)
                                                               ccTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGAc                                  <  1:1839713/101‑1 (MQ=255)
                                                               ccTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGGCGCGGGGCAGAAAACGCGCAAAATCCGCGGAATCCCTCGGGCTCGACCACAGATAGCCCTACCCGAc                                  >  1:723805/1‑101 (MQ=255)
                                                                cTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACa                                 <  1:2713438/101‑1 (MQ=255)
                                                                 tGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCCCAGATAGCCCTGCCCCAccc                                >  1:1496277/1‑99 (MQ=255)
                                                                  gTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATc                               <  2:2997399/101‑1 (MQ=255)
                                                                  gTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATc                               <  2:2934833/101‑1 (MQ=255)
                                                                   tCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGGGCAGAAAAAGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCg                              >  2:752613/1‑101 (MQ=255)
                                                                   tCCCCCGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCg                              <  2:1065152/101‑1 (MQ=255)
                                                                    cccccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGa                             <  1:3266502/101‑1 (MQ=255)
                                                                    cccccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGa                             <  1:2000041/101‑1 (MQ=255)
                                                                     ccccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCgcc                            >  2:4137384/1‑99 (MQ=255)
                                                                      cccGCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACg                           >  1:3636629/1‑101 (MQ=255)
                                                                         gCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGccc                        <  1:3051864/101‑1 (MQ=255)
                                                                         gCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGccc                        >  2:3638744/1‑101 (MQ=255)
                                                                         gCCCTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGccc                        >  1:1905004/1‑101 (MQ=255)
                                                                           ccTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAg                      <  2:2281874/101‑1 (MQ=255)
                                                                           ccTCAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCATGCCCGACATCGACGCCCAg                      <  1:3545595/101‑1 (MQ=255)
                                                                              cAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATc                   <  1:2276486/101‑1 (MQ=255)
                                                                              cAACCGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATc                   <  1:1077316/101‑1 (MQ=255)
                                                                                 ccGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTTCCCGACATCGACGCCCAGATCGAc                >  1:3439546/1‑101 (MQ=255)
                                                                                 ccGCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGAc                <  2:4163180/101‑1 (MQ=255)
                                                                                   gCTGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCa              <  1:863961/101‑1 (MQ=255)
                                                                                    cTGCGGCGGGCGTCGCGGGGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCACATCGACCAg             >  1:1746945/1‑101 (MQ=255)
                                                                                     tGCGGCTGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGc            <  2:848380/101‑1 (MQ=255)
                                                                                     tGCGGCGGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGc            <  1:1590346/101‑1 (MQ=255)
                                                                                      gcggcgGGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCa           >  2:3188788/1‑101 (MQ=255)
                                                                                       cggcggGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCAt          <  2:1921019/101‑1 (MQ=255)
                                                                                       cggcggGCGTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCAt          <  2:2553714/101‑1 (MQ=255)
                                                                                        ggcgggcgTCGCGGAGCCGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATc         <  2:3121245/101‑1 (MQ=255)
                                                                                        ggcgggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATc         <  2:706885/101‑1 (MQ=255)
                                                                                        ggcgggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATc         <  2:2562544/101‑1 (MQ=255)
                                                                                           gggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATcgcg      >  2:1759168/1‑101 (MQ=255)
                                                                                           gggcgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATcgcg      >  2:2953098/1‑101 (MQ=255)
                                                                                              cgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCgcgcgc   >  2:2812319/1‑101 (MQ=255)
                                                                                              cgTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCgcgcgc   >  2:3251237/1‑101 (MQ=255)
                                                                                               gTCGCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCGCGCGCt  <  2:3759616/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
CGGGATTTGACGCCGATCCGAGGCGAACGGCAAGGCGCGAAAGCGCCACAGTCGATTTCCGGTCCTGTCCCCCGCCCTCAACCGCTGCGGCGGGCGTCTCGGAGCAGAAAACGCGCAAAATCCGCGGCATCCATCGGGCTCGACCACAGATAGCCCTGCCCGACATCGACGCCCAGATCGACCAGCATCGCGCGCT  >  NC_014034/2728424‑2728619

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG