Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,806 A→C 4.5% intergenic (+618/‑313) 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,8060A→C 4.5% 25.8 ‑0.9 ‑1.4290intergenic (+618/‑313)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (148/129):  total (153/137)
Fisher's exact test strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.76e-02

Only 100 of 306 total aligned reads displayed.

CCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCC  >  NC_014034/1395706‑1395906
                                                                                                     |                                                                                                      
cccccGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGa                                                                                                        <  1:3181204/101‑1 (MQ=255)
cccccGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGa                                                                                                        <  2:1458257/101‑1 (MQ=255)
   ccGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGtt                                                                                                     >  2:882449/1‑101 (MQ=255)
      cgGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGcc                                                                                                  <  1:5664899/101‑1 (MQ=255)
      cgGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGcc                                                                                                  <  1:5083518/101‑1 (MQ=255)
       gggggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCa                                                                                                 >  1:5560344/1‑101 (MQ=255)
         gggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGt                                                                                               <  2:690836/101‑1 (MQ=255)
         gggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGt                                                                                               <  1:1518457/101‑1 (MQ=255)
         gggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGt                                                                                               <  2:1056214/101‑1 (MQ=255)
          ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                              <  1:704908/101‑1 (MQ=255)
             gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                           >  2:5705846/1‑101 (MQ=255)
             gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                           >  1:493636/1‑101 (MQ=255)
              tgcaTCGAACCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                                                                          <  2:624051/98‑1 (MQ=255)
              tGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                                                                          >  1:1157957/1‑101 (MQ=255)
              tGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAg                                                                                          >  2:3183469/1‑101 (MQ=255)
               ggATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGa                                                                                         >  1:1642303/1‑101 (MQ=255)
                gATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.c                                                                                       >  2:1972324/1‑101 (MQ=255)
                 aTCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.Cg                                                                                      >  2:5748825/1‑101 (MQ=255)
                  tCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGc                                                                                     <  1:3194997/101‑1 (MQ=255)
                      aaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTg                                                                                 >  1:3271256/1‑101 (MQ=255)
                        ccgatctGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                               <  1:450340/95‑1 (MQ=255)
                        cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                               >  2:3091211/1‑101 (MQ=255)
                        cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                               <  1:560281/101‑1 (MQ=255)
                        cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                               <  1:4146816/101‑1 (MQ=255)
                         tccgatctaTCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                             <  1:2326807/93‑1 (MQ=255)
                         tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGTCAGTGATAGA.CGCCGTGGA.c                                                                             <  2:2942923/101‑1 (MQ=255)
                         tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                             <  1:136234/101‑1 (MQ=255)
                         tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                             <  2:257769/101‑1 (MQ=255)
                          ggTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.Cg                                                                            <  1:5763400/101‑1 (MQ=255)
                            tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                          >  2:2439240/1‑101 (MQ=255)
                            tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                          >  2:3242263/1‑101 (MQ=255)
                             aTGGATCTGGCTCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCg                                                                         <  2:4717644/101‑1 (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)
                                 aTCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.cagatcgga                                                                     >  2:2326681/1‑93 (MQ=255)
                                  tCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccc                                                                    <  1:1258332/101‑1 (MQ=255)
                                  tCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccc                                                                    <  1:3793300/101‑1 (MQ=255)
                                   cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                   <  2:3701831/101‑1 (MQ=255)
                                   cTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                   <  1:1736450/101‑1 (MQ=255)
                                    tGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTccccc                                                                  >  2:3783005/1‑101 (MQ=255)
                                     ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                 >  1:339382/1‑101 (MQ=255)
                                     ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                 <  1:467392/101‑1 (MQ=255)
                                     ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                 <  1:798477/101‑1 (MQ=255)
                                         cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                             >  2:5309472/1‑101 (MQ=255)
                                         cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                             <  1:2400696/101‑1 (MQ=255)
                                         cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                             <  1:2653680/101‑1 (MQ=255)
                                         cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                             >  2:3450623/1‑101 (MQ=255)
                                          aaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                                                            >  2:3063034/1‑101 (MQ=255)
                                           aaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                                                           <  2:3139127/101‑1 (MQ=255)
                                            aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgcccggg                                                          <  2:4901552/101‑6 (MQ=255)
                                            aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggggg                                                          <  2:2246817/101‑1 (MQ=255)
                                              ggAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg                                                        <  1:1716095/101‑1 (MQ=255)
                                                ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                       >  1:1554772/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)
                                                  aGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGa                                                    <  1:4392092/101‑1 (MQ=255)
                                                   gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                   <  1:2176042/101‑1 (MQ=255)
                                                   gTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGAt                                                   <  1:664766/101‑1 (MQ=255)
                                                    tCTGTCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATc                                                  <  2:4972843/101‑1 (MQ=255)
                                                        tcgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                              <  2:1706943/101‑1 (MQ=255)
                                                         cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                             >  1:2243835/1‑101 (MQ=255)
                                                          gtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                            >  1:3368231/1‑101 (MQ=255)
                                                             gcgccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                         <  1:2767737/98‑1 (MQ=255)
                                                              cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcgc                                        >  2:4098462/1‑101 (MQ=255)
                                                              cTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACTgggc                                        >  1:5297946/1‑101 (MQ=255)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgcc                                       <  2:1564768/101‑1 (MQ=255)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgcc                                       >  1:4303772/1‑101 (MQ=255)
                                                                ccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccg                                      <  2:2165857/101‑1 (MQ=255)
                                                                   cccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCt                                   <  2:4258656/101‑1 (MQ=255)
                                                                    ccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTa                                  >  2:2376971/1‑101 (MQ=255)
                                                                       cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACTGCGCCGCCTATCg                               >  1:4041182/1‑101 (MQ=255)
                                                                       cgGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGCGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCg                               >  2:1719598/1‑101 (MQ=255)
                                                                       cgGGGCCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGCGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAtgg                               >  1:3892749/1‑99 (MQ=255)
                                                                        gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGc                              >  2:2888802/1‑101 (MQ=255)
                                                                           ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                           >  1:3079630/1‑101 (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:5431009/1‑101 (MQ=255)
                                                                             cGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                         >  2:4055227/1‑101 (MQ=255)
                                                                             cGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCg                         >  1:4750625/1‑101 (MQ=255)
                                                                              gTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCgg                        >  2:3844319/1‑101 (MQ=255)
                                                                               tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGt                       >  2:556470/1‑101 (MQ=255)
                                                                                 gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                     <  1:195161/101‑1 (MQ=255)
                                                                                 gATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCg                     >  2:4623860/1‑101 (MQ=255)
                                                                                   tCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGCNNAg                   >  1:3459976/1‑101 (MQ=255)
                                                                                    cGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgc                  >  2:83980/1‑101 (MQ=255)
                                                                                         ccTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgccgc             <  2:5401096/101‑1 (MQ=255)
                                                                                          cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCTAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCcgcg            >  1:4886796/1‑101 (MQ=255)
                                                                                          cTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCcgcg            <  2:4361027/101‑1 (MQ=255)
                                                                                             cg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAg         <  2:1023460/101‑1 (MQ=255)
                                                                                              g.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAgg        >  2:5252978/1‑100 (MQ=255)
                                                                                              g.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGt        >  2:5117503/1‑101 (MQ=255)
                                                                                                gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTc       <  2:4846176/101‑1 (MQ=255)
                                                                                                 cgTGATTTGCCAGTTATAGC.CGCCGTGGA.CGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg      <  1:3957071/101‑1 (MQ=38)
                                                                                                  gTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGc     >  1:53609/1‑101 (MQ=255)
                                                                                                    gAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTc   >  2:3292887/1‑101 (MQ=255)
                                                                                                     aGTTGCCAGTGATAGG.CCCCGGGGA.CCTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcc  <  1:4077183/101‑1 (MQ=38)
                                                                                                     aGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcc  <  2:1319972/101‑1 (MQ=255)
                                                                                                     |                                                                                                      
CCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCC  >  NC_014034/1395706‑1395906

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