Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,556 C→A 9.2% intergenic (+368/‑563) 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,5560C→A 9.2% 69.7 ‑0.2 ‑0.0358intergenic (+368/‑563)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (16/17):  ref base (152/173):  total (168/190)
Fisher's exact test strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.22e-01

Only 100 of 380 total aligned reads displayed.

GGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395456‑1395656
                                                                                                     |                                                                                                     
 gATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGcc                                                                                                      >  1:4532266/1‑101 (MQ=255)
 gATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGcc                                                                                                      >  1:392516/1‑101 (MQ=255)
 gATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGcc                                                                                                      >  2:2321401/1‑101 (MQ=255)
    cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.c                                                                                                  >  2:4271389/1‑101 (MQ=255)
    cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.c                                                                                                  >  1:138751/1‑101 (MQ=255)
          aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtc                                                                                            >  1:3710328/1‑101 (MQ=255)
           tcgcgtcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcg                                                                                           <  1:4171373/95‑1 (MQ=255)
           tCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcg                                                                                           <  1:2158429/101‑1 (MQ=255)
             gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgtc                                                                                         >  2:1276227/1‑101 (MQ=255)
              ccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgtcg                                                                                        <  1:3479559/101‑1 (MQ=255)
                 cTTTGGCGGTGGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTc                                                                                     <  2:1430966/101‑1 (MQ=255)
                 cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTc                                                                                     <  2:4649496/101‑1 (MQ=255)
                    tGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccc                                                                                  >  1:335522/1‑101 (MQ=255)
                       cGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcg                                                                               >  2:2554847/1‑101 (MQ=255)
                        ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcgc                                                                              >  1:1826125/1‑101 (MQ=255)
                         gTCGACAGGCTGTTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAATCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcgcg                                                                             <  1:1118426/101‑1 (MQ=255)
                             aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCTCggggg                                                                         >  1:2371916/1‑101 (MQ=255)
                                ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGt                                                                      <  1:2910749/101‑1 (MQ=255)
                                ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGt                                                                      <  1:551179/101‑1 (MQ=255)
                                ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGt                                                                      <  2:3903369/101‑1 (MQ=255)
                                  gccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  1:3284377/91‑1 (MQ=0)
                                  cTGGTAACGTGCCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  2:4084897/101‑1 (MQ=255)
                                  agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTgg                                                                    <  1:2688841/91‑1 (MQ=0)
                                   tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   >  2:3412705/1‑101 (MQ=255)
                                   tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGa                                                                   <  1:2088977/101‑1 (MQ=255)
                                    cctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGag                                                                  >  2:386029/9‑100 (MQ=2)
                                     ggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGGGGATc                                                                 >  2:3217497/8‑101 (MQ=0)
                                     ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 <  2:672123/94‑1 (MQ=2)
                                     ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCCCGGGGGCGTGGGTc                                                                 >  1:4022789/8‑101 (MQ=2)
                                     ctcgcacgtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                 <  2:4538886/94‑1 (MQ=2)
                                      tcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                <  1:1877516/95‑1 (MQ=2)
                                      tcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                <  1:813965/95‑1 (MQ=2)
                                      tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                <  1:1276318/101‑1 (MQ=255)
                                       cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCCGGGCCGTGGATCGa                                                               >  2:1383468/6‑101 (MQ=2)
                                       cgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCTGGGATCGa                                                               >  2:3592495/6‑101 (MQ=2)
                                       cgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                               >  1:3631168/6‑101 (MQ=2)
                                        ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              <  2:2992313/97‑1 (MQ=2)
                                        ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                              >  2:2245597/5‑101 (MQ=2)
                                        gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              <  2:3131808/97‑1 (MQ=2)
                                        gcaagtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              <  1:625392/97‑1 (MQ=2)
                                        aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                              <  2:1836224/101‑1 (MQ=255)
                                         gtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                             <  1:3186936/98‑1 (MQ=2)
                                          tggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                            >  1:3543498/3‑101 (MQ=2)
                                           tgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                                           >  2:2840177/1‑101 (MQ=255)
                                           ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                                           <  2:2728273/100‑2 (MQ=2)
                                           ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTACCCCCGCGGGGGCGTGGATCGAAact                                                           >  2:632360/2‑100 (MQ=2)
                                            gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacca                                                          <  1:2474548/101‑1 (MQ=255)
                                            gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTACTGGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAaccc                                                          >  1:4057870/1‑100 (MQ=37)
                                             tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccac                                                         >  1:259179/1‑101 (MQ=255)
                                             tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                                         >  1:748365/1‑98 (MQ=2)
                                             tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                                         >  2:706388/1‑98 (MQ=2)
                                              cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                                        <  1:4343384/101‑1 (MQ=255)
                                              cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                                        >  1:2152533/1‑98 (MQ=2)
                                              cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                                        <  1:1972686/101‑4 (MQ=2)
                                              cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                                        >  1:3459559/1‑97 (MQ=2)
                                               gccccccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCt                                                       <  2:3060836/98‑1 (MQ=255)
                                               gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgggg                                                       >  2:474478/1‑96 (MQ=2)
                                                cTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgct                                                      <  1:4457710/101‑7 (MQ=2)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                                     >  1:2753420/1‑101 (MQ=255)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                                     >  2:849138/1‑101 (MQ=255)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcgg                                                     >  1:4426174/1‑95 (MQ=2)
                                                 cccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcgg                                                     <  2:749696/100‑7 (MQ=2)
                                                  ccccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGTTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGc                                                    >  2:4610531/1‑101 (MQ=255)
                                                   cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                                                   <  2:1724250/101‑10 (MQ=2)
                                                   cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                                                   >  1:1926160/1‑92 (MQ=2)
                                                   cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                                                   <  2:2455227/101‑10 (MQ=2)
                                                    ccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCgg                                                  <  2:335488/101‑1 (MQ=255)
                                                    ccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggcgt                                                  <  1:932470/101‑10 (MQ=2)
                                                        gcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTg                                              >  1:1020859/1‑101 (MQ=255)
                                                          gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                                            <  1:4271636/101‑1 (MQ=255)
                                                            gggCGTGGATCGAAACATCTCGGTGCCGGCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                                          <  1:3811910/101‑1 (MQ=255)
                                                             ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                                         <  2:4369007/101‑1 (MQ=255)
                                                              gCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                                        >  2:4402939/1‑101 (MQ=255)
                                                               ggTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                       <  2:1940258/100‑1 (MQ=255)
                                                                gTGGATCGAAACATCGCGGTGCCGTCGCTCG.TCCAGCCGG.CCCGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                                      <  2:1896781/101‑1 (MQ=255)
                                                                 tGGCTCGAAACATCTCCGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                                     <  2:3767642/101‑1 (MQ=255)
                                                                  ggATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgc                                    >  1:1527343/1‑101 (MQ=255)
                                                                   gATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGGTCCTTTcgcg                                   >  2:163983/1‑101 (MQ=255)
                                                                      cGAAACATCTCGGGGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTTTCACTTTCGCgggg                                >  1:713683/1‑101 (MQ=255)
                                                                        aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg                              >  1:2309426/1‑101 (MQ=255)
                                                                         aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa                             >  2:2168050/1‑101 (MQ=255)
                                                                                cGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTcc                      >  2:3724563/1‑100 (MQ=255)
                                                                                cGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCg                      <  1:3412911/101‑1 (MQ=255)
                                                                                cGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCg                      >  2:3921089/1‑101 (MQ=255)
                                                                                  gTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCt                    <  2:902050/101‑1 (MQ=255)
                                                                                  gTGCCGCCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCCTCCGGATTGTCACTTTTGGGGGGCGAAGCGACgac                    >  2:3292262/1‑99 (MQ=37)
                                                                                       gTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTccccc               <  2:2032542/101‑1 (MQ=255)
                                                                                        tCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcccccc              <  2:608126/101‑1 (MQ=255)
                                                                                        tCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcccccc              <  1:241589/101‑1 (MQ=255)
                                                                                         cGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCcg             >  2:2400110/1‑101 (MQ=255)
                                                                                         cGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCcg             >  2:1245631/1‑101 (MQ=255)
                                                                                           cTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCcgcg           >  2:2494997/1‑101 (MQ=255)
                                                                                           cTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGACGCTCCCCCcgcg           >  1:4079835/1‑101 (MQ=255)
                                                                                              g.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCgggg        <  1:472868/101‑1 (MQ=255)
                                                                                               ggccAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCgggg        <  2:1020780/99‑1 (MQ=255)
                                                                                                tCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCggggg       >  2:1102422/1‑101 (MQ=255)
                                                                                                tCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCggggg       >  1:2530982/1‑101 (MQ=255)
                                                                                                  cAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCg     <  1:306283/101‑1 (MQ=255)
                                                                                                   aGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGt    <  2:1469653/101‑1 (MQ=255)
                                                                                                     ccGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTgg  <  1:292267/101‑1 (MQ=255)
                                                                                                     |                                                                                                     
GGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395456‑1395656

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG