Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,478 G→C 18.3% intergenic (+290/‑641) RCAP_rcc01306 → / → RCAP_rcc01307 Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,395,4780G→C18.3% 227.5 ‑0.6 ‑0.0496intergenic (+290/‑641)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (42/49):  ref base (201/204):  total (243/253)
Fisher's exact test strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

Only 100 of 515 total aligned reads displayed.

CGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCC  >  NC_014034/1395378‑1395578
                                                                                                       |                                                                                                    
cgGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTgg                                                                                                      <  2:8119307/101‑1 (MQ=255)
 gggggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGc                                                                                                     >  1:5773559/1‑101 (MQ=255)
   gggCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCgg                                                                                                   >  1:819161/1‑101 (MQ=255)
        gggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGAc                                                                                              <  1:802163/100‑1 (MQ=255)
         ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGCCa                                                                                             >  2:7829844/1‑101 (MQ=255)
         ggATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACa                                                                                             >  1:376701/1‑101 (MQ=255)
          gATCGAAACTGCCCCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAg                                                                                            <  2:4066409/101‑1 (MQ=255)
           aTCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAgg                                                                                           >  2:6675991/1‑101 (MQ=255)
              gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTg                                                                                        >  2:2899299/1‑101 (MQ=255)
              gAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTg                                                                                        >  1:8533355/1‑101 (MQ=255)
               aaaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTgg                                                                                       <  1:2384678/101‑1 (MQ=255)
                aaCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGt                                                                                      <  1:6777684/101‑1 (MQ=255)
                 accgcccccGCTGGCCGTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTa                                                                                     <  2:8473692/95‑1 (MQ=255)
                 aCTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTa                                                                                     >  1:4535561/1‑101 (MQ=255)
                  cTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTaa                                                                                    <  1:3965134/101‑1 (MQ=255)
                  cTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTaa                                                                                    <  1:863149/101‑1 (MQ=255)
                   tGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAAc                                                                                   <  2:3458749/101‑1 (MQ=255)
                    gCACCCGCTGGCCCTGCCGGTGATCGGTTATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGTAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                  >  1:6525000/1‑101 (MQ=255)
                    gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                  >  1:3497587/1‑101 (MQ=255)
                    gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                  <  2:3694929/101‑1 (MQ=255)
                    gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                  <  1:1221939/101‑1 (MQ=255)
                    gCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACg                                                                                  >  1:5282248/1‑101 (MQ=255)
                     cACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGGAACgt                                                                                 >  1:4129292/1‑101 (MQ=255)
                      ccccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACgtg                                                                                <  1:6321692/100‑1 (MQ=255)
                       cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACgtgt                                                                               <  2:630624/101‑1 (MQ=255)
                       cccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACgtgt                                                                               <  1:18147/101‑1 (MQ=255)
                        ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                              <  2:129209/101‑1 (MQ=255)
                        ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTc                                                                              <  1:1006881/101‑1 (MQ=255)
                        ccGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCCGGGAACgggtg                                                                              >  1:4270898/1‑100 (MQ=255)
                          gCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGc                                                                            >  1:4045051/1‑101 (MQ=255)
                            tGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTc                                                                          <  1:4426694/101‑1 (MQ=255)
                               cccTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccc                                                                       >  2:7607735/1‑101 (MQ=255)
                                 cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                     <  2:1821111/101‑1 (MQ=255)
                                 cTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTcccccc                                                                     <  2:8602628/101‑1 (MQ=255)
                                     cggtgatcggtgatCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGTgg                                                                 >  2:58239/1‑101 (MQ=255)
                                       gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                               >  1:2862823/1‑101 (MQ=255)
                                       gtgatcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCgggg                                                               <  2:9157280/101‑1 (MQ=255)
                                          atcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCg                                                            >  1:7727501/1‑101 (MQ=255)
                                           tcggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGt                                                           >  2:2182117/1‑101 (MQ=255)
                                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                         <  2:4610049/101‑1 (MQ=255)
                                             ggtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTgg                                                         >  2:4823053/1‑101 (MQ=255)
                                             cgctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGa                                                        >  2:5002965/11‑101 (MQ=2)
                                              gtgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGa                                                        <  2:5773324/101‑1 (MQ=255)
                                              gctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                       <  1:902569/92‑1 (MQ=2)
                                               tgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                       <  1:7220204/101‑1 (MQ=255)
                                               tgatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                       <  2:2671251/101‑1 (MQ=255)
                                               ggatcAGCGGTCGCGCCCCGCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGAt                                                       <  1:7387187/100‑1 (MQ=255)
                                               ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      >  1:5335307/9‑101 (MQ=0)
                                               ctggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGGGGATc                                                      >  1:2105526/9‑101 (MQ=2)
                                                gatcAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                      >  1:6036381/1‑101 (MQ=255)
                                                 ggacgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                    >  1:4587074/7‑101 (MQ=2)
                                                   gcgtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  <  2:1596254/97‑1 (MQ=2)
                                                    aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  <  1:6923366/101‑1 (MQ=255)
                                                    aGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  >  2:7278805/1‑101 (MQ=255)
                                                     gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGGTGGAAGCa                                                >  2:4265272/3‑101 (MQ=2)
                                                      cGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                >  2:7699325/1‑101 (MQ=255)
                                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                              <  1:3293409/101‑3 (MQ=2)
                                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACagc                                              >  2:7906430/1‑99 (MQ=2)
                                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                               >  2:3966104/1‑101 (MQ=255)
                                                       ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGGGGCTCGAAACAt                                               >  2:8976857/1‑101 (MQ=255)
                                                        gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtc                                              <  1:1593071/101‑1 (MQ=255)
                                                         gcGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                             <  1:2671250/100‑1 (MQ=255)
                                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                           >  2:1667243/1‑99 (MQ=2)
                                                          cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                            <  2:5466039/101‑1 (MQ=255)
                                                            cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                         >  2:6072326/1‑97 (MQ=2)
                                                            cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                          <  2:3220184/101‑1 (MQ=255)
                                                             tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCggg                                         >  2:9056260/1‑100 (MQ=255)
                                                             tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                         <  2:8738134/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                        <  1:6195585/101‑1 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                        >  2:6793337/1‑101 (MQ=255)
                                                              ccccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                        <  2:1378752/101‑1 (MQ=255)
                                                                ccccGCGGGGGCGTGGATCGAAACTACACCG..GG.GCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaatgg                                     <  1:484645/101‑9 (MQ=2)
                                                                ccccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGcc                                      <  2:2313766/101‑1 (MQ=255)
                                                                  ccGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCGGt                                    <  1:9169591/101‑1 (MQ=255)
                                                                  ccGCGGGGGCGTGGATCGAAACACCCCGGACCTGTCGGCGG.CGA.ACACCGATACCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtcgatgatcaa                                     <  1:6389223/101‑10 (MQ=0)
                                                                     cgGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                 >  1:3019490/1‑101 (MQ=255)
                                                                      gggggCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                >  1:2253455/1‑101 (MQ=255)
                                                                      gggggCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                >  1:2490001/1‑101 (MQ=255)
                                                                        gggCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCg                              >  1:7569269/1‑101 (MQ=255)
                                                                          gCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTc                            <  1:8314668/101‑1 (MQ=255)
                                                                           cGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTcc                           <  1:5122420/101‑1 (MQ=255)
                                                                           cGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCCGCGCTCGTcc                           >  1:4945909/1‑101 (MQ=255)
                                                                            gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                          <  2:7411876/101‑1 (MQ=255)
                                                                            gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                          >  1:548174/1‑101 (MQ=255)
                                                                            gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCa                          >  2:8798040/1‑101 (MQ=255)
                                                                            gTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATAGAAACATCTCGGTGCCGTCGCTCGTCCa                          >  1:7962610/1‑101 (MQ=255)
                                                                              ggATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGc                        <  1:1329052/101‑1 (MQ=255)
                                                                              ggATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTCGGTGCCGTCCCTCGTCCCGc                        >  1:873700/1‑101 (MQ=255)
                                                                                  cGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGc                    <  2:5161117/101‑1 (MQ=255)
                                                                                   gAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGcc                   >  2:1012674/1‑101 (MQ=255)
                                                                                       cATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGcc               >  2:2022417/1‑100 (MQ=255)
                                                                                         tCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcg             <  2:705073/101‑1 (MQ=255)
                                                                                          cgccg..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgt            >  1:358467/1‑101 (MQ=255)
                                                                                          cgccg..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgt            <  1:5878758/101‑1 (MQ=255)
                                                                                           gccg..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc           >  1:138237/1‑101 (MQ=255)
                                                                                           gccg..cT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCtcgtc           <  2:643922/101‑1 (MQ=255)
                                                                                                  t.ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc      >  1:6572113/1‑101 (MQ=255)
                                                                                                  t.ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcc      >  2:3634490/1‑101 (MQ=255)
                                                                                                      ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTccccc   <  1:7871644/101‑1 (MQ=255)
                                                                                                       gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTcccccc  <  1:2797691/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
CGGGGGCGTGGATCGAAACTGCACCCGCTGGCCCTGCCGGTGATCGGTGATCAGCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCG..CT.TTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCGTCCAGCCGGCCTGCTCGTCGCTCCCCCC  >  NC_014034/1395378‑1395578

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG