Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,816 G→C 4.5% intergenic (+628/‑303) 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,8160G→C 4.5% 22.7 ‑1.3 ‑0.9290intergenic (+628/‑303)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (153/123):  total (158/132)
Fisher's exact test strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.15e-01

Only 100 of 304 total aligned reads displayed.

GGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGG  >  NC_014034/1395716‑1395916
                                                                                                     |                                                                                                      
ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                        <  1:2922227/101‑1 (MQ=255)
ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                        <  1:5346330/101‑1 (MQ=255)
ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                        >  2:1086876/1‑101 (MQ=255)
 gCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGa                                                                                                       <  1:1361793/101‑1 (MQ=255)
 gCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGa                                                                                                       >  2:3422218/1‑101 (MQ=255)
   gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                                     >  2:5705846/1‑101 (MQ=255)
    tgcaTCGAACCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                                                                                    <  2:624051/98‑1 (MQ=255)
     ggATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGa                                                                                                   >  1:1642303/1‑101 (MQ=255)
     ggATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGa                                                                                                   >  2:2936542/1‑101 (MQ=255)
       aTCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.Cg                                                                                                >  2:5748825/1‑101 (MQ=255)
       aTCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.Cg                                                                                                <  2:3066557/101‑1 (MQ=255)
          gAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCg                                                                                             >  2:1715996/1‑101 (MQ=255)
          gAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCg                                                                                             >  2:3267762/1‑101 (MQ=255)
            aaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTg                                                                                           >  1:1107059/1‑101 (MQ=255)
            aaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTg                                                                                           <  2:4070834/101‑1 (MQ=255)
              cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                         <  2:2302331/101‑1 (MQ=255)
              cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                         <  1:4146816/101‑1 (MQ=255)
              cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                         <  1:560281/101‑1 (MQ=255)
               tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                       >  1:4858102/1‑101 (MQ=255)
               tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                       <  2:257769/101‑1 (MQ=255)
                ggTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.Cg                                                                                      <  2:2043389/101‑1 (MQ=255)
                  tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                    >  1:3405087/1‑101 (MQ=255)
                  tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                    >  2:2439240/1‑101 (MQ=255)
                     ggATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCt                                                                                 <  2:3295916/101‑1 (MQ=255)
                        tCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccc                                                                              <  1:3793300/101‑1 (MQ=255)
                         cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                             >  1:4361293/1‑101 (MQ=255)
                         cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                             <  2:5274994/101‑1 (MQ=255)
                         cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGATGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                             >  1:4939063/1‑101 (MQ=255)
                          tGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                                                                            >  1:2205142/1‑101 (MQ=255)
                           ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                           <  2:236541/101‑1 (MQ=255)
                           cgatctaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                           <  1:657254/95‑1 (MQ=255)
                             cgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgc                                                                         >  2:3810500/1‑101 (MQ=255)
                             cgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgc                                                                         >  2:839918/1‑101 (MQ=255)
                              gcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcccg                                                                        >  1:3373228/1‑101 (MQ=255)
                                aaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                                                                      >  2:3722761/1‑101 (MQ=255)
                                 aaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                                                                     <  2:3139127/101‑1 (MQ=255)
                                   cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                                   <  2:5070009/101‑1 (MQ=255)
                                      ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                                 <  1:2484643/92‑1 (MQ=0)
                                      ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                                 <  2:4506569/92‑1 (MQ=0)
                                      ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                                 >  1:1554772/10‑101 (MQ=0)
                                      aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                                >  1:1023535/1‑101 (MQ=255)
                                      aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                                <  1:4804818/101‑1 (MQ=255)
                                       aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                               >  1:3121786/1‑101 (MQ=255)
                                       aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                               <  1:5762172/101‑1 (MQ=255)
                                         gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGGGGAt                                                             >  2:5530636/1‑101 (MQ=255)
                                           cTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           <  1:1243645/101‑1 (MQ=255)
                                                gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACt                                                      >  2:1520897/1‑101 (MQ=255)
                                                gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                      <  1:4549792/101‑1 (MQ=255)
                                                  cgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgc                                                    <  1:4695201/101‑1 (MQ=255)
                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTgggc                                                  >  1:5297946/1‑101 (MQ=255)
                                                     tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgcc                                                 >  2:3192613/1‑101 (MQ=255)
                                                     tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgcc                                                 >  1:4303772/1‑101 (MQ=255)
                                                        ccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCgccgcc                                              >  1:2342065/1‑101 (MQ=255)
                                                          ccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTa                                            <  2:196753/101‑1 (MQ=255)
                                                           ctcgGGGTCGTGGATCGAAACCTGCG.GCGTGACTTGCCAGTGATAGC.CGCCGTGGC.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                           <  2:1372979/99‑1 (MQ=37)
                                                            gcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATc                                          <  1:1972435/101‑1 (MQ=255)
                                                             cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAtag                                         >  1:1373067/1‑99 (MQ=255)
                                                             cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATcc                                         >  2:4753159/1‑100 (MQ=255)
                                                             cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCg                                         <  2:3034676/101‑1 (MQ=255)
                                                             cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCGCCGCCTATCg                                         >  1:4041182/1‑101 (MQ=255)
                                                             cgGGGCCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGCGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAtgg                                         >  1:3892749/1‑99 (MQ=255)
                                                              gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGc                                        >  2:2888802/1‑101 (MQ=255)
                                                               ggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCa                                       >  2:5189379/1‑101 (MQ=255)
                                                                gggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCAt                                      >  1:5426897/1‑101 (MQ=255)
                                                                 ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                     <  2:2205020/101‑1 (MQ=255)
                                                                 ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                     >  1:1721435/1‑101 (MQ=255)
                                                                  gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                    <  1:2686957/101‑1 (MQ=255)
                                                                  gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                    >  1:1075071/1‑101 (MQ=255)
                                                                   cGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                                   <  1:1778487/101‑1 (MQ=255)
                                                                    gTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCgg                                  <  2:4926546/101‑1 (MQ=255)
                                                                     tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCggc                                 >  2:3599438/1‑100 (MQ=255)
                                                                      ggATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTc                                >  2:4144951/1‑101 (MQ=255)
                                                                       gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                               <  1:5560128/101‑1 (MQ=255)
                                                                       gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                               >  2:4623860/1‑101 (MQ=255)
                                                                        aTCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGa                              >  2:3107585/1‑101 (MQ=255)
                                                                          cGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgc                            >  2:4270728/1‑101 (MQ=255)
                                                                            aaaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgcgc                          >  2:3892335/1‑101 (MQ=255)
                                                                             aaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgcc                         <  1:3267962/101‑1 (MQ=255)
                                                                               ccTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgccgc                       <  2:2111163/101‑1 (MQ=255)
                                                                               ccTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgccgc                       <  2:5401096/101‑1 (MQ=255)
                                                                                cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCcgcg                      <  2:4361027/101‑1 (MQ=255)
                                                                                cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCcgcg                      >  1:445575/1‑101 (MQ=255)
                                                                                   cg.gcgTGAGTTGCCAGTGATAGT.CGCCGTGGA.CGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAg                   <  2:4865900/101‑1 (MQ=255)
                                                                                    g.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCCGt                  <  1:3299885/101‑1 (MQ=255)
                                                                                      gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGATCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTc                 <  1:695049/101‑1 (MQ=255)
                                                                                       ggTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg                <  1:4922327/100‑1 (MQ=255)
                                                                                       cgTGATTTGCCAGTTATAGC.CGCCGTGGA.CGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg                <  1:3957071/101‑1 (MQ=38)
                                                                                       cgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg                <  1:2248935/101‑1 (MQ=255)
                                                                                         tGAGATGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCCGCCGcc              >  2:1696046/1‑100 (MQ=255)
                                                                                          gAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTc             >  1:3000615/1‑101 (MQ=255)
                                                                                          gAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTc             <  1:3192802/101‑1 (MQ=255)
                                                                                           aGTTGCCAGTGATAGG.CCCCGGGGA.CCTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcc            <  1:4077183/101‑1 (MQ=38)
                                                                                            gTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccc           <  1:3601286/101‑1 (MQ=255)
                                                                                            gTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccc           <  2:1032943/101‑1 (MQ=255)
                                                                                              tGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccccc         >  2:3692734/1‑101 (MQ=255)
                                                                                              tGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccccc         >  2:5582521/1‑101 (MQ=255)
                                                                                               gCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCCCTcccccc        >  2:4971940/1‑101 (MQ=255)
                                                                                                ccAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcg       >  1:772102/1‑101 (MQ=255)
                                                                                                   gTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCgg    >  2:781711/1‑101 (MQ=255)
                                                                                                     gATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCgggg  <  1:839986/101‑1 (MQ=255)
                                                                                                     |                                                                                                      
GGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGGGG  >  NC_014034/1395716‑1395916

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