Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,670 C→G 5.8% intergenic (+482/‑449) 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,6700C→G 5.8% 38.7 ‑0.2 ‑0.0383intergenic (+482/‑449)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/9):  ref base (196/163):  total (210/173)
Fisher's exact test strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.97e-01

Only 100 of 437 total aligned reads displayed.

GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTC  >  NC_014034/1395570‑1395770
                                                                                                    |                                                                                                       
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGGTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                         >  1:5260182/1‑101 (MQ=255)
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                         >  2:864779/1‑101 (MQ=255)
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                         <  1:4078446/101‑1 (MQ=255)
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                         <  1:217769/101‑1 (MQ=255)
 cTCCCCCTGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCa                                                                                                        <  1:2530316/101‑1 (MQ=255)
   ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                      >  2:5368054/1‑101 (MQ=255)
    cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCTCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                     <  1:3448363/101‑1 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                    >  1:714704/1‑101 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                    <  1:1436767/101‑1 (MQ=255)
      atctcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                   <  2:665478/97‑1 (MQ=255)
         gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                                <  2:2473618/101‑1 (MQ=255)
          cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg                                                                                               >  1:6448621/1‑101 (MQ=255)
            ggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcg                                                                                             >  1:2218683/1‑101 (MQ=255)
               gCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.Gc                                                                                        <  1:6128623/101‑1 (MQ=255)
               gCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.Gc                                                                                        <  2:4445359/101‑1 (MQ=255)
               gCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.Gc                                                                                        <  1:4191755/101‑1 (MQ=255)
                 gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAg                                                                                      >  1:5892832/1‑101 (MQ=255)
                 gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAg                                                                                      >  1:116810/1‑101 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGt                                                                                     <  2:3791909/101‑1 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGt                                                                                     >  1:3016505/1‑101 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGt                                                                                     <  1:6673696/101‑1 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGtt                                                                                    >  2:6074073/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGtt                                                                                    >  1:2617612/1‑101 (MQ=255)
                    gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTc                                                                                   <  1:1351721/101‑1 (MQ=255)
                    gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTc                                                                                   <  1:2429467/101‑1 (MQ=255)
                    gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTc                                                                                   <  2:1693857/101‑1 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCg                                                                                  >  2:5082713/1‑101 (MQ=255)
                        gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.Gc                                                                              >  2:2596562/1‑101 (MQ=255)
                         aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.gac                                                                             >  2:3192049/1‑99 (MQ=255)
                          aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCg                                                                            <  2:4743179/101‑1 (MQ=255)
                               ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTc                                                                       >  2:7080853/1‑101 (MQ=255)
                               ccTTGCGGAGTGTCCCTTTCGCGGGGCGAAGCGGCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTc                                                                       <  2:5742182/101‑1 (MQ=38)
                                cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCg                                                                      >  2:4309167/1‑101 (MQ=255)
                                cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCg                                                                      >  2:4660863/1‑101 (MQ=255)
                                   gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                   <  1:4569561/101‑1 (MQ=255)
                                    cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                  >  2:1004178/1‑101 (MQ=255)
                                     gggttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTccc                                                                 >  1:6875181/4‑101 (MQ=255)
                                       aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                               >  2:1520442/1‑101 (MQ=255)
                                        ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                              <  2:3315763/101‑1 (MQ=255)
                                        ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                              >  1:2060978/1‑101 (MQ=255)
                                              cTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                        >  2:1408783/1‑101 (MQ=255)
                                              cTTTCGCGGGACGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                        >  2:356006/1‑101 (MQ=255)
                                               tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                       >  1:1521216/1‑101 (MQ=255)
                                               tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                       >  1:6582516/1‑101 (MQ=255)
                                                ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                      >  2:344877/1‑101 (MQ=255)
                                                  cgcgGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                    <  1:6265978/101‑1 (MQ=255)
                                                      gggCGAGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  1:5028677/101‑1 (MQ=255)
                                                      gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  2:2771258/101‑1 (MQ=255)
                                                       tcactcgagccGCTCCCCCCGGGGGGGCGTGGTGCAAAACATCGGCTTGTCCGCTTCGAA.CCTGGTCGC.GTCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                <  2:2322139/91‑1 (MQ=2)
                                                       ggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  2:5590471/101‑1 (MQ=255)
                                                        gCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                              <  1:192187/101‑1 (MQ=255)
                                                         cGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             <  1:436388/101‑1 (MQ=255)
                                                          gctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             >  1:3129202/4‑101 (MQ=2)
                                                          gctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             <  1:4546228/98‑1 (MQ=2)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  1:105771/1‑101 (MQ=255)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  1:673089/1‑101 (MQ=255)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  2:4055764/1‑101 (MQ=255)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  1:461645/1‑101 (MQ=255)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGCCGCTCCCCGCGCGGGGGCGTGGATcagt                                            >  1:1805909/1‑98 (MQ=255)
                                                            tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCTATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaac                                           >  1:5100689/2‑100 (MQ=2)
                                                            tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                           >  1:6690887/2‑101 (MQ=2)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                          <  1:5658248/101‑1 (MQ=2)
                                                             gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                          <  2:6726750/101‑1 (MQ=2)
                                                              cggcGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCTGGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                        <  2:6816205/98‑1 (MQ=25)
                                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                        >  2:2073894/1‑101 (MQ=255)
                                                                  ccTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTAt                                    <  2:1867012/100‑1 (MQ=255)
                                                                      cccccGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGAt                                <  2:523025/101‑1 (MQ=255)
                                                                         ccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTg                             <  1:6005930/101‑1 (MQ=255)
                                                                          cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTgg                            <  2:4579664/101‑1 (MQ=255)
                                                                            ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                          <  1:4895420/100‑1 (MQ=255)
                                                                            cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                          >  2:1231782/1‑101 (MQ=255)
                                                                             gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                         >  2:6138976/1‑101 (MQ=255)
                                                                             gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcgc                         >  2:4899301/1‑101 (MQ=255)
                                                                              ggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCa                        >  1:1659862/1‑101 (MQ=255)
                                                                                 gCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                     <  1:3238987/101‑1 (MQ=255)
                                                                                 gCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                     <  1:864803/101‑1 (MQ=255)
                                                                                 gCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                     <  2:3035281/101‑1 (MQ=255)
                                                                                   gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                   >  1:1854921/1‑101 (MQ=255)
                                                                                   gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                   >  1:2791315/1‑101 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                  >  2:3894949/1‑101 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                  >  1:1264640/1‑101 (MQ=255)
                                                                                    gggATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                  <  2:6298452/100‑1 (MQ=255)
                                                                                     ggATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                 <  2:6506740/101‑1 (MQ=255)
                                                                                     ggATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaa                 >  2:5874445/1‑101 (MQ=255)
                                                                                         cGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc             <  1:3937991/101‑1 (MQ=255)
                                                                                          gAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt            >  2:2872474/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  1:4041107/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  1:6337534/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           <  1:356010/101‑1 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg           >  2:3987540/1‑101 (MQ=255)
                                                                                            ccgaTCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt          <  2:1195826/98‑1 (MQ=255)
                                                                                            aaCATCTGCATGTAGGCGGGG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt          <  1:1107614/101‑1 (MQ=255)
                                                                                            aaCATCTGCATGTAGGCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt          <  1:4309389/101‑1 (MQ=255)
                                                                                               gtccgCATGTAGCCGGGG.G.GAAGGTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt       <  2:6886098/97‑1 (MQ=25)
                                                                                                tCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc      <  1:2349569/101‑1 (MQ=255)
                                                                                                 cTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGGCTGtcgtcg     >  2:4103895/1‑101 (MQ=255)
                                                                                                 cTGCATGGAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGtcgtcg     <  1:5368326/101‑1 (MQ=255)
                                                                                                  tGCATGTGGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc    <  2:1619603/101‑1 (MQ=255)
                                                                                                  tGCATGTAGGCGGCC.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGc    <  2:1802425/101‑1 (MQ=255)
                                                                                                   gCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt   >  2:1930570/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTC  >  NC_014034/1395570‑1395770

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