Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,538 G→C 5.8% intergenic (+350/‑581) 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,5380G→C 5.8% 46.4 ‑0.9 ‑0.1414intergenic (+350/‑581)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/10):  ref base (189/200):  total (204/210)
Fisher's exact test strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.54e-01

Only 100 of 441 total aligned reads displayed.

CTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCT  >  NC_014034/1395438‑1395638
                                                                                                    |                                                                                                       
cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                                                                                         <  2:3172944/101‑1 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGc                                                                                                      <  1:804068/101‑1 (MQ=255)
    ccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGcc                                                                                                     >  1:3170801/1‑101 (MQ=255)
       cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTc                                                                                                  <  2:201391/101‑1 (MQ=255)
        gcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCACCGCGGGGGCGTGGATCGAAACATCTCGGTGCGGTCg                                                                                                 >  1:5560049/1‑101 (MQ=255)
          gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                                                                               <  2:2649798/101‑1 (MQ=255)
               cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tcc                                                                                         <  1:3384543/101‑1 (MQ=255)
               cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tcc                                                                                         >  2:4726316/1‑101 (MQ=255)
                gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tccc                                                                                        >  1:34793/1‑100 (MQ=255)
                gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                        <  2:3027110/101‑1 (MQ=255)
                gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                        <  2:5410907/101‑1 (MQ=255)
                 tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                       <  1:710469/101‑1 (MQ=255)
                 tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                       <  1:2684305/101‑1 (MQ=255)
                 tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGCGCCGTCGCTCG.Tcccg                                                                                       >  1:92631/1‑99 (MQ=255)
                        aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                               >  1:2930073/1‑101 (MQ=255)
                         aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                             <  2:4889528/101‑1 (MQ=255)
                         aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                             <  2:5751698/101‑1 (MQ=255)
                          aCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCCGCCGG.CCT.Gc                                                                            >  1:3661510/1‑101 (MQ=255)
                           cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCt                                                                           >  1:1446338/1‑101 (MQ=255)
                            aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtc                                                                          <  2:835649/101‑1 (MQ=255)
                             tCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcg                                                                         >  1:2397782/1‑101 (MQ=255)
                               gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcgtc                                                                       <  1:2747905/101‑1 (MQ=255)
                                ccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtcgtcg                                                                      >  1:3349507/1‑101 (MQ=255)
                                   cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTc                                                                   <  2:2944359/101‑1 (MQ=255)
                                    tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcc                                                                  <  1:2138662/101‑1 (MQ=255)
                                       ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTccccc                                                               <  1:3578388/101‑1 (MQ=255)
                                        gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcccccc                                                              <  1:5053436/101‑1 (MQ=255)
                                           gTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcgcg                                                           <  2:4555794/101‑1 (MQ=255)
                                             cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggg                                                         <  1:2234273/101‑1 (MQ=255)
                                             cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggg                                                         >  1:1821655/1‑101 (MQ=255)
                                               aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggggg                                                       >  1:2289480/1‑101 (MQ=255)
                                                 aGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCg                                                     >  2:5664604/1‑101 (MQ=255)
                                                 aGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCGCCCCCCGCGGGGGCg                                                     <  2:3661282/101‑1 (MQ=255)
                                                    gccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTgg                                                  <  2:173591/91‑1 (MQ=0)
                                                    ccggTACCGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCTCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTg                                                   <  1:4350015/99‑1 (MQ=2)
                                                    cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCCCGGGGGCGTgg                                                  >  2:3437753/1‑101 (MQ=255)
                                                    cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTAGCTCG.TCCAGCCGG.CCT.GATCGTCGCTCCCCCCCCGGGGGCGTgg                                                  >  1:4866181/1‑101 (MQ=255)
                                                    aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTgg                                                  >  2:1415721/11‑101 (MQ=2)
                                                     tGGTAGCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGa                                                 <  1:4078514/101‑1 (MQ=255)
                                                     aggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGa                                                 >  2:204898/10‑101 (MQ=0)
                                                      ggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGac                                                >  2:1010127/9‑100 (MQ=0)
                                                       gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCTTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                               <  2:2806348/101‑1 (MQ=255)
                                                         cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             >  2:3521360/6‑101 (MQ=2)
                                                         cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             <  1:1661948/96‑1 (MQ=2)
                                                         cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             <  2:2502632/96‑1 (MQ=2)
                                                         cggcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                             >  2:4672455/6‑101 (MQ=2)
                                                          ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                            >  2:2305697/5‑101 (MQ=2)
                                                           gtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                           <  1:5645155/98‑1 (MQ=2)
                                                           cacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                           <  2:654011/98‑1 (MQ=2)
                                                           cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                           <  2:5019911/101‑1 (MQ=255)
                                                           cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                           >  1:3754592/1‑101 (MQ=255)
                                                            acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAac                                          >  2:4506727/3‑101 (MQ=2)
                                                             ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                         >  2:1591659/2‑100 (MQ=2)
                                                             cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                          <  2:1837263/100‑1 (MQ=2)
                                                              gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacca                                        <  1:156012/101‑1 (MQ=255)
                                                              gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactg                                        >  2:1919656/1‑99 (MQ=2)
                                                              gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacta                                        <  1:4637550/101‑3 (MQ=2)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGTGCGTGGATCGAAacccc                                       >  1:682591/1‑99 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                      <  2:1309263/101‑4 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCt                                     >  2:1524406/1‑101 (MQ=255)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                     >  1:5324871/1‑97 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                     >  1:4472075/1‑97 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgc                                     <  2:2199938/101‑6 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtg                                    >  1:2856248/1‑95 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                   >  1:2072009/1‑94 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                   >  1:2451242/1‑94 (MQ=2)
                                                                    ccccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGc                                  >  1:257802/1‑101 (MQ=255)
                                                                    ccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                                  <  2:5563063/101‑8 (MQ=2)
                                                                    ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                                  >  1:110422/1‑93 (MQ=2)
                                                                    ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                                  >  2:4085895/1‑93 (MQ=2)
                                                                    ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                                  <  1:1324057/101‑9 (MQ=2)
                                                                     cccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                                 >  2:2585378/1‑101 (MQ=255)
                                                                     cccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcgc                                 >  1:4199227/1‑92 (MQ=2)
                                                                     cccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcgccg                                 >  2:766770/1‑93 (MQ=0)
                                                                     cccccGCGCGGGCTTGGGTCGAAGCGTCTCGATGCCGTCGCCCT.TCCAGCCGG.CCT.GCTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                                 <  2:5293530/101‑1 (MQ=25)
                                                                      ccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCgg                                <  1:1291006/101‑1 (MQ=255)
                                                                       cccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGa                               <  1:3930898/101‑1 (MQ=255)
                                                                       cccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGAGAGAAACCACCTTGCGGa                               >  2:5201106/1‑101 (MQ=255)
                                                                         cgcgGGGGCGTGGATCGAAACATCTCGGTGACGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAtt                             >  1:5144688/1‑101 (MQ=255)
                                                                            gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                          >  1:4679000/1‑101 (MQ=255)
                                                                             ggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCa                         <  2:2860748/101‑1 (MQ=255)
                                                                              gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                        <  2:4164540/101‑1 (MQ=255)
                                                                              gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                        <  1:2198476/101‑1 (MQ=255)
                                                                               ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCC.TCCCGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                       <  1:4954920/101‑1 (MQ=255)
                                                                                gCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                      <  1:1006253/101‑1 (MQ=255)
                                                                                    ggATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgc                  >  2:5073312/1‑101 (MQ=255)
                                                                                     gATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg                 >  1:2467323/1‑101 (MQ=255)
                                                                                       tCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCggg               <  2:120903/101‑1 (MQ=255)
                                                                                         gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc             >  2:1098824/1‑101 (MQ=255)
                                                                                         gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc             >  1:4745816/1‑101 (MQ=255)
                                                                                          aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg            >  1:5226566/1‑101 (MQ=255)
                                                                                          aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg            <  2:431107/101‑1 (MQ=255)
                                                                                           aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa           <  2:3839898/101‑1 (MQ=255)
                                                                                            aCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGaa          <  2:92620/101‑1 (MQ=255)
                                                                                              aTCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGc        >  2:4785141/1‑101 (MQ=255)
                                                                                                ctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGt      <  2:2397648/101‑1 (MQ=255)
                                                                                                ctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGt      <  1:1060393/101‑1 (MQ=255)
                                                                                                 tcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCACCTTGCGGATTGTGACTTTCGCGGGGCGAAGCGTc     >  1:1917187/1‑101 (MQ=255)
                                                                                                   ggTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGc   >  1:5761933/1‑101 (MQ=255)
                                                                                                    gTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCt  >  2:4391825/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
CTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCT  >  NC_014034/1395438‑1395638

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