Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,483 G→A 17.6% I708I (ATC→CGA RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4830G→A17.6% 199.3 ‑4.0 ‑0.1410I708I (ATC→ATTRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (21/51):  ref base (151/187):  total (172/238)
Fisher's exact test strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.91e-01

Only 100 of 434 total aligned reads displayed.

GGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCG  >  NC_014034/681384‑681583
                                                                                                   |                                                                                                      
  ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.c                                                                                                   <  1:791302/101‑1 (MQ=255)
  ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.c                                                                                                   <  2:1012207/101‑1 (MQ=255)
    agcGTGTCCGAGCCCAGCCCGGAGAGCCCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AC.CTg                                                                                                 <  2:1359833/99‑1 (MQ=255)
    aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTg                                                                                                 <  1:8707689/101‑1 (MQ=255)
    aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTg                                                                                                 >  1:3096696/1‑101 (MQ=255)
     tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgc                                                                                                <  1:4743009/101‑1 (MQ=255)
       gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.GTgcgc                                                                                              >  1:7328393/1‑101 (MQ=255)
       gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgcgc                                                                                              >  1:8768215/1‑101 (MQ=255)
       gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgcgc                                                                                              <  2:7832554/101‑1 (MQ=255)
       gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgcgc                                                                                              >  1:8935973/1‑101 (MQ=255)
         ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCg                                                                                            >  2:4837446/1‑101 (MQ=255)
         ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCg                                                                                            <  2:8497720/101‑1 (MQ=255)
          gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGa                                                                                           >  2:7334812/1‑101 (MQ=255)
             gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGcc                                                                                        >  2:3558256/1‑101 (MQ=255)
             gAGCCCAGCCCGGAGAGCGCCGTAGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCACCGGGCCG.AT.CTTGGCCGAGca                                                                                        >  2:2963427/1‑100 (MQ=38)
                cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAg                                                                                     <  1:6293354/101‑1 (MQ=255)
                 ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGa                                                                                    <  2:4325195/101‑1 (MQ=255)
                 ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGa                                                                                    <  1:7957038/101‑1 (MQ=255)
                   agccaGGGGAGCGCCGTGGCGCTGATCTGCGCCCGATCCGCCGCTTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGAcg                                                                                  <  2:6989282/101‑1 (MQ=37)
                       cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGTCCAGACGcgcc                                                                              <  1:5735789/101‑1 (MQ=255)
                        ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGAGCCG.AT.CTGCGCCGAGCCCAGACGcgccc                                                                             >  2:4919892/1‑100 (MQ=255)
                          agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCgccgcc                                                                           <  1:4452826/101‑1 (MQ=255)
                          agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCgccgcc                                                                           >  1:1468991/1‑101 (MQ=255)
                           gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCt                                                                          >  2:908490/1‑101 (MQ=255)
                           gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCt                                                                          >  2:8078913/1‑101 (MQ=255)
                           gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCt                                                                          >  1:8500082/1‑101 (MQ=255)
                            agCGGCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                         <  1:2916031/101‑1 (MQ=255)
                             gcgcCGTCGCGCTGATCTGCGCCAGATCCTCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTgc                                                                        <  1:4920052/101‑1 (MQ=255)
                                 cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                    <  2:1879837/101‑1 (MQ=255)
                                   tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                  <  2:1526210/101‑1 (MQ=255)
                                   tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                  <  1:3834428/101‑1 (MQ=255)
                                    cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                                 >  2:7695604/1‑101 (MQ=255)
                                      cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                               <  2:2080694/101‑1 (MQ=255)
                                         tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                            <  2:1799056/101‑1 (MQ=255)
                                         tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                            <  2:4575151/101‑1 (MQ=255)
                                         tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                            >  1:2064066/1‑101 (MQ=255)
                                         tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                            >  1:2340595/1‑101 (MQ=255)
                                         ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                            <  1:5501695/92‑1 (MQ=2)
                                          gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                           >  1:5826934/1‑101 (MQ=255)
                                           aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                          >  2:5707473/1‑101 (MQ=255)
                                           aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                          >  2:6534962/1‑101 (MQ=255)
                                            gccgccccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                                                         <  2:5579992/95‑1 (MQ=2)
                                             cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                        >  2:7263502/1‑101 (MQ=255)
                                              tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                       >  1:4371361/1‑101 (MQ=255)
                                              cgcaccAGATCGGCGGTTTCGAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGGCAc                                                       <  1:8814632/97‑1 (MQ=12)
                                              ccggcCAGATCCGCCGCCTCGAGCGCCGCGCTCTGCGTGCTGGGCAGCGCGCCG.AT.CTGCGCCGAGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                       <  2:8853564/98‑1 (MQ=21)
                                               gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTAGTCAGCGCCACg                                                      <  1:7476520/101‑1 (MQ=255)
                                               gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                      <  1:5068900/98‑1 (MQ=2)
                                               gcaccAGATCGGCGGGTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                      <  1:2961108/98‑1 (MQ=9)
                                                   cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                  >  2:3976649/1‑101 (MQ=255)
                                                   cAGATCCGCCGCCTCACGCGCCGCGATCTGCTGACTGGATAGCGCGCCG.AT.CTGCGCCGAGCCCTGCTGTGCCGTCTGCGGCGTCGTCTGTGCCACGaaca                                                  >  2:972495/1‑98 (MQ=16)
                                                     gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                >  2:7908291/1‑101 (MQ=255)
                                                      aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                               <  1:2453450/101‑1 (MQ=255)
                                                      aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                               <  1:6381518/101‑1 (MQ=255)
                                                       tCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                              <  1:1011637/101‑1 (MQ=255)
                                                        ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                             >  1:7749541/1‑101 (MQ=255)
                                                         cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                            >  2:3791049/1‑101 (MQ=255)
                                                          gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                                           >  2:7538481/1‑101 (MQ=255)
                                                            ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                         <  2:6155462/94‑1 (MQ=9)
                                                            ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                         >  1:7121834/8‑101 (MQ=9)
                                                            cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                         <  2:8371104/101‑1 (MQ=255)
                                                             gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                        <  2:7095589/95‑1 (MQ=12)
                                                             gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                        >  2:223466/1‑101 (MQ=255)
                                                             gccGCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                        <  1:1074174/101‑1 (MQ=255)
                                                              cccccggcgcCGCGCTCTCCGTGCTGGTCAGCCCGCCG.AT.CTGCGCCGAGCCCAGACCCGCCGCCTGCGCCGTCGTCAGCACCACGATGTCATCGGTGCCa                                       <  1:544292/95‑1 (MQ=25)
                                                               cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                                      >  1:3240182/1‑101 (MQ=255)
                                                                tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                                     <  2:2177314/98‑4 (MQ=9)
                                                                tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                     <  2:8558090/101‑1 (MQ=255)
                                                                tCGAGCGCCGCGAGCTGCGTGCTGGGCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCGGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                     <  2:5276642/101‑1 (MQ=255)
                                                                 ccaGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                    >  2:2414731/3‑97 (MQ=2)
                                                                 ccaGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                    >  1:1770656/3‑97 (MQ=2)
                                                                   aGCGCCGCACCCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                  <  1:705644/101‑1 (MQ=9)
                                                                   aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAG.AAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                  >  1:1259597/1‑101 (MQ=9)
                                                                    gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                                 >  2:3094648/1‑100 (MQ=2)
                                                                      gccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                               >  1:6052557/1‑98 (MQ=2)
                                                                       ccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCgg                              <  1:8355118/101‑1 (MQ=255)
                                                                       ccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                              <  2:8317579/101‑5 (MQ=2)
                                                                        cgggATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCCCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa                             <  2:6096918/98‑1 (MQ=255)
                                                                          cgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                           >  2:4492151/1‑101 (MQ=255)
                                                                           gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                          <  2:2682814/101‑1 (MQ=255)
                                                                            aTCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGGGCCACAGCCGCGGCTccc                         >  1:4285097/1‑99 (MQ=255)
                                                                             tCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                        <  2:3768138/101‑1 (MQ=255)
                                                                              cTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGGGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                       <  1:8137418/101‑1 (MQ=255)
                                                                               tGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccgc                      >  1:2557766/1‑101 (MQ=255)
                                                                                gCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                     <  1:1877876/101‑1 (MQ=255)
                                                                                  gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                   <  1:6867195/101‑1 (MQ=255)
                                                                                  gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCAGCCTGAGCCGTCGTCAGCGCAACGATGTAATCGGTGCCAAAGACGCAGATAGCCGCCGt                   >  1:9098958/1‑101 (MQ=25)
                                                                                   tGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg                  >  2:1215133/1‑101 (MQ=255)
                                                                                     cTGGTCGGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc                <  2:7450661/101‑1 (MQ=255)
                                                                                      tGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc               <  1:3609888/101‑1 (MQ=255)
                                                                                      tGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc               <  2:5960864/101‑1 (MQ=255)
                                                                                       ggTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa              >  1:3935235/1‑101 (MQ=255)
                                                                                       ggTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa              <  2:1986255/101‑1 (MQ=255)
                                                                                         tcaccccgcCG.AT.TTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc            <  2:2909326/95‑1 (MQ=255)
                                                                                          cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg           <  1:3777312/101‑1 (MQ=255)
                                                                                          cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg           >  2:642622/1‑101 (MQ=255)
                                                                                              gcgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccgc       <  1:7566144/101‑1 (MQ=255)
                                                                                               cgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCgccgcc      <  2:509074/101‑1 (MQ=255)
                                                                                                  cG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc   <  2:8187661/101‑1 (MQ=255)
                                                                                                  cG.AT.CGGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc   <  1:5317330/101‑1 (MQ=255)
                                                                                                   |                                                                                                      
GGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCG  >  NC_014034/681384‑681583

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG