Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,814 G→T 4.4% intergenic (+626/‑305) 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,8140G→T 4.4% 21.6 ‑0.9 ‑0.1297intergenic (+626/‑305)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (152/132):  total (157/140)
Fisher's exact test strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.15e-01

Only 100 of 305 total aligned reads displayed.

GGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGG  >  NC_014034/1395714‑1395914
                                                                                                     |                                                                                                      
 gggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGt                                                                                                       <  2:1056214/101‑1 (MQ=255)
  ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                      <  2:4592524/101‑1 (MQ=255)
  ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                      <  1:5346330/101‑1 (MQ=255)
   gCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGa                                                                                                     <  1:1361793/101‑1 (MQ=255)
     gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                                   <  2:3951481/101‑1 (MQ=255)
     gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                                   >  1:493636/1‑101 (MQ=255)
     gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                                   >  2:5705846/1‑101 (MQ=255)
      tgcaTCGAACCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                                                                                  <  2:624051/98‑1 (MQ=255)
      tGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                                                                                  >  2:3183469/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)
             aaaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGt                                                                                          >  1:4071257/1‑101 (MQ=255)
              aaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTg                                                                                         >  1:1107059/1‑101 (MQ=255)
                cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                       <  1:560281/101‑1 (MQ=255)
                cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                       <  2:2302331/101‑1 (MQ=255)
                 tccgatctaTCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                     <  1:2326807/93‑1 (MQ=255)
                 tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                     <  2:257769/101‑1 (MQ=255)
                  ggTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.Cg                                                                                    <  2:536454/101‑1 (MQ=255)
                    tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                  >  2:3242263/1‑101 (MQ=255)
                    tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                  <  2:1589281/101‑1 (MQ=255)
                    tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                  >  1:3405087/1‑101 (MQ=255)
                    tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                  >  2:2439240/1‑101 (MQ=255)
                     aTGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCg                                                                                 <  2:3637207/101‑1 (MQ=255)
                      tGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGc                                                                                <  2:2100366/101‑1 (MQ=255)
                       ggATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGaagatcgg                                                                                >  2:450290/1‑94 (MQ=255)
                           cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                           <  1:1736450/101‑1 (MQ=255)
                            tGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                                                                          >  2:3693871/1‑101 (MQ=255)
                             ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                         <  2:236541/101‑1 (MQ=255)
                             ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                         <  1:467392/101‑1 (MQ=255)
                             cgatctaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                         <  1:657254/95‑1 (MQ=255)
                                gcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgcg                                                                      <  2:3751965/101‑1 (MQ=255)
                                gcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcccg                                                                      >  1:3373228/1‑101 (MQ=255)
                                 cAAAGGGAAAGTCTGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                     <  2:5017139/101‑1 (MQ=255)
                                 cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                     <  1:2653680/101‑1 (MQ=255)
                                 cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                     <  1:2400696/101‑1 (MQ=255)
                                 cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                     >  2:3450623/1‑101 (MQ=255)
                                  aaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                                                                    >  1:4420985/1‑101 (MQ=255)
                                   aaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                                                                   <  2:3139127/101‑1 (MQ=255)
                                      ggAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg                                                                <  1:1716095/101‑1 (MQ=255)
                                       gAAAGTCTGGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGt                                                               <  1:4017808/101‑1 (MQ=255)
                                       cggccccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCg                                                                <  1:1637411/91‑1 (MQ=0)
                                        ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                               >  2:5267004/10‑101 (MQ=0)
                                        aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                              <  2:3368930/101‑1 (MQ=255)
                                        aaaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTg                                                              >  1:1023535/1‑101 (MQ=255)
                                         gacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGTg                                                              <  1:2996774/93‑1 (MQ=0)
                                         aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                             <  1:5762172/101‑1 (MQ=255)
                                          aGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGGGGa                                                            >  2:4958963/1‑101 (MQ=255)
                                           gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:2176042/101‑1 (MQ=255)
                                           gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGGGGAt                                                           >  2:5530636/1‑101 (MQ=255)
                                             cTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  1:1243645/101‑1 (MQ=255)
                                             cTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGGTGCCAGTGATAGA.CGCCGTGGC.CGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         >  2:1814693/1‑101 (MQ=255)
                                                tcgtcgCCCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      <  2:5188700/101‑1 (MQ=255)
                                                 cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                     >  1:2243835/1‑101 (MQ=255)
                                                  gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                                    <  1:3183658/101‑1 (MQ=255)
                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                                 >  1:3920659/1‑101 (MQ=255)
                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                                 >  2:4702379/1‑101 (MQ=255)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcgc                                                >  1:4206284/1‑101 (MQ=255)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcgc                                                >  2:4098462/1‑101 (MQ=255)
                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTgggc                                                >  1:5297946/1‑101 (MQ=255)
                                                        ccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccg                                              <  2:4784410/101‑1 (MQ=255)
                                                        ccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccg                                              >  2:2013837/1‑101 (MQ=255)
                                                        ccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccg                                              <  2:2165857/101‑1 (MQ=255)
                                                         cccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccgc                                             >  1:4324690/1‑101 (MQ=255)
                                                           cccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCt                                           <  2:4258656/101‑1 (MQ=255)
                                                              gcgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATc                                        <  2:1478883/101‑1 (MQ=255)
                                                               cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAtag                                       >  1:1373067/1‑99 (MQ=255)
                                                               cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCg                                       >  2:704564/1‑101 (MQ=255)
                                                               cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCg                                       <  2:3034676/101‑1 (MQ=255)
                                                                gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGc                                      >  2:2888802/1‑101 (MQ=255)
                                                                 ggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGGCCGCCTATCGcc                                     >  1:1545620/1‑100 (MQ=255)
                                                                  gggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCAt                                    >  1:5426897/1‑101 (MQ=255)
                                                                   ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                   >  1:3079630/1‑101 (MQ=255)
                                                                    gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                  <  1:3783248/101‑1 (MQ=255)
                                                                    gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                  >  1:4845161/1‑101 (MQ=255)
                                                                    gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                  >  1:1075071/1‑101 (MQ=255)
                                                                     cGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                                 >  1:4750625/1‑101 (MQ=255)
                                                                       tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGt                               >  1:3826068/1‑101 (MQ=255)
                                                                       tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGt                               >  1:4145379/1‑101 (MQ=255)
                                                                         gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                             <  1:195161/101‑1 (MQ=255)
                                                                         gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                             <  1:882525/101‑1 (MQ=255)
                                                                         gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                             <  1:5560128/101‑1 (MQ=255)
                                                                          aTCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGa                            >  2:3107585/1‑101 (MQ=255)
                                                                          aTCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGa                            <  2:1710515/101‑1 (MQ=255)
                                                                            cGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgc                          >  2:83980/1‑101 (MQ=255)
                                                                            cGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgc                          >  1:2839054/1‑101 (MQ=255)
                                                                              aaaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCCCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgcgc                        >  2:3816123/1‑101 (MQ=255)
                                                                               aaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgcc                       <  1:3267962/101‑1 (MQ=255)
                                                                                  cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCTAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCcgcg                    >  1:4886796/1‑101 (MQ=255)
                                                                                  cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCcgcg                    <  2:4361027/101‑1 (MQ=255)
                                                                                    gcg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCa                  <  1:4530041/101‑1 (MQ=255)
                                                                                      g.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCCGt                <  1:3299885/101‑1 (MQ=255)
                                                                                        gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGATCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTc               <  1:695049/101‑1 (MQ=255)
                                                                                         cgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg              <  2:4886515/101‑1 (MQ=255)
                                                                                            gAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTc           >  2:3292887/1‑101 (MQ=255)
                                                                                             aGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcc          <  2:1319972/101‑1 (MQ=255)
                                                                                              gTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccc         <  2:1032943/101‑1 (MQ=255)
                                                                                              gTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccc         <  1:3601286/101‑1 (MQ=255)
                                                                                                tGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccccc       >  2:5582521/1‑101 (MQ=255)
                                                                                                  ccAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcg     >  1:772102/1‑101 (MQ=255)
                                                                                                     |                                                                                                      
GGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGCGG  >  NC_014034/1395714‑1395914

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