Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,473 T→A 20.9% T712S (ACC→TCC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4730T→A20.9% 127.1 ‑3.2 ‑0.3212T712S (ACC→TCC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (16/28):  ref base (91/76):  total (108/104)
Fisher's exact test strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.12e-01

Only 100 of 230 total aligned reads displayed.

CTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGC  >  NC_014034/681374‑681573
                                                                                                   |                                                                                                      
cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                                                       >  2:3110101/1‑101 (MQ=255)
 tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                      <  2:1013222/101‑1 (MQ=255)
   aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                                                    >  1:3685386/1‑101 (MQ=255)
     aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                                                  >  1:4718585/1‑101 (MQ=255)
            cgaTCTTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.Tc                                                                                         <  1:1636810/99‑1 (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)
                  tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.Tcagatcg                                                                                   >  2:1636710/1‑95 (MQ=255)
                   ggCCGAGCCCAGCCCGGAGAGCGTCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCg                                                                                  <  1:811363/101‑1 (MQ=255)
                   ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCg                                                                                  >  2:304737/1‑101 (MQ=255)
                   ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCg                                                                                  >  1:1016974/1‑101 (MQ=255)
                    gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGa                                                                                 <  1:3878016/101‑1 (MQ=255)
                    gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGa                                                                                 >  1:2906909/1‑101 (MQ=255)
                      cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGc                                                                               >  1:1709621/1‑101 (MQ=255)
                      cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGc                                                                               >  2:2485234/1‑101 (MQ=255)
                       gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGcc                                                                              >  1:126515/1‑101 (MQ=255)
                       gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGcc                                                                              >  1:3913319/1‑101 (MQ=255)
                       gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGcc                                                                              >  2:4059230/1‑101 (MQ=255)
                        agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGccc                                                                             <  1:1961160/101‑1 (MQ=255)
                          cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAg                                                                           >  1:4376953/1‑101 (MQ=255)
                                ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACgcgc                                                                     <  1:3848647/101‑1 (MQ=255)
                                   gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGcgccgc                                                                  >  1:1692907/1‑101 (MQ=255)
                                           cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                          <  2:1060608/101‑1 (MQ=255)
                                             tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgcc                                                        <  2:4579189/101‑3 (MQ=255)
                                              cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                       >  2:1517549/1‑101 (MQ=255)
                                              cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                       >  2:4631046/1‑101 (MQ=255)
                                              cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGCGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                       >  1:3789320/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                                                    >  1:4658288/1‑101 (MQ=255)
                                                 gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                    >  1:2616391/1‑101 (MQ=255)
                                                 gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                    >  2:1369879/1‑101 (MQ=255)
                                                 gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                    <  2:3293948/101‑1 (MQ=255)
                                                  gggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAg                                                   <  2:3790946/91‑1 (MQ=2)
                                                   tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                  >  2:359491/1‑101 (MQ=255)
                                                   tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                  >  2:1672570/1‑101 (MQ=255)
                                                   ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                  <  2:3440122/92‑1 (MQ=2)
                                                    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)
                                                     aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                >  2:1740547/1‑101 (MQ=255)
                                                      tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                               <  1:1369973/101‑1 (MQ=255)
                                                         gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                            <  1:19430/98‑1 (MQ=2)
                                                              aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                       >  1:2873626/1‑101 (MQ=255)
                                                               gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                      <  2:3946993/101‑1 (MQ=255)
                                                                  ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                   >  2:2667852/1‑101 (MQ=255)
                                                                   cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                  >  1:1219475/1‑101 (MQ=255)
                                                                     cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                                <  1:2599128/93‑1 (MQ=9)
                                                                     ccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                                >  2:2750991/1‑101 (MQ=255)
                                                                      ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                               <  2:4119962/94‑1 (MQ=9)
                                                                      ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                               <  1:2392327/94‑1 (MQ=9)
                                                                      ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                               >  2:2908168/8‑101 (MQ=9)
                                                                       gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                              >  1:4741115/7‑101 (MQ=9)
                                                                       gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                              >  2:4737815/7‑101 (MQ=9)
                                                                        tttccaGCGCCGCACCCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                             <  2:3333571/96‑2 (MQ=9)
                                                                        tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg                             <  2:3466386/96‑2 (MQ=9)
                                                                        ccacgcgcgcCGCGATATGCGTGGTCGTCAGCGCTCCA.A.TCTGCGCCGAGCCCGGACCCCCCGCCTGGGCCGTCGTCAGCGCCACGATGTCATCGGTTCCa                             <  2:227776/95‑1 (MQ=21)
                                                                        ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                             <  1:2015153/101‑1 (MQ=255)
                                                                         ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                            >  1:1565832/5‑99 (MQ=9)
                                                                         cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                            <  2:1709516/101‑1 (MQ=255)
                                                                          tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                           >  2:2619696/4‑98 (MQ=2)
                                                                          tCGAGCGCCGCGATCTGCGGGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGAAGCCCCCCCTGCGCCGGCCTCGGCGCCACGATGACATCGGTTCCaaa                           >  2:2375335/1‑101 (MQ=25)
                                                                           ccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                          <  1:498510/99‑5 (MQ=9)
                                                                           cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                          <  1:4266240/101‑1 (MQ=255)
                                                                            gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                         <  2:4341185/101‑1 (MQ=255)
                                                                            gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                         >  2:4309115/1‑101 (MQ=255)
                                                                            caGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAA.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                         >  1:954597/2‑96 (MQ=2)
                                                                            caGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAA.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                         >  1:897362/2‑96 (MQ=2)
                                                                            caGCGCCGCAACCGGCGTGCTAGACAGCGC.CAG.AGGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                         <  2:640729/100‑6 (MQ=2)
                                                                             aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                        <  2:4158643/101‑1 (MQ=255)
                                                                             aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                        >  1:4110491/1‑101 (MQ=255)
                                                                              gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                       <  1:2181258/101‑2 (MQ=9)
                                                                                gccgcAACCTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                     <  1:2072614/101‑4 (MQ=2)
                                                                                   gcgATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                  >  2:3148797/1‑101 (MQ=255)
                                                                                   gcgATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                  <  1:2751155/101‑1 (MQ=255)
                                                                                    cgATCTGCGTGGTGGTCGGCGCGCCG.A.TCTGCCTCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                 <  1:4050631/101‑1 (MQ=25)
                                                                                    cgATCTGCGTGCGGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                 <  2:539470/101‑1 (MQ=255)
                                                                                     gATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCGAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                <  2:725484/101‑1 (MQ=255)
                                                                                     gATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                >  1:1175207/1‑101 (MQ=255)
                                                                                     gATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                >  2:1151613/1‑101 (MQ=255)
                                                                                      aTCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc               >  2:2834675/1‑101 (MQ=255)
                                                                                      aTCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc               >  2:4664536/1‑101 (MQ=255)
                                                                                       ccTGCGTGCTCGACAGCGCCAGA.A.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATcgca              >  1:540241/2‑100 (MQ=2)
                                                                                        cTGCGTGCTCGACAGCGC.CAGAA.GCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATcgcat             >  2:1937268/1‑99 (MQ=2)
                                                                                         tGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccgc            <  2:4636776/101‑1 (MQ=255)
                                                                                            gTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt         <  2:4458266/101‑1 (MQ=255)
                                                                                             tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg        >  1:1643301/1‑101 (MQ=255)
                                                                                             tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg        >  1:4438953/1‑101 (MQ=255)
                                                                                             tGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg        >  2:235040/1‑101 (MQ=255)
                                                                                               cTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc      >  1:4320134/1‑101 (MQ=255)
                                                                                                tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc     <  1:2668006/101‑1 (MQ=255)
                                                                                                tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc     <  1:2365085/101‑1 (MQ=255)
                                                                                                tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc     <  2:3685157/101‑1 (MQ=255)
                                                                                                tGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGCGCCAAAGCCGCGGATCGCCGCCGTGccc     >  2:3264522/1‑101 (MQ=255)
                                                                                                 ggTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa    >  2:507516/1‑101 (MQ=255)
                                                                                                 ggTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa    <  1:1740652/101‑1 (MQ=255)
                                                                                                   tCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc  <  1:246737/101‑1 (MQ=255)
                                                                                                   |                                                                                                      
CTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.A.TCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGC  >  NC_014034/681374‑681573

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