Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,502 C→G 21.6% R702P (CGT→CCT)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5020C→G21.6% 172.3 ‑1.8 ‑0.6305R702P (CGT→CCT) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (29/37):  ref base (129/110):  total (158/147)
Fisher's exact test strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.50e-01

Only 100 of 324 total aligned reads displayed.

CAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGC  >  NC_014034/681402‑681602
                                                                                                      |                                                                                                    
cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGAc                                                                                                      <  1:6295212/101‑1 (MQ=255)
     cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgcc                                                                                                 <  1:2626774/101‑1 (MQ=255)
      ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgccg                                                                                                >  2:5052794/1‑101 (MQ=255)
       gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgccgc                                                                                               >  1:732015/1‑101 (MQ=255)
       gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgccgc                                                                                               >  2:1009445/1‑101 (MQ=255)
        agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCgccgcc                                                                                              >  1:6106615/1‑101 (MQ=255)
          agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                                            <  1:3795200/101‑1 (MQ=255)
          agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                                            <  1:4129236/101‑1 (MQ=255)
          agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                                            <  2:3736554/101‑1 (MQ=255)
             gcCGTCGCGCTGATCTGCGCCGGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTgcgc                                                                                         <  2:5858508/101‑1 (MQ=255)
              ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                                                                                        <  1:3138634/101‑1 (MQ=255)
               cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                                       >  1:1503299/1‑101 (MQ=255)
                 gcgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                                     <  2:542454/100‑1 (MQ=255)
                       tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                                               >  2:4119247/1‑101 (MQ=255)
                       tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                                               >  2:5329462/1‑101 (MQ=255)
                       ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                                               <  2:4450011/92‑1 (MQ=2)
                       ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                                               <  1:5921145/92‑1 (MQ=2)
                        gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                                              >  1:5960451/1‑101 (MQ=255)
                          tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                                            >  2:4017278/1‑101 (MQ=255)
                          gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                                                                            <  1:1785488/95‑1 (MQ=2)
                           ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                           <  2:2200242/96‑1 (MQ=2)
                           cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                           <  1:3963508/101‑1 (MQ=255)
                           cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                           <  2:2086076/101‑1 (MQ=255)
                            tGCGCCAGACCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                          <  1:6584598/101‑1 (MQ=255)
                            cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                          <  2:5526735/97‑1 (MQ=2)
                            cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                          <  1:231204/97‑1 (MQ=2)
                              caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                                                        >  1:1105011/3‑101 (MQ=2)
                               cccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                                       <  2:1850332/100‑1 (MQ=9)
                               accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                                       <  2:1744351/100‑1 (MQ=9)
                               accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                                       <  1:939404/100‑1 (MQ=9)
                                ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGGCAGCGCCACGATg                                                                      >  2:7057178/1‑101 (MQ=255)
                                  aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGATGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                                                    >  1:3713877/1‑101 (MQ=255)
                                   gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                   >  2:6004962/1‑101 (MQ=255)
                                   gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                   >  1:4411883/1‑101 (MQ=255)
                                   gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                   >  1:3313505/1‑101 (MQ=255)
                                    aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCCCGATGTCAt                                                                  >  2:4343735/1‑101 (MQ=255)
                                     tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                                 <  1:7017639/101‑1 (MQ=255)
                                     tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                                 <  1:6444442/101‑1 (MQ=255)
                                     tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                                 <  1:3596932/101‑1 (MQ=255)
                                       cgcggcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                                               <  2:5009225/101‑1 (MQ=255)
                                       cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                                               >  1:2695396/1‑101 (MQ=255)
                                       cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                                               >  1:1704400/1‑101 (MQ=255)
                                        gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCGTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                                                              >  1:4740623/1‑101 (MQ=255)
                                         cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                                                             >  2:3578396/9‑101 (MQ=9)
                                          ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                            <  2:4929670/94‑1 (MQ=9)
                                          cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                            >  1:3030099/1‑101 (MQ=255)
                                           gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                                           <  1:5462613/95‑1 (MQ=9)
                                            tttccaGCGCCGCAGCCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                                                          <  2:1109432/96‑2 (MQ=9)
                                            tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                                                          <  1:6054484/96‑2 (MQ=9)
                                            ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                                                          >  2:2987071/1‑101 (MQ=255)
                                            ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                                                          >  1:2363580/1‑101 (MQ=255)
                                             ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                                                         >  2:2714753/5‑99 (MQ=2)
                                              tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                                                        <  2:4850458/98‑4 (MQ=9)
                                              tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                        >  2:2378513/1‑101 (MQ=255)
                                               cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                                       <  2:3254667/101‑1 (MQ=255)
                                                gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                                      >  2:6453958/1‑101 (MQ=255)
                                                gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                                      <  1:6724316/101‑1 (MQ=255)
                                                caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                                      >  1:1293356/2‑96 (MQ=9)
                                                caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                                      <  2:6281244/100‑6 (MQ=9)
                                                caGCGCCGCAACCTGCGTGCTCGACAGCGC.CAG.AAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                                      >  1:3190736/2‑96 (MQ=2)
                                                 aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGGCAGCGCCACGCTGTCATCGGCGCCAAAGcc                                                     >  1:556632/1‑101 (MQ=255)
                                                 aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                                     >  1:3492027/1‑101 (MQ=2)
                                                   cgccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgc                                                   >  2:5596131/1‑101 (MQ=255)
                                                   cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                                   <  1:5089339/101‑3 (MQ=9)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                  >  1:1261454/1‑98 (MQ=2)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                  <  1:2233930/101‑4 (MQ=2)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                  <  1:2847488/101‑4 (MQ=9)
                                                     ccgcAACCTGCGTGCTCGCCAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                                                 <  2:6839919/101‑5 (MQ=2)
                                                     ccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                                                 <  1:68503/101‑5 (MQ=2)
                                                        cgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                                              >  1:4355056/1‑101 (MQ=255)
                                                        cgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                                              >  2:1395423/1‑101 (MQ=255)
                                                         tgtCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                             <  2:6922615/99‑1 (MQ=255)
                                                         gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                             >  2:4611648/1‑101 (MQ=255)
                                                         gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                             >  1:4946297/1‑101 (MQ=255)
                                                         gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                             >  1:4527112/1‑101 (MQ=255)
                                                           tCTGCGTGCGGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                                           <  1:5698647/101‑1 (MQ=255)
                                                             tGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccgc                                         <  2:3443149/101‑1 (MQ=255)
                                                              gCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGATCCCATACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCTCGGATCgccgcc                                        <  1:6329794/101‑1 (MQ=255)
                                                                gggCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                                      <  1:131929/99‑1 (MQ=255)
                                                                gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                                      >  1:2746274/1‑101 (MQ=255)
                                                                   cTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCNNNcn                                   >  1:6741340/1‑100 (MQ=255)
                                                                   cTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc                                   >  1:1834766/1‑101 (MQ=255)
                                                                      gTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg                                <  2:1503239/101‑1 (MQ=255)
                                                                      gTCAGCCCCCCG.AC.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg                                <  2:2893659/101‑1 (MQ=255)
                                                                          gcgcgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgcc                            <  1:2498362/101‑1 (MQ=255)
                                                                          gcgcgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgcc                            <  2:1972002/101‑1 (MQ=255)
                                                                             cgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCgccgcc                         <  2:5744931/101‑1 (MQ=255)
                                                                              gcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                        <  2:632945/101‑1 (MQ=255)
                                                                              gcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                        >  2:2575603/1‑101 (MQ=255)
                                                                                cG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc                      <  1:4003680/101‑1 (MQ=255)
                                                                                 g.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg                     >  2:5191397/1‑101 (MQ=255)
                                                                                 g.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg                     <  2:3313389/101‑1 (MQ=255)
                                                                                 g.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg                     >  2:122567/1‑101 (MQ=255)
                                                                                      cTGCGCCGAGCCCAGACGCGCCGCCTGTGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCg                  >  2:338960/1‑101 (MQ=255)
                                                                                              aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc          >  2:3831990/1‑101 (MQ=255)
                                                                                              aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc          <  2:832142/101‑1 (MQ=255)
                                                                                               gCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCa         <  2:5960147/101‑1 (MQ=255)
                                                                                                cccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAg        >  2:6397489/1‑101 (MQ=255)
                                                                                                 ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc       >  1:5289112/1‑101 (MQ=255)
                                                                                                   aGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgcgc     >  1:2008147/1‑101 (MQ=255)
                                                                                                      |                                                                                                    
CAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGC  >  NC_014034/681402‑681602

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG