Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,698 Δ1 bp 5.4% intergenic (+510/‑421) 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,6980A→. 5.4% 54.4 ‑0.7 0.0409intergenic (+510/‑421)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/9):  ref base (192/195):  total (205/204)
Fisher's exact test strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 440 total aligned reads displayed.

CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798
                                                                                                       |                                                                                                       
ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCa                                                                                                         <  1:134279/101‑1 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCg                                                                                                     >  2:4660863/1‑101 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCg                                                                                                     <  1:6530957/101‑1 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCg                                                                                                     <  1:1846292/101‑1 (MQ=255)
       gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                                                  >  2:5014717/1‑101 (MQ=255)
       gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTc                                                                                                  <  2:2548758/101‑1 (MQ=255)
        gggATTGTCACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                                 <  2:5892530/100‑1 (MQ=255)
         ggATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTccc                                                                                                >  1:2916539/1‑101 (MQ=255)
           aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                                                              >  1:2164004/1‑101 (MQ=255)
            ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                                             <  2:476159/101‑1 (MQ=255)
            ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                                             <  2:1891281/101‑1 (MQ=255)
            ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                                                             <  2:3452774/101‑1 (MQ=255)
              gTCGCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                                           >  1:6825014/1‑101 (MQ=255)
                  cTTTCGCGGGACGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgggg                                                                                       >  2:356006/1‑101 (MQ=255)
                   tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                                      >  1:4673702/1‑101 (MQ=255)
                    ttCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGGCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                                     <  1:5344912/101‑1 (MQ=255)
                      cgcgGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGt                                                                                   <  1:6265978/101‑1 (MQ=255)
                        ggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                                  >  1:5853814/10‑101 (MQ=2)
                        cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  1:4270277/101‑1 (MQ=255)
                          gggCGAGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  1:5028677/101‑1 (MQ=255)
                          gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  1:4010549/101‑1 (MQ=255)
                           tcactcgagccGCTCCCCCCGGGGGGGCGTGGTGCAAAACATCGGCTTGTCCGCTTCGAA.CCTGGTCGC.GTCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  2:2322139/91‑1 (MQ=2)
                           ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:2457405/101‑1 (MQ=255)
                           ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGGAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:2216234/101‑1 (MQ=255)
                           gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               >  1:4318858/7‑101 (MQ=2)
                           gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  2:4224277/95‑1 (MQ=2)
                            gCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  1:192187/101‑1 (MQ=255)
                            ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:5974509/96‑1 (MQ=2)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  2:1940605/1‑101 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  2:4055764/1‑101 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  1:105771/1‑101 (MQ=255)
                               aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                          >  2:1652002/1‑101 (MQ=255)
                                tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG.GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                          >  2:3942573/2‑101 (MQ=2)
                                  cggcGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCTGGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      <  2:6816205/98‑1 (MQ=25)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      >  2:2073894/1‑101 (MQ=255)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gg                                                                     <  1:6899952/101‑1 (MQ=255)
                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCATTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                    <  1:4585699/101‑1 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   <  2:2804532/101‑1 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   >  1:2554573/1‑101 (MQ=255)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   <  2:142234/101‑1 (MQ=255)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGCGGATCGGAACT.GGTAt                                                                  >  2:5088065/1‑101 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCTTGTCCGCATCGAA.GCTGGTCAC.ATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcggcgta                                                                   <  1:4558574/101‑7 (MQ=2)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                 <  2:3714038/101‑1 (MQ=255)
                                        tCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGCCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                <  1:832510/101‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGa                                                               <  1:144901/101‑1 (MQ=255)
                                          cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                              <  2:1161367/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            >  2:1574722/1‑101 (MQ=255)
                                                ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgcg                                                        <  1:4895420/100‑1 (MQ=255)
                                                 gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgcgc                                                       >  2:4899301/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgcgc                                                       >  2:6138976/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgcgc                                                       >  2:3158682/1‑101 (MQ=255)
                                                   gggGGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCT.GGTATGGATCTGGCGCaa                                                     >  1:5621309/1‑101 (MQ=255)
                                                     gCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAAc                                                   <  1:159662/101‑1 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  1:6641212/1‑101 (MQ=255)
                                                        tGGATCAAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..a                                              <  2:1219762/101‑1 (MQ=255)
                                                        gggATCGAAACATCGGCATGGAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..a                                              <  1:1677431/100‑1 (MQ=255)
                                                         ggATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aa                                             <  2:6506740/101‑1 (MQ=255)
                                                         ggATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aa                                             <  1:4514480/101‑1 (MQ=255)
                                                          gATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..aaa                                            <  2:6263602/101‑1 (MQ=255)
                                                             cGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTc                                         >  2:1790938/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  2:1296177/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  2:2872474/1‑101 (MQ=255)
                                                               aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTg                                       <  1:356010/101‑1 (MQ=255)
                                                               aaaCATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTg                                       >  1:6337534/1‑101 (MQ=255)
                                                               aaaCATCTCCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTg                                       <  1:3062298/101‑1 (MQ=255)
                                                                accttCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      <  1:4079946/97‑1 (MQ=255)
                                                                aaCATCTGCATGTAGTCGGGG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      <  1:2632981/101‑1 (MQ=255)
                                                                aaCATCTGCATGTAGGCGGGG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      <  1:1107614/101‑1 (MQ=255)
                                                                 aCATCTGCATGTAGTCGGCG.G.GCAGTGCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtc                                     <  1:5130700/101‑1 (MQ=255)
                                                                  cATTTGCATGTAGTCGGCG.G.GCAGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    <  1:3293740/101‑1 (MQ=255)
                                                                  cATCTGCATGTGGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    <  2:2590545/101‑1 (MQ=255)
                                                                   aTCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  2:5259916/101‑1 (MQ=255)
                                                                   aTCTGCATGTAGGCGGGG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  1:6077143/101‑1 (MQ=255)
                                                                       gCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtagatcgg                               >  1:1195874/1‑94 (MQ=255)
                                                                       gCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCt                               >  1:3228292/1‑101 (MQ=255)
                                                                       gCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCt                               >  1:6718741/1‑101 (MQ=255)
                                                                       gCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCt                               >  2:5623958/1‑101 (MQ=255)
                                                                          tGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccc                            >  1:5242305/1‑101 (MQ=255)
                                                                             agccgggg.g.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         <  2:6594548/98‑1 (MQ=255)
                                                                             aGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         >  2:6751660/1‑101 (MQ=255)
                                                                             aGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         <  2:4557053/101‑1 (MQ=255)
                                                                                 ggcg.g.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCgg                     >  1:3919131/1‑101 (MQ=255)
                                                                                         cAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGt               >  1:3737528/1‑101 (MQ=255)
                                                                                         cAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGt               >  2:2974812/1‑101 (MQ=255)
                                                                                              tcaagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTg              >  2:4157613/9‑101 (MQ=2)
                                                                                              cGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATc          <  2:6337210/101‑1 (MQ=255)
                                                                                               caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTgg             <  2:3354468/94‑1 (MQ=2)
                                                                                                 agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGAt           <  2:6875521/96‑1 (MQ=2)
                                                                                                 agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGAt           >  2:707062/6‑101 (MQ=2)
                                                                                                  gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  1:234772/101‑1 (MQ=255)
                                                                                                  gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GATATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:6030080/101‑1 (MQ=255)
                                                                                                   ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  2:6698815/101‑1 (MQ=255)
                                                                                                   attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         >  1:98597/4‑101 (MQ=2)
                                                                                                   attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  1:7136359/98‑1 (MQ=2)
                                                                                                     tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       >  2:1834046/2‑101 (MQ=2)
                                                                                                     tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:80669/100‑1 (MQ=2)
                                                                                                     gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc    <  1:3442421/101‑1 (MQ=255)
                                                                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     >  2:6917555/1‑101 (MQ=2)
                                                                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg  <  2:5771389/101‑1 (MQ=255)
                                                                                                       aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg  <  1:1270707/101‑1 (MQ=255)
                                                                                                       |                                                                                                       
CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG.G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG