Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,651 C→G 7.8% V652V (GTG→GGG RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6510C→G 7.8% 30.1 ‑2.3 ‑0.4192V652V (GTG→GTCRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (11/4):  ref base (85/92):  total (96/96)
Fisher's exact test strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.28e-01

Only 100 of 202 total aligned reads displayed.

CGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCT  >  NC_014034/681551‑681751
                                                                                                    |                                                                                                    
cgcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGcc                                                                                                      >  2:2216707/1‑101 (MQ=255)
        gccgccGTCCCCAGCGCCGCCCCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                                                                                              <  2:2137810/101‑1 (MQ=255)
           gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCTGCCACCTCAGGGcg                                                                                           >  2:2003843/1‑100 (MQ=255)
           gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCt                                                                                           >  2:2269186/1‑101 (MQ=255)
             cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                                         <  2:1177050/101‑1 (MQ=255)
             aaggCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                                         >  2:444987/4‑101 (MQ=255)
              gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                                                                                        <  1:917754/101‑1 (MQ=255)
              gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                                                                                        <  1:1382821/101‑1 (MQ=255)
              gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                                                                                        >  2:3148326/1‑101 (MQ=255)
                  ccAGCGCCGCCCCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgc                                                                                    <  2:4219558/101‑1 (MQ=255)
                  ccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgc                                                                                    <  1:403878/101‑1 (MQ=255)
                       gccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                               >  1:3127803/1‑101 (MQ=255)
                       gccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                               >  2:3808363/1‑101 (MQ=255)
                       gccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                               >  1:408965/1‑101 (MQ=255)
                       cccccccccGTCGCCGAAGCCAGCGCCGCGATCTGCCCCGTCCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                               <  1:2031088/94‑1 (MQ=38)
                        ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                                                                              >  2:1103381/1‑101 (MQ=255)
                           ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc                                                                           >  1:50991/1‑101 (MQ=255)
                            cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                          <  1:27567/101‑1 (MQ=255)
                            cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                          >  1:126622/1‑101 (MQ=255)
                            cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                          <  1:3915824/101‑1 (MQ=255)
                             cccGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTg                                                                         <  2:3542951/100‑1 (MQ=255)
                             aCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTg                                                                         >  2:1903887/1‑101 (MQ=255)
                               cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                       >  2:905999/1‑101 (MQ=255)
                               cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                       <  1:2216846/101‑1 (MQ=255)
                               cGTCGCCGAAGCCAGCGCCGAGCTCTGCGCCGGGCTCAGCGCCGCCGCCTGCACCGGGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                       <  1:315709/101‑1 (MQ=255)
                                gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                      >  1:4349562/1‑101 (MQ=255)
                                gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGACACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                      >  1:3793731/1‑101 (MQ=255)
                                 tCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCg                                                                     >  2:4272759/1‑101 (MQ=255)
                                 tCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTAAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCg                                                                     <  2:485518/101‑1 (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:2151101/91‑1 (MQ=2)
                                    cggtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTg                                                                  >  1:2926653/11‑101 (MQ=2)
                                    cggtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTg                                                                  <  1:4265343/91‑1 (MQ=2)
                                      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                                                            <  1:996473/101‑1 (MQ=255)
                                           cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                           <  2:2434592/101‑1 (MQ=255)
                                           cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                           >  1:4208000/1‑101 (MQ=255)
                                           aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                           >  1:189717/4‑99 (MQ=2)
                                            atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgact                                                          <  2:2505831/99‑4 (MQ=2)
                                            aGCGCCGGGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                          <  2:2840740/101‑1 (MQ=255)
                                             gcgcCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTg                                                         <  2:576832/101‑1 (MQ=255)
                                             gcgcCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTg                                                         <  2:3942211/101‑1 (MQ=255)
                                              cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                                        >  1:1046694/1‑101 (MQ=255)
                                              cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                                        >  1:1539755/1‑101 (MQ=255)
                                               gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                                                       >  1:373812/1‑95 (MQ=2)
                                               gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                                                       >  1:1764012/1‑95 (MQ=2)
                                               gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                                                       >  1:443971/1‑95 (MQ=2)
                                                ccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCg                                                      <  1:1370123/101‑1 (MQ=255)
                                                ccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCg                                                      >  1:3032558/1‑101 (MQ=255)
                                                 cgcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGa                                                     >  2:2362687/1‑101 (MQ=255)
                                                 cgcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGa                                                     >  1:2678484/1‑101 (MQ=255)
                                                    gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                  <  2:3640480/101‑1 (MQ=255)
                                                     acttGCTCCGTGCTCAGTTTCGCCGCCTGCCCGGGGCCGTGCCCGGCGACCTCAGGTCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTTCCGATCTg                                                 <  2:3031618/98‑1 (MQ=255)
                                                     aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                 >  2:2959959/1‑101 (MQ=255)
                                                       cTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCg                                               >  1:3334101/1‑101 (MQ=255)
                                                        tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                                              <  1:3229438/101‑1 (MQ=255)
                                                         gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                             <  1:1003173/101‑1 (MQ=255)
                                                           gcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgc                                           <  1:4142411/101‑1 (MQ=255)
                                                           gcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgc                                           >  2:3154001/1‑101 (MQ=255)
                                                            ccGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgcg                                          >  2:22907/1‑101 (MQ=255)
                                                             cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCgcgc                                         >  2:1509645/1‑101 (MQ=255)
                                                             cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCgcgc                                         >  2:3071680/1‑101 (MQ=255)
                                                             cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCGCAGCTTGCCGATCTGCGTCgcgc                                         >  1:4374512/1‑101 (MQ=255)
                                                                gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa                                      <  2:2678312/101‑1 (MQ=255)
                                                                     gcgcCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGccc                                 <  1:445027/101‑1 (MQ=255)
                                                                     gcgcCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGccc                                 >  2:3303760/1‑101 (MQ=255)
                                                                      cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                <  2:126601/101‑1 (MQ=255)
                                                                       gccgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGa                               >  2:1614982/1‑101 (MQ=255)
                                                                        ccgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGaa                              <  2:2801893/101‑1 (MQ=255)
                                                                         cgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGaaa                             <  1:1855188/101‑1 (MQ=255)
                                                                          cccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAc                            <  2:4349291/100‑1 (MQ=255)
                                                                           ccgccGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                           <  2:3054181/101‑1 (MQ=255)
                                                                            cgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccg                          >  1:2417080/1‑101 (MQ=255)
                                                                             gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                         >  2:2162072/1‑101 (MQ=255)
                                                                             gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                         <  2:2038838/101‑1 (MQ=255)
                                                                               cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccg                       >  2:3602740/1‑101 (MQ=255)
                                                                               cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccc                       >  1:550393/1‑100 (MQ=255)
                                                                                 gCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCcgcg                     >  2:2755164/1‑101 (MQ=255)
                                                                                 gCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCcgcg                     >  1:3248747/1‑101 (MQ=255)
                                                                                   actgAGCCAAGCCCGGCCTCCGTAAGGCTGGAGAGCGCGCCGACCTGATCGGTACTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                   <  2:758709/98‑1 (MQ=255)
                                                                                   aCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                   <  2:408925/101‑1 (MQ=255)
                                                                                    ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                  <  1:906084/101‑1 (MQ=255)
                                                                                    ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                  <  1:716615/101‑1 (MQ=255)
                                                                                       aGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGccc               <  1:1087679/101‑1 (MQ=255)
                                                                                         ccAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAg             <  2:610323/101‑1 (MQ=255)
                                                                                         ccAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAg             <  1:1966985/101‑1 (MQ=255)
                                                                                          cAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGc            >  1:3856095/1‑101 (MQ=255)
                                                                                           aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc           <  2:2971293/101‑1 (MQ=255)
                                                                                           aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc           >  1:3403185/1‑101 (MQ=255)
                                                                                            aGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGccc          >  2:2576306/1‑101 (MQ=255)
                                                                                            aGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGccc          <  2:50978/101‑1 (MQ=255)
                                                                                            aGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGccc          <  1:2293837/101‑1 (MQ=255)
                                                                                                 ggCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGa     <  1:4029455/101‑1 (MQ=255)
                                                                                                    cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt  <  1:3634223/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCT  >  NC_014034/681551‑681751

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