Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,660 G→A 7.1% S649S (AGC→GTA)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6600G→A 7.1% 32.0 ‑2.7 ‑0.2210S649S (AGC→AGT) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (11/4):  ref base (93/102):  total (104/106)
Fisher's exact test strand distribution p-value = 6.46e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.86e-01

Only 100 of 216 total aligned reads displayed.

GCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTG  >  NC_014034/681562‑681759
                                                                                                  |                                                                                                   
gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCTGCCACCTCAGGGcg                                                                                                   >  2:2003843/1‑100 (MQ=255)
 ccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTg                                                                                                  >  2:427042/1‑101 (MQ=255)
  cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                                                 <  2:1177050/101‑1 (MQ=255)
   gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                                                                                                <  1:1382821/101‑1 (MQ=255)
   gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                                                                                                <  1:917754/101‑1 (MQ=255)
       ccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgc                                                                                            <  1:403878/101‑1 (MQ=255)
          gcgcCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcg                                                                                         <  2:3550795/101‑1 (MQ=255)
            gccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                                       >  1:408965/1‑101 (MQ=255)
             ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                                                                                      >  1:1058294/1‑101 (MQ=255)
             ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                                                                                      >  2:1103381/1‑101 (MQ=255)
                 cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                                  >  1:126622/1‑101 (MQ=255)
                 cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                                  <  1:27567/101‑1 (MQ=255)
                  aCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTg                                                                                 >  2:1903887/1‑101 (MQ=255)
                  aCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACATg                                                                                 >  2:1370026/1‑101 (MQ=255)
                    cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                               <  1:3880890/101‑1 (MQ=255)
                    cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                               <  1:2216846/101‑1 (MQ=255)
                    cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                               >  2:905999/1‑101 (MQ=255)
                     gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                              >  1:4349562/1‑101 (MQ=255)
                       cGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCgg                                                                            <  2:185091/101‑1 (MQ=255)
                        gCCGAGGCCAGCGCCGCGATCTGCGCCGTCCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCCCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGt                                                                           <  1:636889/101‑1 (MQ=255)
                        gCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGt                                                                           >  1:1652830/1‑101 (MQ=255)
                         cggtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTg                                                                          <  1:4265343/91‑1 (MQ=2)
                         cggtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTg                                                                          <  1:2151101/91‑1 (MQ=2)
                         ccGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTg                                                                          <  2:2831848/101‑1 (MQ=255)
                           gtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGcc                                                                        >  2:1615103/9‑100 (MQ=2)
                           gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCt                                                                        >  2:2799611/1‑101 (MQ=255)
                           gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCAGACCTGATCGGTGCt                                                                        >  2:2997463/1‑101 (MQ=255)
                           gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGAAAGGCTGGAGAGCGCGCCGACCTGACCGGTgtt                                                                        >  2:1183313/1‑99 (MQ=255)
                             aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                      <  1:736875/101‑1 (MQ=255)
                             aGCCAGCGCCCCCATCCCCCCCCTGCTCACCCCCTCCGCCGCCACCGAGCCAAGCCCCGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                      <  2:1727055/101‑1 (MQ=255)
                              gCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAg                                                                     <  1:3590424/101‑1 (MQ=255)
                               ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGc                                                                    >  2:2295135/1‑101 (MQ=255)
                               ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGc                                                                    <  1:996473/101‑1 (MQ=255)
                               caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                    >  2:507259/5‑100 (MQ=2)
                                cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                   >  1:3153217/1‑101 (MQ=255)
                                cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                   >  2:4142150/1‑101 (MQ=255)
                                cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGGGCGCAGCt                                                                   >  2:3633970/1‑101 (MQ=255)
                                aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                   >  1:189717/4‑99 (MQ=2)
                                  gcgcCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTg                                                                 <  2:3942211/101‑1 (MQ=255)
                                  gcgcCGCGATCTGCGCCGCGCCCAGCCCCGGGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactg                                                                 <  1:3101673/101‑5 (MQ=2)
                                   cgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactgg                                                                >  1:1012766/1‑96 (MQ=2)
                                    gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                                                               >  1:373812/1‑95 (MQ=2)
                                    gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                                                               >  1:1764012/1‑95 (MQ=2)
                                     ccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCg                                                              >  1:2689548/1‑101 (MQ=255)
                                     ccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCg                                                              <  1:1370123/101‑1 (MQ=255)
                                        cgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATc                                                           >  2:3424947/1‑101 (MQ=255)
                                          gtCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                         <  2:2621157/100‑1 (MQ=255)
                                          acttGCTCCGTGCTCAGTTTCGCCGCCTGCCCGGGGCCGTGCCCGGCGACCTCAGGTCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTTCCGATCTg                                                         <  2:3031618/98‑1 (MQ=255)
                                          aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                         >  2:2959959/1‑101 (MQ=255)
                                          aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                         >  2:2535764/1‑101 (MQ=255)
                                           tCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGc                                                        <  2:521314/101‑1 (MQ=255)
                                            cTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCg                                                       >  1:3334101/1‑101 (MQ=255)
                                             tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                                                      <  1:3229438/101‑1 (MQ=255)
                                               cgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcg                                                    <  1:2520702/101‑1 (MQ=255)
                                                gcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgc                                                   >  2:3154001/1‑101 (MQ=255)
                                                gcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgc                                                   <  1:4142411/101‑1 (MQ=255)
                                                  cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCGCAGCTTGCCGATCTGCGTCgcgc                                                 >  1:4374512/1‑101 (MQ=255)
                                                     gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa                                              <  2:2678312/101‑1 (MQ=255)
                                                      cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                             <  2:206391/101‑1 (MQ=255)
                                                      cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                             <  2:3779298/101‑1 (MQ=255)
                                                           cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                        <  2:126601/101‑1 (MQ=255)
                                                           cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                        <  1:486286/101‑1 (MQ=255)
                                                              cgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGaaa                                     <  1:1855188/101‑1 (MQ=255)
                                                                ccgccGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                   <  2:3054181/101‑1 (MQ=255)
                                                                 cgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccg                                  >  2:2381945/1‑101 (MQ=255)
                                                                  gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                 >  2:2162072/1‑101 (MQ=255)
                                                                   ccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgcc                                <  1:2024220/101‑1 (MQ=255)
                                                                    cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccg                               >  2:2238031/1‑101 (MQ=255)
                                                                    cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccg                               <  2:578595/101‑1 (MQ=255)
                                                                    cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccc                               >  1:550393/1‑100 (MQ=255)
                                                                      gCACCGAGCTAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCTGCTTGCCGATCTGCTTCGCGCTCAGGCCCGAAACCGCcgcg                             >  2:4421338/1‑101 (MQ=255)
                                                                      gCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCcgcg                             >  1:3248747/1‑101 (MQ=255)
                                                                        actgAGCCAAGCCCGGCCTCCGTAAGGCTGGAGAGCGCGCCGACCTGATCGGTACTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                           <  2:758709/98‑1 (MQ=255)
                                                                        aCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                           >  2:3063797/1‑101 (MQ=255)
                                                                        aCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                           <  2:408925/101‑1 (MQ=255)
                                                                         ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCTCCGCGCTg                          <  1:2960180/101‑1 (MQ=255)
                                                                         ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                          <  1:906084/101‑1 (MQ=255)
                                                                         ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                          <  1:716615/101‑1 (MQ=255)
                                                                             gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                      >  1:4444587/1‑101 (MQ=255)
                                                                              ccAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAg                     <  1:1966985/101‑1 (MQ=255)
                                                                               cAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGc                    <  1:2348044/101‑1 (MQ=255)
                                                                                aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc                   <  2:2971293/101‑1 (MQ=255)
                                                                                aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc                   >  1:3403185/1‑101 (MQ=255)
                                                                                  gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                 <  2:4339725/101‑1 (MQ=255)
                                                                                      ggCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTGAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGa             >  2:2276898/1‑101 (MQ=255)
                                                                                      ggCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGa             <  1:4029455/101‑1 (MQ=255)
                                                                                       gCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGAt            >  2:3143176/1‑101 (MQ=255)
                                                                                         cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt          <  1:3634223/101‑1 (MQ=255)
                                                                                         cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGCTGGTCAGCTTGCCGATCTGCGTTGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt          <  1:3808616/101‑1 (MQ=255)
                                                                                          aCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTg         >  2:2707788/1‑101 (MQ=255)
                                                                                          aCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTg         >  1:2913446/1‑101 (MQ=255)
                                                                                          aCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTg         <  2:654525/101‑1 (MQ=255)
                                                                                           ccGCAAGGCTGGCGAGCGGGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgc        <  2:1458253/101‑1 (MQ=255)
                                                                                           ccGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgc        >  1:4226784/1‑101 (MQ=255)
                                                                                            cGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgcg       <  1:2128271/101‑1 (MQ=255)
                                                                                             gCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgcgc      >  2:3264777/1‑101 (MQ=255)
                                                                                             gCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgcgc      >  1:1000559/1‑101 (MQ=255)
                                                                                              cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc     <  1:1150259/101‑1 (MQ=255)
                                                                                              cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc     <  1:1541733/101‑1 (MQ=255)
                                                                                               aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCAAGCCCCGCGATCTGCGCCg    >  2:1563440/1‑101 (MQ=255)
                                                                                                  |                                                                                                   
GCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTG  >  NC_014034/681562‑681759

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