Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,536 C→A 9.9% intergenic (+348/‑583) 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,5360C→A 9.9% 67.7 ‑0.6 ‑0.0302intergenic (+348/‑583)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/17):  ref base (134/138):  total (147/155)
Fisher's exact test strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 315 total aligned reads displayed.

CGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTC  >  NC_014034/1395436‑1395635
                                                                                                    |                                                                                                      
cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                                                                        <  1:5830/101‑1 (MQ=255)
cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                                                                        <  1:2853850/101‑1 (MQ=255)
   tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                                                                                     >  1:2066742/1‑101 (MQ=255)
    ccccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGTCCGACAGGCCGTTATAGGGGCGGTCGCCCCGCGGGGCCGGGGTTCGAAACTACTTGGTg                                                                                                    >  2:833080/1‑101 (MQ=16)
       cccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCg                                                                                                 <  1:3261659/101‑1 (MQ=255)
         cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTc                                                                                               >  1:554238/1‑101 (MQ=255)
           cgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                                                                             <  1:158392/101‑1 (MQ=255)
            gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                                                                            >  2:348213/1‑101 (MQ=255)
            gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGGCGCt                                                                                            >  1:1890986/1‑101 (MQ=255)
             ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGGCACTc                                                                                           >  2:1153172/1‑101 (MQ=255)
               ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTCGGTGCCGTCGCTCgc                                                                                         >  1:422710/1‑100 (MQ=255)
                gCGTTGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tc                                                                                       <  1:4161914/101‑1 (MQ=255)
                 cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tcc                                                                                      <  1:1947437/101‑1 (MQ=255)
                  gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                     >  2:1042176/1‑101 (MQ=255)
                   tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                    <  1:1341215/101‑1 (MQ=255)
                      aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCg                                                                                 >  2:3624506/1‑101 (MQ=255)
                      aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCg                                                                                 >  2:3402862/1‑101 (MQ=255)
                        cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCCG.c                                                                              <  1:3128442/101‑1 (MQ=255)
                          aagcATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                            >  1:2418477/4‑101 (MQ=255)
                             cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCt                                                                        >  2:2007715/1‑101 (MQ=255)
                              aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtc                                                                       <  2:1626572/101‑1 (MQ=255)
                              aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtc                                                                       <  2:1567323/101‑1 (MQ=255)
                              aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtc                                                                       >  1:178262/1‑101 (MQ=255)
                               tCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcg                                                                      <  1:2916134/101‑1 (MQ=255)
                                cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcgt                                                                     >  1:2983715/1‑101 (MQ=255)
                                      tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcc                                                               >  1:3211688/1‑101 (MQ=255)
                                       ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTccc                                                              <  2:4207246/101‑1 (MQ=255)
                                        tGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCACCCGG.CCT.GCTCGTCGCTcccc                                                             <  1:2411508/101‑1 (MQ=255)
                                           cGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcg                                                          >  1:3543589/1‑101 (MQ=255)
                                            ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcgc                                                         >  1:2643189/1‑101 (MQ=255)
                                            ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcgc                                                         >  2:4149808/1‑101 (MQ=255)
                                               cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GGTCGTCGCTTCCCCCGCggg                                                      >  2:3205056/1‑101 (MQ=255)
                                               cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggg                                                      >  2:765501/1‑101 (MQ=255)
                                                gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                     >  1:3010406/1‑101 (MQ=255)
                                                gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                     <  2:2151676/101‑1 (MQ=255)
                                                 aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggggg                                                    >  2:2081737/1‑101 (MQ=255)
                                                  cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGc                                                   >  2:1083280/1‑101 (MQ=255)
                                                  cAGGCTGGTAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGc                                                   <  1:3058214/101‑1 (MQ=255)
                                                      agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTgg                                               <  2:2285179/91‑1 (MQ=0)
                                                      aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTgg                                               >  1:1035369/11‑101 (MQ=2)
                                                        ggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGAt                                             <  2:1984778/93‑1 (MQ=2)
                                                        ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                             >  1:391596/1‑101 (MQ=255)
                                                         ataccgcgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                             <  1:3497960/94‑1 (MQ=0)
                                                          tcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                           >  1:3733577/7‑101 (MQ=2)
                                                          tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                           >  2:546495/1‑101 (MQ=255)
                                                           cgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                          <  2:1713020/96‑1 (MQ=2)
                                                           aaCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                          >  1:1096869/1‑101 (MQ=255)
                                                            ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  2:3191313/97‑1 (MQ=2)
                                                            gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         >  1:1743777/5‑101 (MQ=2)
                                                             gtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                        >  2:3552237/4‑101 (MQ=2)
                                                               ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                      <  2:961627/100‑1 (MQ=2)
                                                               ggtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAAAACCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                      <  1:1378143/100‑1 (MQ=2)
                                                               cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                      <  2:499947/100‑2 (MQ=2)
                                                                ggcGCGCCCCCCGGGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccc                                     <  2:1806895/99‑2 (MQ=0)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTCGATCGAAaccac                                    <  1:2293636/101‑1 (MQ=255)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                    >  1:1988317/1‑98 (MQ=2)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                    <  2:58969/101‑4 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                   >  2:2358172/1‑98 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                   >  1:3801984/1‑97 (MQ=2)
                                                                  cGCGCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                   <  2:2554911/101‑5 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCt                                  >  2:3322624/1‑101 (MQ=255)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgg                                  >  2:3956885/1‑96 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                  >  2:899049/1‑97 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgc                                  >  1:4004352/1‑96 (MQ=2)
                                                                    cTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgct                                 >  1:3318772/1‑95 (MQ=2)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                >  2:664447/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                >  1:1799319/1‑101 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgggc                                >  1:812489/1‑94 (MQ=2)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                <  1:1225729/101‑8 (MQ=2)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                >  1:2133854/1‑94 (MQ=2)
                                                                     cccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                <  1:3625092/100‑8 (MQ=2)
                                                                      ccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccgcgtgcg                               >  1:444584/1‑94 (MQ=2)
                                                                      ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                               >  1:3345775/1‑93 (MQ=2)
                                                                      ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                               >  1:896318/1‑93 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGNATCGAAactatgctggc                              >  2:561563/1‑92 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                              <  2:2618532/101‑10 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                              <  1:277709/101‑10 (MQ=2)
                                                                        ccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgggcgcg                             >  2:4119558/1‑91 (MQ=2)
                                                                        ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                             <  2:2518090/101‑11 (MQ=2)
                                                                        ccccGCGGGGGCGTGGAGCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcgcg                             <  2:303094/101‑11 (MQ=0)
                                                                         cccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCACCTTGCggg                            >  2:276987/1‑100 (MQ=255)
                                                                          ccGCGGGGGCGTGGATCGAAACATCTCGGCGCCGTCGCTCG.TCCCGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAt                           <  1:2483203/101‑1 (MQ=255)
                                                                             cgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGt                        <  2:1447839/101‑1 (MQ=255)
                                                                              gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTggc                       >  2:1908889/1‑99 (MQ=255)
                                                                              gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                       >  2:1310057/1‑101 (MQ=255)
                                                                              gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                       >  2:2415076/1‑101 (MQ=255)
                                                                              gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTGGCGGATTGTc                       >  1:965997/1‑101 (MQ=255)
                                                                                gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                     <  1:3678978/101‑1 (MQ=255)
                                                                                gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                     <  1:2726450/101‑1 (MQ=255)
                                                                                  gCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                   <  2:3092441/101‑1 (MQ=255)
                                                                                   cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                  <  2:391545/101‑1 (MQ=255)
                                                                                    gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                 >  2:4001114/1‑101 (MQ=255)
                                                                                    gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                 <  1:117110/101‑1 (MQ=255)
                                                                                     tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                <  1:3405409/101‑1 (MQ=255)
                                                                                     tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                >  2:2885313/1‑101 (MQ=255)
                                                                                           gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc          >  2:3873211/1‑101 (MQ=255)
                                                                                            aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcgagatc         >  2:177818/1‑96 (MQ=255)
                                                                                             aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa        >  2:3812846/1‑101 (MQ=255)
                                                                                             aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa        <  2:2478036/101‑1 (MQ=255)
                                                                                                 tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCg    >  1:1151060/1‑101 (MQ=255)
                                                                                                    |                                                                                                      
CGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTC  >  NC_014034/1395436‑1395635

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG