Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,471 G→C 21.0% T712T (ACC→ACG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4710G→C21.0% 100.2 ‑3.2 ‑0.2207T712T (ACC→ACG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (16/27):  ref base (89/73):  total (106/101)
Fisher's exact test strand distribution p-value = 4.14e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.29e-01

Only 100 of 229 total aligned reads displayed.

CGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCA  >  NC_014034/681371‑681571
                                                                                                    |                                                                                                      
   cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                                                     >  2:3110101/1‑101 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                    <  2:341932/101‑1 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                    <  2:1013222/101‑1 (MQ=255)
     gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                                                   >  1:2196410/1‑101 (MQ=255)
      aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                                                  <  1:1960564/101‑1 (MQ=255)
      aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                                                  <  1:4642062/101‑1 (MQ=255)
        aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                                                >  1:4718585/1‑101 (MQ=255)
              gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.t                                                                                        <  1:2003158/101‑1 (MQ=255)
               cgaTCTTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.Tc                                                                                       <  1:1636810/99‑1 (MQ=255)
               ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.Tc                                                                                       >  1:2815165/1‑101 (MQ=255)
                cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCt                                                                                      <  1:4342142/101‑1 (MQ=255)
                cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCt                                                                                      <  1:4662578/101‑1 (MQ=255)
                 aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTg                                                                                     >  2:3594988/1‑101 (MQ=255)
                 aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTg                                                                                     >  1:1130143/1‑101 (MQ=255)
                  tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.agatcg                                                                                     <  2:4024827/101‑6 (MQ=255)
                  tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTgc                                                                                    <  2:670652/101‑1 (MQ=255)
                    gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTgcgc                                                                                  >  2:4254391/1‑101 (MQ=255)
                    gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTgcgc                                                                                  >  1:3258451/1‑101 (MQ=255)
                      ggCCGAGCCCAGCCCGGAGAGCGTCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCg                                                                                <  1:811363/101‑1 (MQ=255)
                      ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCg                                                                                >  1:1016974/1‑101 (MQ=255)
                       gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGa                                                                               >  1:2906909/1‑101 (MQ=255)
                       gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGa                                                                               <  1:3878016/101‑1 (MQ=255)
                         cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGc                                                                             >  1:1709621/1‑101 (MQ=255)
                         cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGc                                                                             >  1:2726182/1‑101 (MQ=255)
                          gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGcc                                                                            >  1:126515/1‑101 (MQ=255)
                            gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCa                                                                          <  1:1746313/101‑1 (MQ=255)
                             cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAg                                                                         >  1:4376953/1‑101 (MQ=255)
                                  cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGAcgcg                                                                    >  2:3316963/1‑101 (MQ=255)
                                   ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACgcgc                                                                   <  1:3848647/101‑1 (MQ=255)
                                          gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTgc                                                            <  1:901284/101‑1 (MQ=255)
                                              cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                        <  2:1060608/101‑1 (MQ=255)
                                                tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgcc                                                      <  2:4579189/101‑3 (MQ=255)
                                                 cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                     <  1:3983329/101‑1 (MQ=255)
                                                 cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                     >  2:4631046/1‑101 (MQ=255)
                                                   cggtGATTTGCGCCAGTTCCGCCGCCTCGAGCGCCGCTATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgttgtc                                                   <  2:2435431/98‑1 (MQ=21)
                                                   cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                   >  1:2637116/1‑101 (MQ=255)
                                                   cgcTGATCGGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                   <  2:2196286/101‑1 (MQ=255)
                                                    gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                  >  2:1369879/1‑101 (MQ=255)
                                                     gggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAg                                                 <  2:3790946/91‑1 (MQ=2)
                                                      tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                >  2:1672570/1‑101 (MQ=255)
                                                       gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                                               <  2:1264514/93‑1 (MQ=2)
                                                       gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                                               <  1:3384721/93‑1 (MQ=2)
                                                       gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                               >  2:3848803/1‑101 (MQ=255)
                                                         tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                             <  1:1369973/101‑1 (MQ=255)
                                                             caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                         <  2:3302928/99‑1 (MQ=9)
                                                              accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                        <  1:984012/100‑1 (MQ=2)
                                                                cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                      >  1:2452113/1‑101 (MQ=255)
                                                                  gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                    <  2:3946993/101‑1 (MQ=255)
                                                                   aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTAAt                                   >  2:1646819/1‑101 (MQ=255)
                                                                     ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                 >  2:2667852/1‑101 (MQ=255)
                                                                      cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                >  1:1219475/1‑101 (MQ=255)
                                                                        cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                              >  2:3532720/9‑101 (MQ=9)
                                                                         ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                             <  2:4119962/94‑1 (MQ=9)
                                                                         ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                             <  1:2090106/94‑1 (MQ=9)
                                                                         ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                             >  2:2908168/8‑101 (MQ=9)
                                                                          gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                            >  2:4737815/7‑101 (MQ=9)
                                                                          gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                            >  1:4741115/7‑101 (MQ=9)
                                                                           tttccaGCGCCGCAGCCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                           <  2:3093756/96‑2 (MQ=9)
                                                                           tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                           <  2:3466386/96‑2 (MQ=9)
                                                                            ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                          >  2:2309584/5‑99 (MQ=9)
                                                                            ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                          >  1:1565832/5‑99 (MQ=9)
                                                                            cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                          <  2:1709516/101‑1 (MQ=255)
                                                                            cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                          >  1:4282807/1‑101 (MQ=255)
                                                                             tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                         >  2:2619696/4‑98 (MQ=2)
                                                                             tCGAGCGCCGCGATCTGCGGGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGAAGCCCCCCCTGCGCCGGCCTCGGCGCCACGATGACATCGGTTCCaaa                         >  2:2375335/1‑101 (MQ=25)
                                                                             cggatcatggcgATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                         >  2:1176306/10‑101 (MQ=255)
                                                                              cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                        >  2:2768087/1‑101 (MQ=255)
                                                                              cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                        <  1:3833270/101‑1 (MQ=255)
                                                                               gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                       >  2:4309115/1‑101 (MQ=255)
                                                                               gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                       <  2:4341185/101‑1 (MQ=255)
                                                                               caGCGCCGCAACCGGCGTGCTAGACAGCGC.CAG.AGGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                       <  2:640729/100‑6 (MQ=2)
                                                                                aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                      >  1:4110491/1‑101 (MQ=255)
                                                                                aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                      <  2:4158643/101‑1 (MQ=255)
                                                                                  cgccgcAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                    <  1:648886/101‑3 (MQ=9)
                                                                                    ccgcAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                  <  1:1694970/101‑5 (MQ=2)
                                                                                      gcgATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                >  2:3148797/1‑101 (MQ=255)
                                                                                       cgATCTGCGTGCGGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc               <  2:539470/101‑1 (MQ=255)
                                                                                        gATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCGAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg              <  2:725484/101‑1 (MQ=255)
                                                                                        gATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg              <  1:2570581/101‑1 (MQ=255)
                                                                                        gATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg              >  1:1175207/1‑101 (MQ=255)
                                                                                         aTCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc             >  2:4664536/1‑101 (MQ=255)
                                                                                         aTCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc             >  2:2834675/1‑101 (MQ=255)
                                                                                           cTGCGTGCTCGACAGCGC.CAGAA.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATcgcat           >  2:1937268/1‑99 (MQ=2)
                                                                                            tGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccgc          <  2:4636776/101‑1 (MQ=255)
                                                                                             gCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc         >  2:1794634/1‑101 (MQ=255)
                                                                                              cGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg        >  2:4086839/1‑101 (MQ=255)
                                                                                              cGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg        <  1:509934/101‑1 (MQ=255)
                                                                                               gTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt       <  2:4458266/101‑1 (MQ=255)
                                                                                               gTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt       >  2:1555431/1‑101 (MQ=255)
                                                                                                tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGTGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCAAAAGCCGCGGATCGCCGCCGTg      <  1:4059472/101‑1 (MQ=255)
                                                                                                tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg      >  1:1643301/1‑101 (MQ=255)
                                                                                                tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg      >  1:3727914/1‑101 (MQ=255)
                                                                                                tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg      >  2:235040/1‑101 (MQ=255)
                                                                                                   tGGTCAGCGCGCCG.A.TCTGCTCCGAGCCCAGACGCACCGCCTGGGCGGTCGCCAGCGCCACGATGTCATCGGTGGCAAAGCCTTGGATCGCCGCCGTGccc   <  2:250019/101‑1 (MQ=25)
                                                                                                   tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc   <  1:2365085/101‑1 (MQ=255)
                                                                                                   tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc   <  1:2668006/101‑1 (MQ=255)
                                                                                                   tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGCGCCAAAGCCGCGGATCGCCGCCGTGccc   >  2:3264522/1‑101 (MQ=255)
                                                                                                    ggTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa  <  1:2062363/101‑1 (MQ=255)
                                                                                                    ggTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa  >  2:507516/1‑101 (MQ=255)
                                                                                                    ggTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa  >  1:2608189/1‑101 (MQ=255)
                                                                                                    |                                                                                                      
CGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCA  >  NC_014034/681371‑681571

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