Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,675 A→T 8.9% intergenic (+487/‑444) 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,6750A→T 8.9% 52.0 ‑1.8 ‑0.6237intergenic (+487/‑444)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/13):  ref base (120/94):  total (129/108)
Fisher's exact test strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.38e-01

Only 100 of 263 total aligned reads displayed.

CCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCC  >  NC_014034/1395575‑1395775
                                                                                                    |                                                                                                       
 cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                        <  1:4438630/101‑1 (MQ=255)
  ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                       >  1:4588781/1‑101 (MQ=255)
     cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCGACATGTGGTcgg                                                                                                    >  2:642142/1‑101 (MQ=38)
      gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                                   >  2:2679407/1‑101 (MQ=255)
        gggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcgg                                                                                                 >  1:3949122/1‑101 (MQ=255)
            gTGGATCGAAGCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAg                                                                                           >  1:4364701/1‑101 (MQ=255)
            gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGGAGTCGGCGG..GCAg                                                                                           >  2:2115600/1‑101 (MQ=255)
             tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                          >  1:970804/1‑101 (MQ=255)
             tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                          >  2:2087524/1‑101 (MQ=255)
             tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGt                                                                                          <  2:3924753/101‑1 (MQ=255)
             tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAgg                                                                                          <  2:2370631/101‑2 (MQ=255)
              ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGtt                                                                                         <  2:1509675/101‑1 (MQ=255)
              ggATCGAAAACACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGtt                                                                                         >  2:3704904/1‑101 (MQ=255)
                   gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gc                                                                                   >  2:3124007/1‑101 (MQ=255)
                    aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.Gcc                                                                                  <  2:725296/101‑1 (MQ=255)
                     aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCg                                                                                 <  1:2624974/101‑1 (MQ=255)
                      accaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGc                                                                                >  2:4399477/1‑101 (MQ=255)
                          ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTc                                                                            <  2:5004547/101‑1 (MQ=255)
                          ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTc                                                                            <  1:3118157/101‑1 (MQ=255)
                              gCGGATTGTCTCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                        <  1:2165062/101‑1 (MQ=255)
                              gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                        >  1:2694637/1‑101 (MQ=255)
                              gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTc                                                                        >  2:4212090/1‑101 (MQ=255)
                               cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcc                                                                       <  1:2733478/101‑1 (MQ=255)
                                 gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccc                                                                     >  1:2536429/1‑101 (MQ=255)
                                 gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccc                                                                     >  2:1759941/1‑101 (MQ=255)
                                   ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                   >  2:4413192/1‑101 (MQ=255)
                                   ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                   >  2:837391/1‑101 (MQ=255)
                                    tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                  >  2:1144047/1‑101 (MQ=255)
                                    tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                  <  2:1076791/101‑1 (MQ=255)
                                    tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                  <  2:2587669/101‑1 (MQ=255)
                                     gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                 <  2:4026298/101‑1 (MQ=255)
                                     gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                 >  1:325436/1‑101 (MQ=255)
                                      tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                                >  1:4180695/1‑101 (MQ=255)
                                           ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                           >  1:1020281/1‑101 (MQ=255)
                                             cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                         >  1:4545643/1‑101 (MQ=255)
                                              tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                         <  1:293838/91‑1 (MQ=2)
                                              gcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                        <  2:2207154/101‑1 (MQ=255)
                                               gggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                       <  2:2146594/100‑1 (MQ=255)
                                               gggggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                       <  1:1100674/100‑1 (MQ=255)
                                                gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                       <  1:2882074/93‑1 (MQ=2)
                                                ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                      <  2:2230151/101‑1 (MQ=255)
                                                 gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                     >  1:4273876/1‑101 (MQ=255)
                                                  ggggAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                    <  1:837493/98‑1 (MQ=255)
                                                  ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATTTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                    <  2:3336160/101‑1 (MQ=255)
                                                  ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                    >  1:3128158/1‑101 (MQ=255)
                                                  ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                    <  1:1498504/101‑1 (MQ=255)
                                                   gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                   >  1:3793015/1‑101 (MQ=255)
                                                   gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                   <  1:720017/101‑1 (MQ=255)
                                                    cGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                  <  1:3419938/101‑1 (MQ=255)
                                                      aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                <  1:2306476/101‑1 (MQ=255)
                                                       gtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.TACACCGGGGCTCTGGTG..ACGGTTTGG.GGCGGAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcg                                              >  2:3242960/4‑100 (MQ=2)
                                                       aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                               >  2:3435013/1‑101 (MQ=255)
                                                       aGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                               >  1:2233681/1‑101 (MQ=255)
                                                        gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt                                              <  1:1695066/101‑1 (MQ=255)
                                                        gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                               <  2:2927946/101‑1 (MQ=2)
                                                         cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                             >  2:4326350/1‑101 (MQ=255)
                                                           tcccTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                           <  1:3689995/98‑1 (MQ=255)
                                                           gcGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                           <  2:2724383/100‑1 (MQ=255)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCGGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcatcga                                          <  1:4308045/101‑6 (MQ=2)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGC.GAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGGGGGGGTCGAAAcacccatg                                        >  1:1714131/1‑94 (MQ=2)
                                                               cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                       <  2:1599850/100‑1 (MQ=255)
                                                                ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGCGGCAGGGTTTCCA.GCGGC.GTCCCTCCCCCCCCGGGGGCGGGGATCGAAAcccgccccg                                       >  1:2049403/1‑93 (MQ=0)
                                                                 cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                     <  1:2087672/101‑1 (MQ=255)
                                                                   cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                                   <  1:2805292/101‑1 (MQ=255)
                                                                     cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                                 >  2:2836856/1‑101 (MQ=255)
                                                                      gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                                >  1:1235156/1‑101 (MQ=255)
                                                                      gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                                <  1:3567444/101‑1 (MQ=255)
                                                                        gggggGCTGGACCGAAATATCTGCGTGTAGTCGGGGG..GCTGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                              <  2:3835088/101‑1 (MQ=25)
                                                                          gggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaa                            <  2:351360/101‑1 (MQ=255)
                                                                           ggggTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                           <  1:3068081/98‑1 (MQ=255)
                                                                           ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                           >  1:728828/1‑101 (MQ=255)
                                                                             cGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                         >  2:4930133/1‑101 (MQ=255)
                                                                              gTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                        >  1:1530876/1‑101 (MQ=255)
                                                                              gTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACTGGTATGGATCTGGCGCAAACgg                        >  2:479927/1‑101 (MQ=255)
                                                                               tGGCTCGAAACATCTGCTTGTAGGCGGCGG..GGAGTGGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAGCTGGCGCAAACGGa                       <  1:334782/101‑1 (MQ=25)
                                                                                ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                      >  1:4133036/1‑101 (MQ=255)
                                                                                ggATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                      >  2:170361/1‑101 (MQ=255)
                                                                                 ggtCGAAGCATCTGGTGGTGGGCGGCGG..GCGGTTGGA.GCCGGAGTGGCGCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                     <  1:3611064/99‑1 (MQ=21)
                                                                                 gATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                     >  1:138248/1‑101 (MQ=255)
                                                                                 gATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                     >  2:3603901/1‑101 (MQ=255)
                                                                                 gATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa                     >  2:3665253/1‑101 (MQ=255)
                                                                                  aTCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg                    <  2:1310250/101‑1 (MQ=255)
                                                                                   tCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                   >  1:2593528/1‑101 (MQ=255)
                                                                                   tCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt                   >  1:2426981/1‑101 (MQ=255)
                                                                                    cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                  >  2:4527218/1‑101 (MQ=255)
                                                                                    cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                  >  1:516782/1‑101 (MQ=255)
                                                                                    cGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc                  <  2:2886316/101‑1 (MQ=255)
                                                                                     gAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt                 >  2:1896613/1‑101 (MQ=255)
                                                                                     gAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt                 >  2:1940270/1‑101 (MQ=255)
                                                                                      aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                <  1:89249/101‑1 (MQ=255)
                                                                                      aaaCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg                <  1:220889/101‑1 (MQ=255)
                                                                                       aaCATCTGCATGTGGTGGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt               <  1:2106180/101‑1 (MQ=255)
                                                                                        aCATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc              >  1:1928921/1‑101 (MQ=255)
                                                                                         cATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg             >  2:4155464/1‑101 (MQ=255)
                                                                                         cATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg             <  1:4909450/101‑1 (MQ=255)
                                                                                                 tGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTccc     >  2:228840/1‑101 (MQ=255)
                                                                                                  gTAGTCGGGGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc    <  1:1161287/101‑1 (MQ=255)
                                                                                                  gTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc    <  2:4599327/101‑1 (MQ=255)
                                                                                                  gTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccc    >  1:4734729/1‑101 (MQ=255)
                                                                                                    aGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTcccccc  >  2:55441/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
CCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG..GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCC  >  NC_014034/1395575‑1395775

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