Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,672 T→A 6.0% intergenic (+484/‑447) 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,6720T→A 6.0% 21.8 ‑0.6 ‑0.0183intergenic (+484/‑447)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (7/4):  ref base (91/81):  total (98/85)
Fisher's exact test strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.56e-01

Only 100 of 203 total aligned reads displayed.

TCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCC  >  NC_014034/1395572‑1395771
                                                                                                    |                                                                                                         
tCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCAt                                                                                                           <  2:660415/101‑1 (MQ=255)
  cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                         <  2:1799200/101‑1 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                      <  1:346838/101‑1 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                      <  1:3522220/101‑1 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                      <  2:789832/101‑1 (MQ=255)
      cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                                     <  1:3356290/101‑1 (MQ=255)
        cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg                                                                                                   >  2:511482/1‑101 (MQ=255)
         gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg.c                                                                                                 >  1:787532/1‑101 (MQ=255)
         gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg.c                                                                                                 >  1:1613343/1‑101 (MQ=255)
         gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg.c                                                                                                 <  1:1941426/101‑1 (MQ=255)
            ggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.Cg..g.g                                                                                           <  1:2507234/101‑1 (MQ=255)
            ggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.Cg..g.g                                                                                           >  1:362220/1‑101 (MQ=255)
              cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCa                                                                                         <  2:2090907/101‑1 (MQ=255)
              cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCa                                                                                         <  1:630707/101‑1 (MQ=255)
                tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGt                                                                                       <  2:2142036/101‑1 (MQ=255)
                tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGt                                                                                       <  1:3496503/101‑1 (MQ=255)
                 ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGtt                                                                                      >  1:890488/1‑101 (MQ=255)
                 ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGtt                                                                                      >  2:4113940/1‑101 (MQ=255)
                  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTc                                                                                     <  1:3322825/101‑1 (MQ=255)
                  gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTc                                                                                     <  1:60362/101‑1 (MQ=255)
                    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGa                                                                                   <  1:4022781/101‑1 (MQ=255)
                    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGa                                                                                   <  1:1083374/101‑1 (MQ=255)
                    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGa                                                                                   >  2:3549874/1‑101 (MQ=255)
                    tCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGa                                                                                   <  1:548626/101‑1 (MQ=255)
                      gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.Gc                                                                                >  1:2611421/1‑101 (MQ=255)
                       aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.Gcc                                                                               >  1:3930954/1‑101 (MQ=255)
                        aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCg                                                                              >  1:2944103/1‑101 (MQ=255)
                          ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.Ca                                                                           >  2:3412179/1‑101 (MQ=255)
                            accTTGCGGATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGt                                                                         <  2:440112/101‑1 (MQ=255)
                             ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTc                                                                        >  1:479137/1‑101 (MQ=255)
                             ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTc                                                                        <  1:1691289/101‑1 (MQ=255)
                                ggCGGATTGCCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCt                                                                     <  2:3358922/100‑1 (MQ=255)
                                  cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTcc                                                                   >  1:1327968/1‑101 (MQ=255)
                                    gATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTcccc                                                                 <  1:664524/101‑1 (MQ=255)
                                    gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTcccc                                                                 <  2:174846/101‑1 (MQ=255)
                                    gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTcccc                                                                 <  1:1557197/101‑1 (MQ=255)
                                     aTTGTCACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTccccc                                                                <  2:3389022/101‑1 (MQ=255)
                                      ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTcccccc                                                               >  2:1159018/1‑101 (MQ=255)
                                       tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCcg                                                              >  2:2614954/1‑101 (MQ=255)
                                             tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCggggg                                                        >  2:429380/1‑101 (MQ=255)
                                             tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCggggg                                                        >  1:2091141/1‑101 (MQ=255)
                                              ttCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGc                                                       <  2:48514/101‑1 (MQ=255)
                                                  ggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG.CG..GCGAACACCGATACCG.C.GTCGCTCCCCCCGCGGGGGCGTg                                                    >  2:3429367/10‑101 (MQ=2)
                                                  ggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG.CG..GCGAACACCGATACCG.C.GTCGCTCCCCCCGCGGGGGCGTg                                                    >  1:2920285/10‑101 (MQ=2)
                                                   ggggCGAAGCGTCGCTCCCCCGGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGa                                                  >  1:3114965/1‑101 (MQ=255)
                                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGa                                                  >  2:3654634/1‑101 (MQ=255)
                                                   ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGa                                                  >  1:2549836/1‑101 (MQ=255)
                                                   gatcccgccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCCTGACCGCAT.CGA.A.GCTGGTCAC.ATCG.CCGTCGCTCCCCCCGCGGGGGCGTgg                                                   <  2:2197809/91‑1 (MQ=2)
                                                     gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG.CG..GCGAACACCGATACCG.C.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                 <  2:2750288/95‑1 (MQ=2)
                                                         gtgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAG.TGTCA.GGTTTTCCA.GCGG.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                             <  1:3235548/99‑1 (MQ=2)
                                                         ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG.CG..GCGAACACCGATACCG.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                             <  1:4104645/99‑1 (MQ=2)
                                                         catcGTCGCTCCCCCCGCGGGGACGTGGATCGAATCCTATGAGGGCTGGCGGCCC..G.GTCGTTGAT.CCCGCCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              <  2:3165394/98‑1 (MQ=0)
                                                         aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGaac                                            >  2:1892291/1‑100 (MQ=255)
                                                          tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG.CG..GCGAACACCGATACCG.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                            >  1:2362919/2‑101 (MQ=2)
                                                           gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                          >  1:2687011/1‑101 (MQ=255)
                                                            cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGG.CG..GCGAACACCGATACCG.C.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaaa                                          >  2:3298710/1‑99 (MQ=2)
                                                            cGTCGCCCCCCCCGCGGGGGCGGGTATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTCCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                         <  2:2250636/101‑1 (MQ=37)
                                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACTgga                                       >  1:1031621/1‑100 (MQ=255)
                                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                       <  1:2778662/101‑1 (MQ=255)
                                                              ccGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                       <  1:1816451/100‑1 (MQ=255)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTa                                      >  1:3056867/1‑101 (MQ=255)
                                                                 cttcccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGGAACTGGTATg                                    >  2:875290/4‑101 (MQ=255)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                   >  1:2966343/1‑101 (MQ=255)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                   >  2:1603454/1‑101 (MQ=255)
                                                                  cccccccGGGGGGGGGGGGATCAAAACACCTGCTTGTAGTGGG.CG..G.GCGGTCGGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                   <  1:3873470/100‑1 (MQ=21)
                                                                      cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                               <  2:3153209/101‑1 (MQ=255)
                                                                      cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                               <  1:4199048/101‑1 (MQ=255)
                                                                         ggggggcCGTGAATCGAAACATCTGGTTGTGGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                            <  2:1048267/99‑1 (MQ=37)
                                                                         gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGGTCTGgc                            >  1:594519/1‑101 (MQ=255)
                                                                         gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                            >  1:1403156/1‑101 (MQ=255)
                                                                         gcgGGGGCGTGGACCGAAACATCGGCGGGGAGTCGG.CG..G.GCAGTTGGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                            <  2:4017738/101‑1 (MQ=37)
                                                                           gggggCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCT.CAGTCGATTCCCCCGCGGGGGCGTGGATCGAAACTGTTATGGATGTgactc                          >  2:2072258/1‑97 (MQ=37)
                                                                           gggggCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                          >  2:3200016/1‑101 (MQ=255)
                                                                             gggCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaa                        >  2:2959342/1‑101 (MQ=255)
                                                                               gCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                      >  1:2737661/1‑101 (MQ=255)
                                                                                cGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                     >  2:3742661/1‑101 (MQ=255)
                                                                                cGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                     >  1:1783396/1‑101 (MQ=255)
                                                                                 gTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                    <  1:3185690/101‑1 (MQ=255)
                                                                                 gTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                    >  1:640301/1‑101 (MQ=255)
                                                                                   gggtCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                  >  1:2416128/4‑101 (MQ=255)
                                                                                   ggATCGAAACATCTGCAGGTGGTCGG.CG..G.GCAGTGCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                  <  2:1818926/101‑1 (MQ=255)
                                                                                   ggATCGAAACAGCTGCATGTAGGCGG.GG..G.CCAGTGCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                  <  2:2476630/101‑1 (MQ=37)
                                                                                    gATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                 >  2:616699/1‑101 (MQ=255)
                                                                                    gATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                 >  2:1187649/1‑101 (MQ=255)
                                                                                    gATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAATGGTATCGGTCTGGCGCAAACGGaag                 >  2:3092418/1‑100 (MQ=255)
                                                                                       cGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc              >  1:964882/1‑101 (MQ=255)
                                                                                       cGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc              >  2:554227/1‑101 (MQ=255)
                                                                                        gAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt             <  2:2609382/101‑1 (MQ=255)
                                                                                        gAAACACCGGCAGGTAGGCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt             <  2:3901031/101‑1 (MQ=37)
                                                                                             aTCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt        <  2:591493/101‑1 (MQ=255)
                                                                                              tCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc       >  2:1643354/1‑101 (MQ=255)
                                                                                              tCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc       >  2:3507926/1‑101 (MQ=255)
                                                                                              tCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTggcgtc       >  1:919872/1‑101 (MQ=255)
                                                                                                 gCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt    >  2:1465308/1‑101 (MQ=255)
                                                                                                 gCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt    >  2:1925088/1‑101 (MQ=255)
                                                                                                 gCATGTAGCCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt    <  2:1876006/101‑1 (MQ=255)
                                                                                                  cATGGAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTc   <  1:4114201/101‑1 (MQ=255)
                                                                                                   gtGTGGTCGG.GG..G.GGAGTGGGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcc  <  1:2382559/100‑1 (MQ=37)
                                                                                                   aTGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcc  <  2:2686805/101‑1 (MQ=255)
                                                                                                   aTGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcc  >  2:1290095/1‑101 (MQ=255)
                                                                                                    |                                                                                                         
TCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGG.CG..G.GCAGTTCGA.GCCG.CAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCC  >  NC_014034/1395572‑1395771

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG