Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,022 C→G 21.6% intergenic (+834/‑97) 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,0220C→G21.6% 138.1 ‑0.1 ‑1.2241intergenic (+834/‑97)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (25/27):  ref base (94/95):  total (119/122)
Fisher's exact test strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.82e-02

Only 100 of 267 total aligned reads displayed.

GATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395922‑1396122
                                                                                                    |                                                                                                       
 aTCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACg                                                                                                        <  2:577179/101‑1 (MQ=255)
  tCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGc                                                                                                       >  1:3537587/1‑101 (MQ=255)
   cGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGcc                                                                                                      >  2:1004476/1‑101 (MQ=255)
    gAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCt                                                                                                     >  2:3729097/1‑101 (MQ=255)
            gCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                                                                          <  2:2312661/101‑1 (MQ=255)
            gCGTTCATGGCCTACGTCGGGCAGGGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTc                                                                                          <  2:879863/101‑1 (MQ=255)
                 cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                                                                                     <  2:275834/101‑1 (MQ=255)
                 cATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcc                                                                                     >  1:4407172/1‑101 (MQ=255)
                  aTGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTccc                                                                                    >  2:40314/1‑101 (MQ=255)
                   tGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccc                                                                                   >  2:2472438/1‑101 (MQ=255)
                     gCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTcccccc                                                                                 >  1:4721398/1‑101 (MQ=255)
                      ccTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCcg                                                                                <  2:4750750/101‑1 (MQ=255)
                          cGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCggg                                                                            >  2:1118716/1‑101 (MQ=255)
                             cGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGc                                                                         >  2:4475097/1‑101 (MQ=255)
                                gCAGTGGCCTGGTCGCTCCCCCCGCGGGGCCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTg                                                                      <  2:2509259/101‑1 (MQ=255)
                                gCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCGCCCTCGATCCGCGCTTGCCTCACCCCTGC.CGC..GGCGCTCCCCCCGCggggggcgg                                                                      >  1:850513/1‑97 (MQ=38)
                                  aGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGa                                                                    <  2:832172/101‑1 (MQ=255)
                                   tctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                    >  1:392571/8‑101 (MQ=2)
                                   tctgcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                    <  2:241647/94‑1 (MQ=2)
                                   gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGTAt                                                                   <  2:3089516/101‑1 (MQ=255)
                                   gTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   >  1:3149491/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)
                                     tgcgcggTCGCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  2:3113748/96‑1 (MQ=2)
                                     ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 >  2:2995570/1‑101 (MQ=255)
                                     ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  1:2328784/101‑1 (MQ=255)
                                     ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 >  2:3490451/1‑101 (MQ=255)
                                     ggACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  1:829343/101‑1 (MQ=255)
                                      gtcTGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCCCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGa                                                                <  2:850439/99‑1 (MQ=255)
                                      gcgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  2:43652/97‑1 (MQ=9)
                                      gACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGCCGTGGATCGa                                                                >  2:2153862/1‑101 (MQ=255)
                                       cggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                <  2:3611578/97‑1 (MQ=2)
                                       cgcggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                >  1:3755318/4‑101 (MQ=9)
                                       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)
                                         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                                                            <  2:2222146/101‑1 (MQ=255)
                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                            <  1:3054694/101‑2 (MQ=2)
                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                           <  2:1814144/101‑1 (MQ=255)
                                           gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAg                                                           <  1:1331592/101‑1 (MQ=255)
                                            tCGCTCCCCCCGCGGGGGCGTTGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActt                                                           >  1:4545124/1‑99 (MQ=2)
                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActt                                                           >  1:291247/1‑99 (MQ=2)
                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtg                                                          >  2:39766/1‑98 (MQ=2)
                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGgcg                                                        >  2:2019108/1‑101 (MQ=255)
                                               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:1626684/101‑7 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                       <  2:666053/101‑7 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatg                                                       >  1:752793/1‑95 (MQ=2)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                      >  1:4230754/1‑101 (MQ=255)
                                                tCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCa                                                      >  2:322147/1‑101 (MQ=255)
                                                 ccccccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCAt                                                     <  1:1826724/101‑1 (MQ=255)
                                                   ccccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacc                                                    <  1:2395349/101‑10 (MQ=2)
                                                    cccGCGGGGGCGTGGATCGAAACTCGGCCTCGATCCG.GTCAGCCTC.GTATTGCGCGCTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAActtgatgacca                                                   >  1:1695958/1‑91 (MQ=2)
                                                    cccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTg                                                  <  2:1942291/101‑1 (MQ=255)
                                                    cccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTg                                                  >  1:1989488/1‑101 (MQ=255)
                                                     ccGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGa                                                 <  1:1794447/101‑1 (MQ=255)
                                                        cgGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGt                                              <  1:4647008/101‑1 (MQ=255)
                                                         gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                             <  2:2619756/101‑1 (MQ=255)
                                                         gggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGta                                             >  1:582953/1‑101 (MQ=255)
                                                          ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                            <  2:1075932/101‑1 (MQ=255)
                                                          ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                            >  2:1367225/1‑101 (MQ=255)
                                                          ggggCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGtat                                            <  1:2731202/101‑1 (MQ=255)
                                                             gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                         <  1:388947/101‑1 (MQ=255)
                                                             gCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATg                                         <  1:848214/101‑1 (MQ=255)
                                                              cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                        >  1:3305346/1‑101 (MQ=255)
                                                              cGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGc                                        <  2:3346026/101‑1 (MQ=255)
                                                               gTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGcc                                       >  1:2129515/1‑101 (MQ=255)
                                                                 ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                                     >  1:866564/1‑101 (MQ=255)
                                                                 ggATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCg                                     >  1:405523/1‑101 (MQ=255)
                                                                  gATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCgg                                    <  1:308101/101‑1 (MQ=255)
                                                                   aTCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGt                                   <  2:1173173/101‑1 (MQ=255)
                                                                    tCGAAGCTCGCCCTGGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTg                                  <  1:2284985/101‑1 (MQ=255)
                                                                      gacaCTCCCCCTCGATCCGCGCTTGCCTCGCGCCTGC.CGC..GTCTCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                <  1:336498/98‑1 (MQ=38)
                                                                      gAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTgag                                >  1:2657605/1‑101 (MQ=255)
                                                                       aaaCTCGCCCTGGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACAGGTGCATCTGATGATTATGTCCGGTgaga                               >  2:3149029/1‑101 (MQ=37)
                                                                        aaCTCGCCCTCCATCCGCGCTTGCCCCACGCCAGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGAc                              <  1:2504419/101‑1 (MQ=255)
                                                                         ccTCGCCCCCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                             <  1:3703983/100‑1 (MQ=255)
                                                                         aCTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACg                             <  1:1118802/101‑1 (MQ=255)
                                                                             ggccTCGATCCGCGCTTGGCTCCCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTg                         <  2:1440003/99‑1 (MQ=255)
                                                                               ccTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTc                       <  1:1367319/101‑1 (MQ=255)
                                                                                  cGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCt                    <  2:2723795/101‑1 (MQ=255)
                                                                                   gATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTc                   >  2:1061649/1‑101 (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)
                                                                                        gcgcTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTcccccc              >  1:2089826/1‑101 (MQ=255)
                                                                                         cccTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCcg             <  2:1202749/99‑1 (MQ=255)
                                                                                            ttGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCgg          >  2:2903557/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)
                                                                                                  cACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGGGAGCCGGCGGTCGCCCCCCCCGCGGGGCCgg    >  2:4388490/1‑100 (MQ=37)
                                                                                                   aCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTg   >  2:4681798/1‑101 (MQ=255)
                                                                                                   aCGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTg   <  2:2367257/101‑1 (MQ=255)
                                                                                                    |                                                                                                       
GATCGAAACGATGCGTTCATGGCCTACGTCGGGCAGTGGACTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGCCCTCGATCCGCGCTTGCCTCACGCCTGC.CGC..GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAGGCGCATCTGATGTATATGCCCGGTGAGACGGCTGTCGCTCCCCCCGCGGGGGCGTGG  >  NC_014034/1395922‑1396122

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