Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,024 C→A 21.0% intergenic (+836/‑95) 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,396,0240C→A21.0% 152.1 ‑0.1 ‑0.1244intergenic (+836/‑95)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (25/26):  ref base (91/101):  total (116/128)
Fisher's exact test strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.78e-01

Only 100 of 272 total aligned reads displayed.

TCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAT  >  NC_014034/1395924‑1396124
                                                                                                    |                                                                                                       
  gAAACGATGCGTTCATGGGCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                                       <  1:4731998/101‑1 (MQ=255)
  gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCCCGCCt                                                                                                       >  1:2060740/1‑101 (MQ=255)
  gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                                       >  1:2890823/1‑101 (MQ=255)
  gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                                       >  2:3729097/1‑101 (MQ=255)
      cGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.cg                                                                                                  >  1:3112429/1‑101 (MQ=255)
          gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                                                                            <  2:2312661/101‑1 (MQ=255)
          gCGTTCATGGCCTACGTCGGGCAGGGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                                                                            <  2:879863/101‑1 (MQ=255)
           cGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCg                                                                                           <  2:1088789/101‑1 (MQ=255)
             ttCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCt                                                                                         <  2:286958/101‑1 (MQ=255)
              tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                                                                                        >  1:3875343/1‑101 (MQ=255)
              tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                                                                                        >  1:2293382/1‑101 (MQ=255)
              tCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTc                                                                                        >  1:4516713/1‑101 (MQ=255)
                aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                                                                                      >  2:40314/1‑101 (MQ=255)
                 tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                                                                                     >  2:2472438/1‑101 (MQ=255)
                       aCGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCgg                                                                               >  2:1053877/1‑101 (MQ=255)
                          tCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCggggg                                                                            <  2:2184017/101‑1 (MQ=255)
                           cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc                                                                           <  2:2525020/101‑1 (MQ=255)
                             ggCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCgg                                                                         >  1:2619908/1‑100 (MQ=255)
                              gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGCCCTCGATCCGCGCTTGCCTCACCCCTGC.CGC..GGCGCTCCCCCCGCggggggcgg                                                                        >  1:850513/1‑97 (MQ=38)
                                aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      <  2:832172/101‑1 (MQ=255)
                                aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      >  2:484475/1‑101 (MQ=255)
                                 tctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                      >  1:392571/8‑101 (MQ=2)
                                 tctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                      <  2:241647/94‑1 (MQ=2)
                                 gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     >  2:3939381/1‑101 (MQ=255)
                                 gTGGACTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     <  1:1004555/101‑1 (MQ=255)
                                  ctgcgcggTCGCTCCCCCCGCGGGTGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     <  1:1872634/95‑1 (MQ=2)
                                   tgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                    <  2:3794143/96‑1 (MQ=9)
                                   ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   <  1:829343/101‑1 (MQ=255)
                                   ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   >  2:3490451/1‑101 (MQ=255)
                                   ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   >  2:2995570/1‑101 (MQ=255)
                                    gtcTGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCCCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGa                                                                  <  2:850439/99‑1 (MQ=255)
                                    gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   <  2:3100275/97‑1 (MQ=9)
                                     gcTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCGCCCGCGGGGGCGTGGATCGaa                                                                 <  1:2689300/100‑1 (MQ=255)
                                     cggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                  <  2:3611578/97‑1 (MQ=2)
                                     aCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGATACTCGCCCTAGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                 <  1:4122714/101‑1 (MQ=255)
                                     aCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                 <  2:1552378/101‑1 (MQ=255)
                                     aCTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                 <  2:962473/101‑1 (MQ=255)
                                     aCTGGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                 <  1:2488023/101‑1 (MQ=255)
                                       cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                <  1:4429291/100‑1 (MQ=2)
                                        ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                               >  1:1548499/1‑101 (MQ=9)
                                        ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                               >  1:3727439/1‑101 (MQ=9)
                                        ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                              >  1:3021149/1‑101 (MQ=255)
                                        ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                              <  2:2222146/101‑1 (MQ=255)
                                         gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                             <  1:1331592/101‑1 (MQ=255)
                                           cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtg                                                            >  1:3519283/1‑98 (MQ=2)
                                            gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcg                                                          <  1:606474/101‑1 (MQ=255)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACTGGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCGCCCGCGGGGGCGCGGATCGAAACAgggga                                                         >  2:3121550/1‑98 (MQ=38)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgat                                                          >  1:1275756/1‑96 (MQ=2)
                                             cTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcgc                                                         <  1:656913/101‑1 (MQ=255)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                         >  2:4324036/1‑95 (MQ=2)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                         <  2:666053/101‑7 (MQ=2)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                         <  2:1626684/101‑7 (MQ=2)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                        >  2:1073545/1‑101 (MQ=255)
                                              tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                        >  1:4230754/1‑101 (MQ=255)
                                               ccccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatga                                                        >  1:3268973/1‑94 (MQ=2)
                                               ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                                       >  2:1980772/1‑101 (MQ=255)
                                                cccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgac                                                       >  1:1927630/1‑93 (MQ=2)
                                                  cccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacca                                                     <  2:1834208/101‑11 (MQ=2)
                                                  cccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTg                                                    >  1:3638546/1‑101 (MQ=255)
                                                   ccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGa                                                   <  1:1794447/101‑1 (MQ=255)
                                                    cgcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGAt                                                  >  2:2905250/1‑101 (MQ=255)
                                                    cgcgGGGGCGTGGATCGAAACCCGCCCTCTATCCGCGCTTGCCCCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGAt                                                  <  2:4399515/101‑1 (MQ=38)
                                                     gcgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATg                                                 <  2:3229130/101‑1 (MQ=255)
                                                      cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                                <  1:4647008/101‑1 (MQ=255)
                                                      cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                                <  2:1143967/101‑1 (MQ=255)
                                                      cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                                <  2:4381624/101‑1 (MQ=255)
                                                       gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                               <  2:4721142/101‑1 (MQ=255)
                                                       gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                               <  2:2619756/101‑1 (MQ=255)
                                                           gggTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                           <  1:3056332/99‑1 (MQ=255)
                                                           gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                           <  1:776149/101‑1 (MQ=255)
                                                           gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                           <  1:848214/101‑1 (MQ=255)
                                                            cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                          >  2:4649610/1‑101 (MQ=255)
                                                             gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                                         >  1:2129515/1‑101 (MQ=255)
                                                               ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                                       >  1:866564/1‑101 (MQ=255)
                                                                  tCGAAGCTCGCCCTGGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTg                                    <  1:2284985/101‑1 (MQ=255)
                                                                    gacaCTCCCCCTCGATCCGCGCTTGCCTCGCGCCTGC.CGC..GTCTCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                  <  1:336498/98‑1 (MQ=38)
                                                                     aaaCTCGCCCTGGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACAGGTGCATCTGATGATTATGTCCGGTgaga                                 >  2:3149029/1‑101 (MQ=37)
                                                                       aCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                               <  1:1118802/101‑1 (MQ=255)
                                                                         tCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGc                             >  1:4035273/1‑101 (MQ=255)
                                                                           ggccTCGATCCGCGCTTGGCTCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTg                           <  2:1440003/99‑1 (MQ=255)
                                                                           gCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTg                           >  1:4237785/1‑101 (MQ=255)
                                                                             ccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                         <  1:1367319/101‑1 (MQ=255)
                                                                                cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt                      <  2:2723795/101‑1 (MQ=255)
                                                                                  aTCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcc                    <  1:177652/101‑1 (MQ=255)
                                                                                   tCCGCGCTTGCCTCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTccc                   <  2:742785/101‑1 (MQ=255)
                                                                                   tCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTNAGACGGCTGTCGCTcnc                   >  2:1508748/1‑101 (MQ=255)
                                                                                      gcgcTTGCCTCACGCCTGC.CGC..GTCGGTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc                <  1:3989001/101‑1 (MQ=255)
                                                                                       cgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcg               >  2:3389751/1‑101 (MQ=255)
                                                                                       cccTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcg               <  2:1202749/99‑1 (MQ=255)
                                                                                         cTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcgcg             <  1:3736485/101‑1 (MQ=255)
                                                                                           tGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCggg           >  2:1354183/1‑101 (MQ=255)
                                                                                              cTCACGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGc        <  1:1053956/101‑1 (MQ=255)
                                                                                               tCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATGTGACTTAGATGCCCGTATGGAGGACGGGCGGTCCCCCCCCggggggg       >  2:2267115/1‑99 (MQ=12)
                                                                                                cACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGt      >  2:172877/1‑101 (MQ=255)
                                                                                                 aCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTg     <  2:2367257/101‑1 (MQ=255)
                                                                                                   gCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGa   <  1:1825089/101‑1 (MQ=255)
                                                                                                   gCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGa   <  2:2129394/101‑1 (MQ=255)
                                                                                                   gCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGa   <  2:1807739/101‑1 (MQ=255)
                                                                                                    ccTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAt  <  1:2372193/101‑1 (MQ=255)
                                                                                                    ccTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAt  >  1:4634411/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
TCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGGAT  >  NC_014034/1395924‑1396124

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