Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,480 G→C 21.7% G709G (GGC→GGG RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4800G→C21.7% 176.2 ‑2.7 ‑1.0303G709G (GGC→GGGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (27/38):  ref base (129/105):  total (159/144)
Fisher's exact test strand distribution p-value = 6.77e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.11e-01

Only 100 of 323 total aligned reads displayed.

CCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCA  >  NC_014034/681381‑681580
                                                                                                   |                                                                                                      
cccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgcgc                                                                                                       >  2:2358682/1‑101 (MQ=255)
 ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGcc                                                                                                      <  1:4395944/101‑1 (MQ=255)
  cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                                                                     >  1:175494/1‑101 (MQ=255)
  cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                                                                     >  1:562620/1‑101 (MQ=255)
    gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.At                                                                                                  <  1:2263511/101‑1 (MQ=255)
    gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.At                                                                                                  <  1:4650493/101‑1 (MQ=255)
      cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                               >  1:3939324/1‑101 (MQ=255)
      cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                               <  1:81504/101‑1 (MQ=255)
       aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTg                                                                                              >  2:5894934/1‑101 (MQ=255)
          gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgcgc                                                                                           >  1:60321/1‑101 (MQ=255)
           tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGcc                                                                                          <  2:6977635/101‑1 (MQ=255)
              ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAg                                                                                       >  2:4714688/1‑101 (MQ=255)
              ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAg                                                                                       <  1:4209453/101‑1 (MQ=255)
               cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGc                                                                                      >  1:3396797/1‑101 (MQ=255)
                 agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGccc                                                                                    <  1:7063158/101‑1 (MQ=255)
                 agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGccc                                                                                    >  2:1239321/1‑101 (MQ=255)
                   cccagccgGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAg                                                                                  <  2:4852400/101‑1 (MQ=255)
                   cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAg                                                                                  <  1:1044410/101‑1 (MQ=255)
                       gcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGAcgc                                                                              >  2:6435758/1‑101 (MQ=255)
                             agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCgccgcc                                                                        <  2:6879587/101‑1 (MQ=255)
                               agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                      <  2:3736554/101‑1 (MQ=255)
                                gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTgc                                                                     <  1:6691583/101‑1 (MQ=255)
                                    cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                 >  1:1503299/1‑101 (MQ=255)
                                     gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgt                                                                >  1:1665798/1‑101 (MQ=255)
                                      tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                               <  2:311680/101‑1 (MQ=255)
                                         cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                            >  1:3846750/1‑101 (MQ=255)
                                          gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                           <  2:6743848/101‑1 (MQ=255)
                                          gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                           >  2:5641927/1‑101 (MQ=255)
                                            tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                         <  1:6094901/101‑1 (MQ=255)
                                             gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCGGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                        <  1:1267621/101‑1 (MQ=255)
                                             gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                        >  2:1588922/1‑101 (MQ=255)
                                             cagccgcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                        >  2:5820085/6‑101 (MQ=255)
                                              aTCTGCGCCAGAACCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                       >  1:5079241/1‑101 (MQ=255)
                                               tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                      >  2:4017278/1‑101 (MQ=255)
                                               gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                                                      <  1:256566/95‑1 (MQ=2)
                                               gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                                                      <  1:1785488/95‑1 (MQ=2)
                                                ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                                                     <  2:2200242/96‑1 (MQ=2)
                                                cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                     >  1:6321273/1‑101 (MQ=255)
                                                 tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                    <  2:747290/101‑1 (MQ=255)
                                                 tGCGCCAGACCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                    <  1:6584598/101‑1 (MQ=255)
                                                 cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                    <  1:231204/97‑1 (MQ=2)
                                                  gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCGACg                                                   >  2:2115011/1‑101 (MQ=255)
                                                  gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                   <  2:5024420/98‑1 (MQ=2)
                                                   caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                                  >  1:592505/3‑101 (MQ=2)
                                                    cccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                 <  2:4086836/100‑1 (MQ=9)
                                                    accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                 <  2:1744351/100‑1 (MQ=9)
                                                       aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                              >  1:4105284/1‑101 (MQ=255)
                                                       aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                              >  1:2919747/1‑101 (MQ=255)
                                                       aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                              <  1:2743221/101‑1 (MQ=255)
                                                       aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGATGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                              >  1:3713877/1‑101 (MQ=255)
                                                        gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                             >  1:3313505/1‑101 (MQ=255)
                                                        gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                             >  1:4411883/1‑101 (MQ=255)
                                                         aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCCCGATGTCAt                                            >  2:4343735/1‑101 (MQ=255)
                                                         aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                            <  1:5278396/101‑1 (MQ=255)
                                                          tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                           <  1:3596932/101‑1 (MQ=255)
                                                          tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                           >  2:552712/1‑101 (MQ=255)
                                                          cccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                           <  1:991507/100‑1 (MQ=255)
                                                            ggcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                         >  2:2729381/11‑101 (MQ=9)
                                                            cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                         >  2:616016/1‑101 (MQ=255)
                                                              cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                                       >  1:2628802/9‑101 (MQ=9)
                                                               ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                      <  2:4929670/94‑1 (MQ=9)
                                                               ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                      >  2:2907028/8‑101 (MQ=9)
                                                               cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                      >  1:3030099/1‑101 (MQ=255)
                                                                gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                     <  2:1879065/95‑1 (MQ=12)
                                                                gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                     >  2:3307890/7‑101 (MQ=9)
                                                                gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                     <  1:3981151/95‑1 (MQ=9)
                                                                gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACTCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                     <  2:3394130/101‑1 (MQ=255)
                                                                 tttccaGCGCCGCAGCCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                                    <  2:1109432/96‑2 (MQ=9)
                                                                 tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                                    >  2:5107543/6‑100 (MQ=2)
                                                                 tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                                    >  2:1315635/6‑100 (MQ=2)
                                                                 ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                                    >  2:5983090/1‑101 (MQ=255)
                                                                 ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                                    >  1:2363580/1‑101 (MQ=255)
                                                                  ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                                   >  1:5463988/5‑99 (MQ=2)
                                                                   tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                  <  1:6253503/101‑1 (MQ=255)
                                                                     gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                <  1:2767272/101‑1 (MQ=255)
                                                                     gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                <  1:5895809/101‑1 (MQ=255)
                                                                     caGCGCCGCAACCTGCGTGCTCGACAGCGC.CAG.AAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                >  1:3190736/2‑96 (MQ=2)
                                                                      aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGC.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                               <  2:5259344/101‑1 (MQ=9)
                                                                      aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                               >  2:6387363/1‑101 (MQ=9)
                                                                      aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                               >  2:3525382/1‑101 (MQ=2)
                                                                          ccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                           <  1:68503/101‑5 (MQ=2)
                                                                          ccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                           <  1:6994464/101‑5 (MQ=2)
                                                                              gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                       >  1:6738841/1‑101 (MQ=255)
                                                                              gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                       >  1:4527112/1‑101 (MQ=255)
                                                                              gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                       >  1:4946297/1‑101 (MQ=255)
                                                                               aTCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                      >  2:2449926/1‑101 (MQ=255)
                                                                                tCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                     <  2:3007016/101‑1 (MQ=255)
                                                                                tCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                     <  2:540994/101‑1 (MQ=255)
                                                                                tCTGCGTGCGGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                     <  1:5698647/101‑1 (MQ=255)
                                                                                   gCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                  <  2:1441887/101‑1 (MQ=255)
                                                                                     gggCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                <  1:131929/99‑1 (MQ=255)
                                                                                       gCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc              >  2:3483571/1‑101 (MQ=255)
                                                                                        cTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc             >  1:1834766/1‑101 (MQ=255)
                                                                                          ggTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa           <  1:63946/101‑1 (MQ=255)
                                                                                             cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAggc        >  2:6316033/1‑99 (MQ=255)
                                                                                             cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg        >  2:4130423/1‑101 (MQ=255)
                                                                                             cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg        >  2:4473767/1‑101 (MQ=255)
                                                                                             cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGACAAAGCAGCGGATCGCCGCCGTGccagatc        >  1:1222562/1‑96 (MQ=255)
                                                                                               gcgcgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgcc      <  2:1972002/101‑1 (MQ=255)
                                                                                                   gcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa  >  2:2575603/1‑101 (MQ=255)
                                                                                                   |                                                                                                      
CCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCA  >  NC_014034/681381‑681580

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