Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,669 G→C 5.3% intergenic (+481/‑450) 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,6690G→C 5.3% 21.0 ‑1.4 ‑0.6224intergenic (+481/‑450)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/4):  ref base (103/109):  total (111/113)
Fisher's exact test strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.46e-01

Only 100 of 251 total aligned reads displayed.

CGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCT  >  NC_014034/1395569‑1395769
                                                                                                    |                                                                                                       
cgCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTg                                                                                                         <  1:2760293/101‑1 (MQ=255)
 gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                        >  1:1071265/1‑101 (MQ=255)
 gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                        >  2:3632111/1‑101 (MQ=255)
  cTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCa                                                                                                       <  2:3769001/101‑1 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                     >  1:3264555/1‑101 (MQ=255)
      ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCCTGTa                                                                                                   >  2:4003822/1‑101 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                  <  1:2772661/101‑1 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                  >  1:157611/1‑101 (MQ=255)
       cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                  <  1:618687/101‑1 (MQ=255)
         cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATCTGCATGTAGTc                                                                                                >  1:954095/1‑101 (MQ=255)
         cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                                >  2:3172023/1‑101 (MQ=255)
         cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                                <  1:228663/101‑1 (MQ=255)
          gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                               <  1:2104600/101‑1 (MQ=255)
          gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCTGCATGTAGTcg                                                                                               >  2:2137601/1‑101 (MQ=255)
           cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg                                                                                              >  1:809245/1‑101 (MQ=255)
             ggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcg                                                                                            >  1:68600/1‑101 (MQ=255)
                 cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCa                                                                                      >  1:3908300/1‑101 (MQ=255)
                  gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAg                                                                                     >  2:2535801/1‑101 (MQ=255)
                   tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGC.G.GCact                                                                                    >  1:453205/1‑99 (MQ=255)
                   tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGt                                                                                    >  1:3747002/1‑101 (MQ=255)
                    ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGtt                                                                                   >  2:2352148/1‑101 (MQ=255)
                    ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGtt                                                                                   >  1:1837522/1‑101 (MQ=255)
                     gATCGAAACCACCTTGCGGATTGGCACTTTCGCGGTGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTc                                                                                  <  2:527116/101‑1 (MQ=255)
                       gcGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGa                                                                                <  2:1487304/100‑1 (MQ=255)
                        cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.g                                                                              <  1:2218143/101‑1 (MQ=255)
                         gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.Gc                                                                             >  2:3483763/1‑101 (MQ=255)
                         gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.Gc                                                                             >  1:3982282/1‑101 (MQ=255)
                          aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.Gcc                                                                            <  1:3050159/101‑1 (MQ=255)
                           aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCg                                                                           <  1:741645/101‑1 (MQ=255)
                             cccccTTGCGGATTGTCCCTTTCGCGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCa                                                                         <  1:1325651/98‑1 (MQ=255)
                             ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCa                                                                         >  2:1269782/1‑101 (MQ=255)
                             ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCa                                                                         <  2:2897632/101‑1 (MQ=255)
                                ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTc                                                                      <  1:2441117/101‑1 (MQ=255)
                                 cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCg                                                                     >  1:578040/1‑101 (MQ=255)
                                  ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGc                                                                    <  2:3519710/101‑1 (MQ=255)
                                   tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCt                                                                   <  1:3440399/101‑1 (MQ=255)
                                    gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                  >  2:594554/1‑101 (MQ=255)
                                    ccGGATTGTCCCTTGCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                  <  2:4190437/100‑1 (MQ=255)
                                     cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                 >  1:4014569/1‑101 (MQ=255)
                                     cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                 <  1:3713683/101‑1 (MQ=255)
                                       gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                               <  2:3064234/101‑1 (MQ=255)
                                       gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                               >  1:2606808/1‑101 (MQ=255)
                                       gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccc                                                               <  1:427003/101‑1 (MQ=255)
                                        aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                              <  1:3627139/101‑1 (MQ=255)
                                          tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                            <  1:3636560/101‑1 (MQ=255)
                                            tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                          <  1:500662/101‑1 (MQ=255)
                                             cccTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                         <  2:4054374/99‑1 (MQ=255)
                                             cACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                         <  1:2981603/101‑1 (MQ=255)
                                               cTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                       <  1:2352298/101‑1 (MQ=255)
                                                tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                      >  2:712163/1‑101 (MQ=255)
                                                tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                      >  2:1719752/1‑101 (MQ=255)
                                                tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTACCCCCGCggggg                                                      >  2:2015694/1‑101 (MQ=255)
                                                    tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                   >  1:3853430/11‑101 (MQ=2)
                                                     cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                 <  2:2745645/101‑1 (MQ=255)
                                                      ggggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGa                                                >  1:1943772/1‑101 (MQ=255)
                                                       gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:3974134/101‑1 (MQ=255)
                                                       gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:896881/101‑1 (MQ=255)
                                                        gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                               >  2:3392779/7‑101 (MQ=2)
                                                         gCGAGGCGTCCCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                             <  2:1581173/101‑1 (MQ=255)
                                                           gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                           >  1:1869936/1‑101 (MQ=255)
                                                           cgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAA.GCTGGTCCC.GTCGCCGTCGGTCCCCCCGCGGGGGCGTGGATCGa                                            <  2:154366/96‑1 (MQ=2)
                                                            aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  2:472437/1‑101 (MQ=255)
                                                            aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  2:4071815/1‑101 (MQ=255)
                                                             tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  1:2140397/2‑101 (MQ=2)
                                                             aGCGTCGGCCCCCCCGGGGGGGCGTGGATCGAAACAGCTGCATGTAGTGGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         <  2:2459554/101‑1 (MQ=37)
                                                              gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                        <  2:331992/101‑1 (MQ=255)
                                                               cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                       >  1:311163/1‑101 (MQ=255)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                     <  1:1269875/101‑1 (MQ=255)
                                                                 gcGCCCCCCCCGGCGGGGGGTGGATCAACACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                     <  1:32959/100‑1 (MQ=25)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                  <  1:2681678/101‑1 (MQ=255)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                  <  2:1071177/101‑1 (MQ=255)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                  >  2:3039144/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                 >  2:3225306/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                 >  2:2939734/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                 >  1:4068018/1‑101 (MQ=255)
                                                                      ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                >  1:3098573/1‑101 (MQ=255)
                                                                      ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                <  2:1389129/101‑1 (MQ=255)
                                                                         cccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTGGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                             <  1:2566959/101‑1 (MQ=255)
                                                                            gcgGGGGCGTGGATCGAAACATCTGCAGGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                          <  2:1837398/101‑1 (MQ=255)
                                                                              gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                        >  2:771128/1‑101 (MQ=255)
                                                                               ggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCa                       <  2:3908662/101‑1 (MQ=255)
                                                                                  gCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                    >  1:3368259/1‑101 (MQ=255)
                                                                                  gCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                    <  1:2661211/101‑1 (MQ=255)
                                                                                   ggTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                   <  2:4245022/100‑1 (MQ=255)
                                                                                   cGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                   >  1:602825/1‑101 (MQ=255)
                                                                                    gtagATGGAAACATATGCATGTGGCCGGGG.G.GGAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                  <  1:1193797/98‑1 (MQ=25)
                                                                                    gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                  >  2:1547835/1‑101 (MQ=255)
                                                                                    gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                  >  1:2037351/1‑101 (MQ=255)
                                                                                     tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                 <  2:1743662/101‑1 (MQ=255)
                                                                                     tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                 >  2:229294/1‑101 (MQ=255)
                                                                                         tCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt             >  1:1122675/1‑101 (MQ=255)
                                                                                           gAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt           <  2:1869812/101‑1 (MQ=255)
                                                                                           gAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt           >  1:1438693/1‑101 (MQ=255)
                                                                                            aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg          >  2:1691194/1‑101 (MQ=255)
                                                                                              aCACCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc        <  2:2702594/101‑1 (MQ=255)
                                                                                               cATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCAGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg       >  1:3958474/1‑101 (MQ=255)
                                                                                                ctcggCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt      <  2:4141045/97‑1 (MQ=255)
                                                                                                aTCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt      >  2:1172522/1‑101 (MQ=255)
                                                                                                 tCTGGAGGTGGTGGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc     <  2:4341927/101‑1 (MQ=37)
                                                                                                    gCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt  >  1:102695/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
CGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCT  >  NC_014034/1395569‑1395769

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG