Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,481 C→A 18.2% G709V (GGC→GTA)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4810C→A18.2% 195.9 ‑4.4 ‑0.3395G709V (GGC→GTC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (21/51):  ref base (148/175):  total (169/226)
Fisher's exact test strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.31e-01

Only 100 of 425 total aligned reads displayed.

CCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAC  >  NC_014034/681382‑681581
                                                                                                   |                                                                                                      
  ggCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCAGCGATCTGCGTGCTGGTCAGCGCGCCG.a                                                                                                    >  2:7034659/1‑101 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                                <  2:8467935/101‑1 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                                <  1:803676/101‑1 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                                >  1:6115362/1‑101 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                                <  1:5374175/101‑1 (MQ=255)
      aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTg                                                                                               >  2:6105077/1‑101 (MQ=255)
       tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgc                                                                                              <  1:8739075/101‑1 (MQ=255)
        cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgcg                                                                                             >  2:1636318/1‑101 (MQ=255)
         gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgcgc                                                                                            >  2:2745817/1‑101 (MQ=255)
         gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGANNNNNGTGCTGGTCAGCGCGCCG.AT.CTgcgc                                                                                            >  1:4628099/1‑101 (MQ=255)
          tGGCCGAGCCCAGCCCGGAGCGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGcc                                                                                           <  1:8820805/101‑1 (MQ=255)
            gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGa                                                                                         >  1:8371467/1‑101 (MQ=255)
               gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGcc                                                                                      <  1:780080/101‑1 (MQ=255)
               gAGCCCAGCCCGGAGAGCGCCGTAGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCACCGGGCCG.AT.CTTGGCCGAGca                                                                                      >  2:2963427/1‑100 (MQ=38)
                agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGccc                                                                                     >  2:8137069/1‑101 (MQ=255)
                 gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCa                                                                                    <  2:3145642/101‑1 (MQ=255)
                 gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCa                                                                                    <  2:5813213/101‑1 (MQ=255)
                    cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCTCGCCG.AT.CTGCGCCGAGCCCAGAc                                                                                 <  2:4382498/101‑1 (MQ=255)
                    cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGAc                                                                                 <  1:7951066/101‑1 (MQ=255)
                       gcgggAGAGCGCAGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGAcgcg                                                                              <  2:1022608/98‑1 (MQ=255)
                         cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgcc                                                                            >  1:4181556/1‑101 (MQ=255)
                         cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgcc                                                                            >  2:2264079/1‑101 (MQ=255)
                          ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgccg                                                                           <  1:8563801/101‑1 (MQ=255)
                           gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgccgc                                                                          >  1:5779110/1‑101 (MQ=255)
                              agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                       >  2:138578/1‑101 (MQ=255)
                              agCGCCGTCGCGCTGATCTGCGCCAGACCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.GTGCGCCGTGCCCAGACGCGCCGCCTg                                                                       <  2:2518776/101‑1 (MQ=255)
                               gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTgc                                                                      <  2:5484785/101‑1 (MQ=255)
                                   cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                  <  2:1879837/101‑1 (MQ=255)
                                      cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                               <  1:7093524/101‑1 (MQ=255)
                                       gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcgt                                                              >  2:7767872/1‑101 (MQ=255)
                                       gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcgc                                                              >  2:9027387/1‑100 (MQ=255)
                                        cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                             <  2:2080694/101‑1 (MQ=255)
                                         gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                            >  1:1413430/1‑101 (MQ=255)
                                         gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                            >  1:633313/1‑101 (MQ=255)
                                         gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTAGTAAGCGCGCCG.AT.CTTCACCGAGCCCAGACGCGCCGCCTGCGCCgtcgtcg                                                            >  1:7355287/1‑100 (MQ=25)
                                          cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAg                                                           <  1:4334396/101‑1 (MQ=255)
                                           tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                          <  2:1799056/101‑1 (MQ=255)
                                           ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                          <  2:5705442/92‑1 (MQ=2)
                                            gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                         >  2:6877291/1‑101 (MQ=255)
                                              tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                       <  1:204451/101‑1 (MQ=255)
                                               ccgcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                      <  2:4549563/99‑1 (MQ=255)
                                                cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                     <  2:7170955/97‑1 (MQ=2)
                                                cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                     <  1:7128055/97‑1 (MQ=2)
                                                ccggcCAGATCCGCCGCCTCGAGCGCCGCGCTCTGCGTGCTGGGCAGCGCGCCG.AT.CTGCGCCGAGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                     <  2:8853564/98‑1 (MQ=21)
                                                   cccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                  <  1:3679598/100‑1 (MQ=255)
                                                    ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                 <  2:7280903/101‑1 (MQ=255)
                                                    ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                 >  1:7486580/1‑101 (MQ=255)
                                                     cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                >  1:5166113/1‑101 (MQ=255)
                                                     cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                >  1:4947848/1‑101 (MQ=255)
                                                     cAGATCCGCCGCCTCACGCGCCGCGATCTGCTGACTGGATAGCGCGCCG.AT.CTGCGCCGAGCCCTGCTGTGCCGTCTGCGGCGTCGTCTGTGCCACGaaca                                                >  2:972495/1‑98 (MQ=16)
                                                       gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                              <  2:4627964/101‑1 (MQ=255)
                                                       gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                              >  2:3798358/1‑101 (MQ=255)
                                                        aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                             <  1:2453450/101‑1 (MQ=255)
                                                         tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                            <  1:3091636/101‑1 (MQ=255)
                                                         tCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                            <  1:1011637/101‑1 (MQ=255)
                                                          ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                           <  2:6049280/101‑1 (MQ=255)
                                                          ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                           >  1:9171611/1‑101 (MQ=255)
                                                           tgccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                          <  2:767235/100‑1 (MQ=255)
                                                           ggcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                          <  2:2177711/91‑1 (MQ=9)
                                                            gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgga                                         >  1:767248/1‑100 (MQ=255)
                                                              cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                       <  2:8371104/101‑1 (MQ=255)
                                                               gccGCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                      <  1:1074174/101‑1 (MQ=255)
                                                                tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                                     <  1:5693567/96‑2 (MQ=9)
                                                                  tCGAGCGCCGCGCTCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                   <  1:2951491/101‑1 (MQ=255)
                                                                  tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                   >  1:3537636/1‑101 (MQ=255)
                                                                  tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                   <  2:8558090/101‑1 (MQ=255)
                                                                   ccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                  <  2:5174006/99‑5 (MQ=9)
                                                                   ccaGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                  >  1:1770656/3‑97 (MQ=2)
                                                                    gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                 <  1:4425093/101‑1 (MQ=255)
                                                                     aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                <  1:6956006/101‑1 (MQ=9)
                                                                     aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAG.AAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                >  1:1259597/1‑101 (MQ=9)
                                                                      gccgcgtgATCTGCGTGCCGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcg                               <  2:2923858/94‑1 (MQ=255)
                                                                        gccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                             >  2:7461390/1‑101 (MQ=255)
                                                                        gccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                             >  2:2015913/1‑101 (MQ=255)
                                                                         ccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCgg                            >  2:637027/1‑101 (MQ=255)
                                                                         ccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCgg                            <  2:5084072/101‑1 (MQ=255)
                                                                         ccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                            <  1:6526339/101‑5 (MQ=2)
                                                                         ccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                            <  2:8317579/101‑5 (MQ=2)
                                                                           gcgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                          <  1:8185859/101‑1 (MQ=255)
                                                                           gcgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                          <  1:2544415/101‑1 (MQ=255)
                                                                           gcgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                          >  2:7974949/1‑101 (MQ=255)
                                                                            cgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                         <  2:6432780/101‑1 (MQ=255)
                                                                               tCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                      <  2:3768138/101‑1 (MQ=255)
                                                                                  gCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                   <  1:1877876/101‑1 (MQ=255)
                                                                                   cGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg                  >  1:9020362/1‑101 (MQ=255)
                                                                                    gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                 <  1:4855742/101‑1 (MQ=255)
                                                                                    gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                 >  2:7370179/1‑101 (MQ=255)
                                                                                    gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                 >  1:2328972/1‑101 (MQ=255)
                                                                                     tGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg                >  2:1215133/1‑101 (MQ=255)
                                                                                      gCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGCTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc               >  1:3197011/1‑101 (MQ=255)
                                                                                       cTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc              >  1:5951844/1‑101 (MQ=255)
                                                                                       cTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCCCCGCCTGCGCCGTCGTCATCGCCACGATATCATCGGTGCCAAACCCACCGAGcaccgcgccgcc              >  2:4257220/1‑95 (MQ=18)
                                                                                        tGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc             >  1:2799621/1‑101 (MQ=255)
                                                                                         ggccAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa            <  2:3574253/98‑1 (MQ=255)
                                                                                         ggTCAGCGCGCCA.AT.CTGCGCAGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa            <  2:4884992/101‑1 (MQ=255)
                                                                                          gTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg           >  2:7877737/1‑101 (MQ=255)
                                                                                          gTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg           >  1:3874839/1‑101 (MQ=255)
                                                                                            cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg         >  1:8788796/1‑101 (MQ=255)
                                                                                            cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg         <  1:3786879/101‑1 (MQ=255)
                                                                                              gcgcgcCG.AT.CTGCGCCGAGCCCGGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcggc       <  1:5483337/101‑3 (MQ=255)
                                                                                                   |                                                                                                      
CCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAC  >  NC_014034/681382‑681581

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