Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,627 C→G 7.4% A660A (GCG→GCC RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6270C→G 7.4% 24.7 ‑3.2 ‑0.2175A660A (GCG→GCCRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (11/2):  ref base (86/76):  total (97/78)
Fisher's exact test strand distribution p-value = 3.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.02e-01

Only 100 of 209 total aligned reads displayed.

CCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCG  >  NC_014034/681527‑681727
                                                                                                    |                                                                                                    
ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                                                                                      <  1:4007466/101‑1 (MQ=255)
    gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                                                                                  >  1:576900/1‑101 (MQ=255)
     aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                                                                                                 <  2:1226531/101‑1 (MQ=255)
     aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                                                                                                 >  1:1458359/1‑101 (MQ=255)
     aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                                                                                                 >  2:2582160/1‑101 (MQ=255)
      tGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGc                                                                                                >  2:3958104/1‑101 (MQ=255)
      tGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGAACCGCCCCCGTCGCCGAAGCCAGCGCCGCGCTCTCCGCCGTGCTCAGCGCCGCCGCCtcc                                                                                                >  1:2751457/1‑99 (MQ=38)
       gTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCa                                                                                               >  2:1747965/1‑101 (MQ=255)
       gTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCa                                                                                               <  1:1784879/101‑1 (MQ=255)
            cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                                                                                          >  2:2520539/1‑101 (MQ=255)
            cGGTCCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                                                                                          <  2:2009030/101‑1 (MQ=255)
             ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                                                                         >  2:1585657/1‑101 (MQ=255)
             ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                                                                         >  2:549522/1‑101 (MQ=255)
             ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                                                                         <  2:3406215/101‑1 (MQ=255)
             ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                                                                         >  2:3590176/1‑101 (MQ=255)
               tGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCa                                                                                       >  2:2176252/1‑101 (MQ=255)
                gCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCaa                                                                                      <  2:4219384/101‑1 (MQ=255)
                 ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                                                                                     >  2:27554/1‑101 (MQ=255)
                 ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                                                                                     <  1:2309057/101‑1 (MQ=255)
                    aaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGccc                                                                                  >  2:1910086/1‑101 (MQ=255)
                    aaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGccc                                                                                  <  1:3648153/101‑1 (MQ=255)
                     aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg                                                                                 <  1:416315/101‑1 (MQ=255)
                        cgcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGcc                                                                              >  2:2216707/1‑101 (MQ=255)
                         gcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCa                                                                             >  2:461714/1‑101 (MQ=255)
                              tCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCa                                                                        <  1:4379167/101‑1 (MQ=255)
                                gccgccGTCCCCAGCGCCGCCCCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                                                                      <  2:2137810/101‑1 (MQ=255)
                                    ccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTg                                                                  >  2:427042/1‑101 (MQ=255)
                                     cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                 <  2:1177050/101‑1 (MQ=255)
                                     cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                 <  2:71673/101‑1 (MQ=255)
                                     aaggCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                 >  2:444987/4‑101 (MQ=255)
                                      gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                                                                <  1:917754/101‑1 (MQ=255)
                                          ccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgc                                                            <  1:403878/101‑1 (MQ=255)
                                               gccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                       >  1:408965/1‑101 (MQ=255)
                                                ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                                                      <  2:1977927/101‑1 (MQ=255)
                                                ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                                                      >  2:1103381/1‑101 (MQ=255)
                                                 cgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGa                                                     >  1:2159851/1‑101 (MQ=255)
                                                    cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                  >  1:126622/1‑101 (MQ=255)
                                                    cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                  <  1:3915824/101‑1 (MQ=255)
                                                     cccGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTg                                                 <  2:3542951/100‑1 (MQ=255)
                                                     aCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACATg                                                 >  2:1370026/1‑101 (MQ=255)
                                                      ccGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAACCCGTCCACCGCAAGGCTGGAGAGCGCGCCGACCTGa                                                <  2:2029206/101‑1 (MQ=255)
                                                       cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                               <  1:4356763/101‑1 (MQ=255)
                                                       cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                               <  1:2216846/101‑1 (MQ=255)
                                                       cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                               <  1:1361768/101‑1 (MQ=255)
                                                       cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                               >  2:905999/1‑101 (MQ=255)
                                                        gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGACACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                              >  1:3793731/1‑101 (MQ=255)
                                                           gCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGt                                           >  1:2802099/1‑101 (MQ=255)
                                                           gCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGt                                           >  1:1652830/1‑101 (MQ=255)
                                                            cggtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTg                                          >  1:2926653/11‑101 (MQ=2)
                                                             cGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGc                                         >  2:791643/1‑101 (MQ=255)
                                                              gtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGcc                                        >  2:1615103/9‑100 (MQ=2)
                                                              gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCAGACCTGATCGGTGCt                                        >  2:2997463/1‑101 (MQ=255)
                                                              gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGAAAGGCTGGAGAGCGCGCCGACCTGACCGGTgtt                                        >  2:1183313/1‑99 (MQ=255)
                                                                aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                      <  1:736875/101‑1 (MQ=255)
                                                                 gCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAg                                     <  1:3590424/101‑1 (MQ=255)
                                                                  caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                    >  2:507259/5‑100 (MQ=2)
                                                                   cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                   >  2:4142150/1‑101 (MQ=255)
                                                                   cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                   >  1:4208000/1‑101 (MQ=255)
                                                                   cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                   <  2:2434592/101‑1 (MQ=255)
                                                                   cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                   >  1:3153217/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                                 <  1:1379215/101‑1 (MQ=255)
                                                                     gcgcCGCGATCTGCGCCGCGCCCAGCCCCGGGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactg                                 <  1:3101673/101‑5 (MQ=2)
                                                                      cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                >  1:1539755/1‑101 (MQ=255)
                                                                      cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                >  2:2909260/1‑101 (MQ=255)
                                                                      cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                >  1:1046694/1‑101 (MQ=255)
                                                                       gccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGcc                               <  1:963248/101‑1 (MQ=255)
                                                                       gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                               >  1:373812/1‑95 (MQ=2)
                                                                       gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                               >  1:443971/1‑95 (MQ=2)
                                                                       gccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggc                               >  1:1764012/1‑95 (MQ=2)
                                                                         cgcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGa                             >  2:2362687/1‑101 (MQ=255)
                                                                         cgcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGa                             >  1:2678484/1‑101 (MQ=255)
                                                                            gagcTGCCCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCGCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                          <  1:64894/98‑1 (MQ=255)
                                                                            gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                          <  2:3640480/101‑1 (MQ=255)
                                                                                tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                      <  1:3229438/101‑1 (MQ=255)
                                                                                tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                      >  2:2966976/1‑101 (MQ=255)
                                                                                 gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                     <  1:1003173/101‑1 (MQ=255)
                                                                                 gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                     <  1:4303636/101‑1 (MQ=255)
                                                                                  cgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcg                    <  1:2131335/101‑1 (MQ=255)
                                                                                  cgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcg                    <  1:2520702/101‑1 (MQ=255)
                                                                                   gcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgc                   >  2:3154001/1‑101 (MQ=255)
                                                                                    ccGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgcg                  >  2:22907/1‑101 (MQ=255)
                                                                                     cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCgcgc                 >  2:3071680/1‑101 (MQ=255)
                                                                                     cGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCgcgc                 >  2:1509645/1‑101 (MQ=255)
                                                                                      gggcccAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCt                <  2:3934055/96‑1 (MQ=255)
                                                                                        gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa              <  2:2678312/101‑1 (MQ=255)
                                                                                        gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa              <  1:3380581/101‑1 (MQ=255)
                                                                                         cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg             <  2:3779298/101‑1 (MQ=255)
                                                                                         atCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg             <  2:100113/100‑1 (MQ=255)
                                                                                            aGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGcc          >  1:3564982/1‑101 (MQ=255)
                                                                                             gcgcCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGccc         <  1:549574/101‑1 (MQ=255)
                                                                                             gcgcCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGccc         >  2:3303760/1‑101 (MQ=255)
                                                                                              cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg        <  1:486286/101‑1 (MQ=255)
                                                                                              cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg        >  2:1855064/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)
                                                                                                   ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc   >  2:942226/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCG  >  NC_014034/681527‑681727

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