Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,549 C→A 26.5% intergenic (+361/‑570) 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,5490C→A26.5% 265.9 ‑1.9 ‑1.3355intergenic (+361/‑570)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (35/59):  ref base (121/139):  total (157/198)
Fisher's exact test strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.08e-02

Only 100 of 380 total aligned reads displayed.

GGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGG  >  NC_014034/1395449‑1395649
                                                                                                    |                                                                                                       
ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                                                                                         <  2:1107056/101‑1 (MQ=255)
ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                                                                                         >  1:222972/1‑101 (MQ=255)
    cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tcc                                                                                                    <  2:3931898/101‑1 (MQ=255)
      tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                                  <  2:3637813/101‑1 (MQ=255)
       ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGc                                                                                                 >  2:2635514/1‑101 (MQ=255)
           cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGGCGTCGCTCG.TCCAGCCGG.c                                                                                            >  2:1765393/1‑101 (MQ=255)
            gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                                           >  2:3227081/1‑101 (MQ=255)
             aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                                          >  2:151340/1‑101 (MQ=255)
                cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCt                                                                                      >  2:2980976/1‑101 (MQ=255)
                 aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.CCCAGCCGG.CCT.GCtc                                                                                     >  1:1076904/1‑101 (MQ=255)
                       gcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCt                                                                               >  1:3657427/1‑101 (MQ=255)
                        cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTc                                                                              >  1:2256477/1‑101 (MQ=255)
                        cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.CCCAGCCGG.CCT.GCTCGTCGCTc                                                                              >  1:1832802/1‑101 (MQ=255)
                          ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCCCTccc                                                                            >  1:4972627/1‑101 (MQ=255)
                            ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTccccc                                                                          >  2:4443925/1‑101 (MQ=255)
                             gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcccccc                                                                         <  2:2148512/101‑1 (MQ=255)
                               ggTCGGCGGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcgc                                                                       <  1:3934919/101‑1 (MQ=255)
                               ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcgc                                                                       >  2:2762534/1‑101 (MQ=255)
                                  cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggg                                                                    <  2:34000/101‑1 (MQ=255)
                                    aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggggg                                                                  >  1:3094561/1‑101 (MQ=255)
                                          tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  1:3727403/1‑101 (MQ=255)
                                          tGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            <  2:2372520/101‑1 (MQ=255)
                                          gcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCACCCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            <  2:1661797/92‑1 (MQ=0)
                                          gcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            <  2:1302225/92‑1 (MQ=0)
                                          gcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            <  1:477855/92‑1 (MQ=0)
                                           ggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:28874/93‑1 (MQ=2)
                                           ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           >  1:2300646/1‑101 (MQ=255)
                                           cctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGGGGAt                                                           >  1:3464217/9‑101 (MQ=2)
                                            gTAACGTTTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  2:2247891/1‑101 (MQ=255)
                                            gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGTATc                                                          <  2:4745741/101‑1 (MQ=255)
                                            gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  1:4092647/101‑1 (MQ=255)
                                            ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  1:1558154/94‑1 (MQ=2)
                                             tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  1:2836702/101‑1 (MQ=255)
                                             tAACGTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         >  2:3358588/1‑101 (MQ=255)
                                             tAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  2:5693/101‑1 (MQ=255)
                                             ccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  1:3778561/95‑1 (MQ=2)
                                              ctgcggtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  1:1754332/96‑1 (MQ=2)
                                              cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  1:490431/96‑1 (MQ=2)
                                              accGTGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  2:2435406/99‑1 (MQ=255)
                                              aaCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  1:3342490/101‑1 (MQ=255)
                                               ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGGGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       >  2:1945034/5‑101 (MQ=2)
                                               gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                       >  1:4421915/5‑101 (MQ=2)
                                                cgcgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       <  2:4565951/98‑1 (MQ=2)
                                                cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCT.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  1:1535899/1‑101 (MQ=255)
                                                 tggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                     <  2:2170493/99‑1 (MQ=2)
                                                 gtgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAac                                                     >  1:895680/1‑101 (MQ=255)
                                                  cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAGCAATGCGGCGCGGGC.CTAT.GGCTGCCGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                                    <  1:1650046/100‑2 (MQ=2)
                                                  cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                     <  2:369475/100‑1 (MQ=2)
                                                   gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAacca                                                   >  2:3888596/1‑101 (MQ=255)
                                                   gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccc                                                   <  2:3134417/101‑2 (MQ=2)
                                                   gtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactg                                                   <  1:4639615/101‑3 (MQ=2)
                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacccc                                                  <  1:4543318/101‑3 (MQ=2)
                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                                  <  1:2016563/101‑4 (MQ=2)
                                                    gcGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                                  <  2:1057343/100‑4 (MQ=2)
                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                                 >  2:2544048/1‑97 (MQ=2)
                                                      gccccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                                <  1:4462868/98‑6 (MQ=2)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCACCtt                                               >  1:1927380/1‑101 (MQ=255)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccccg                                               >  2:908102/1‑99 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                              >  1:549801/1‑94 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              >  1:1176760/1‑94 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              <  1:216669/101‑8 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              <  1:2438192/101‑8 (MQ=2)
                                                        cccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                              <  2:4506471/100‑8 (MQ=2)
                                                         ccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcg                                             <  1:1300738/101‑9 (MQ=2)
                                                         ccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                                             >  2:4398618/1‑94 (MQ=2)
                                                         ccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                                             >  1:780972/1‑94 (MQ=2)
                                                          cccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                                            <  1:2983954/101‑1 (MQ=255)
                                                          cccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                                            >  2:1718024/1‑101 (MQ=255)
                                                          cccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggcg                                            >  2:3079507/1‑93 (MQ=0)
                                                           ccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcgcg                                           <  2:2661674/101‑11 (MQ=2)
                                                           ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                                           >  1:3863349/1‑91 (MQ=2)
                                                           ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                                           <  2:1590595/101‑11 (MQ=2)
                                                            cccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGa                                          <  2:5011985/101‑1 (MQ=255)
                                                            cccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggcgtc                                          <  2:3830578/101‑11 (MQ=2)
                                                              cgcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAtt                                        <  1:336927/101‑1 (MQ=255)
                                                              cgcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAtt                                        <  1:4838087/101‑1 (MQ=255)
                                                                tgggggCGCGAATCAAAACATCTCCGAACCGTCCCACG.TCCAGCCGG.CCT.ACTCGTCGCACCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGt                                      <  2:2859752/100‑1 (MQ=25)
                                                                 gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCCTGCGGATTGtg                                     >  2:2080098/1‑100 (MQ=255)
                                                                 gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                                     <  1:3227332/101‑1 (MQ=255)
                                                                  ggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.GCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCa                                    <  1:608531/101‑1 (MQ=255)
                                                                    ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                                  <  2:149210/101‑1 (MQ=255)
                                                                    ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                                  >  1:684189/1‑101 (MQ=255)
                                                                      ggTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                <  1:4015395/100‑1 (MQ=255)
                                                                      cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                <  1:1018856/101‑1 (MQ=255)
                                                                      cGTGGATCAAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                <  1:1507009/101‑1 (MQ=255)
                                                                       gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                               <  1:1693341/101‑1 (MQ=255)
                                                                        tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                              <  1:2132918/101‑1 (MQ=255)
                                                                           aTCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCgg                           >  1:2587850/1‑101 (MQ=255)
                                                                              gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCAAAACCACCTTGCGGATTGTCACTTTCGCGGGGc                        >  1:4332161/1‑101 (MQ=255)
                                                                                   aTCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGc                   <  1:4267675/101‑1 (MQ=255)
                                                                                    tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGGGAAGCg                  >  1:4566164/1‑101 (MQ=255)
                                                                                    tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAggg                  >  2:2998641/1‑99 (MQ=255)
                                                                                     ctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGt                 >  2:4029655/1‑101 (MQ=255)
                                                                                       cGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCg               <  1:4176050/101‑1 (MQ=255)
                                                                                           gCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcc           >  2:3567174/1‑101 (MQ=255)
                                                                                             cGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcccc         >  1:352900/1‑101 (MQ=255)
                                                                                               tCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcccccc       <  2:4214012/101‑1 (MQ=255)
                                                                                               tCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcccccc       <  2:4951489/101‑1 (MQ=255)
                                                                                               tCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcccccc       <  2:684100/101‑1 (MQ=255)
                                                                                                cGCTCG.CCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCcg      <  2:1832669/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
GGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGG  >  NC_014034/1395449‑1395649

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