Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,673 C→G 11.8% G645A (GGC→GCC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6730C→G11.8% 48.4 ‑2.1 ‑0.0194G645A (GGC→GCC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (16/7):  ref base (86/85):  total (102/92)
Fisher's exact test strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.98e-01

Only 100 of 212 total aligned reads displayed.

CGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAG  >  NC_014034/681573‑681773
                                                                                                    |                                                                                                    
 gccgccCCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                                                     <  2:446730/101‑1 (MQ=255)
 gccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTAAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGcc                                                                                                     >  2:1327279/1‑101 (MQ=255)
  ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                                                                                                    >  2:1819719/1‑101 (MQ=255)
   cgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGa                                                                                                   <  1:1270101/101‑1 (MQ=255)
   cgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGa                                                                                                   <  1:1572991/101‑1 (MQ=255)
    gccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAc                                                                                                  >  1:4571084/1‑101 (MQ=255)
    gccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAc                                                                                                  >  2:3524556/1‑101 (MQ=255)
     ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc                                                                                                 <  1:3242083/101‑1 (MQ=255)
     ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc                                                                                                 >  1:1388146/1‑101 (MQ=255)
      ccccGCCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                                                <  2:2250995/99‑1 (MQ=255)
      cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                                                >  1:382138/1‑101 (MQ=255)
          gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                                            <  1:1023784/101‑1 (MQ=255)
          gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                                            <  2:4173304/101‑1 (MQ=255)
          gTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATc                                                                                            <  1:3502447/101‑1 (MQ=255)
           tCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCg                                                                                           <  2:4553139/101‑1 (MQ=255)
           tCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCg                                                                                           >  2:3171159/1‑101 (MQ=255)
           tCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCg                                                                                           >  2:3657775/1‑101 (MQ=255)
              cggtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTg                                                                                        >  2:2179470/11‑101 (MQ=2)
               gggggcaaccgccgcGTTCTGCCCCGCGCTCAGCGCCGGATTCTGCGCGGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGc                                                                                       <  2:3712235/92‑1 (MQ=255)
               cGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGc                                                                                       >  1:3705223/1‑101 (MQ=255)
                 tggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGccc                                                                                     >  1:3407786/8‑99 (MQ=2)
                    ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGc                                                                                  >  1:2027272/1‑101 (MQ=255)
                    ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGc                                                                                  >  1:4231961/1‑101 (MQ=255)
                    caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                                  >  2:789765/5‑100 (MQ=2)
                     cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                                 >  1:899127/1‑101 (MQ=255)
                     cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                                 >  2:376650/1‑101 (MQ=255)
                     cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                                 <  1:1431692/101‑1 (MQ=255)
                     aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                                 >  1:3651990/4‑99 (MQ=2)
                      atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgact                                                                                >  2:342457/3‑98 (MQ=2)
                      atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgact                                                                                >  1:1116513/3‑98 (MQ=2)
                      atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgact                                                                                >  1:1746847/3‑98 (MQ=2)
                      agggccgcGATCTGCGCCGTGCTCAGCTCCGCCGCCTGCACCGAGCCAAGCCCGGCCCCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                                <  1:2780262/98‑1 (MQ=255)
                      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                                <  2:4081486/101‑1 (MQ=255)
                       gcgcCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTg                                                                               >  2:3236072/1‑101 (MQ=255)
                        cgccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGc                                                                              <  2:2080561/101‑1 (MQ=255)
                        cgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactgg                                                                              >  1:3687893/1‑96 (MQ=2)
                             cgATCTGCGCCGTGCTCAGCGCGGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATc                                                                         <  2:2018014/101‑1 (MQ=255)
                             cgATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggccagc                                                                         >  2:2785727/1‑91 (MQ=2)
                              gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                                        <  2:2332850/101‑1 (MQ=255)
                              gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                                        <  1:4384166/101‑1 (MQ=255)
                               aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                       <  1:2433645/101‑1 (MQ=255)
                               aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                                       >  1:4097415/1‑101 (MQ=255)
                                tCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGc                                                                      <  2:510778/101‑1 (MQ=255)
                                tCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGc                                                                      >  1:1095313/1‑101 (MQ=255)
                                 cTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGccaaatcgga                                                                     <  2:183316/101‑9 (MQ=255)
                                 cTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCg                                                                     <  2:899221/101‑1 (MQ=255)
                                 cTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCg                                                                     >  2:3612449/1‑101 (MQ=255)
                                  tGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCAAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGt                                                                    >  2:2203460/1‑101 (MQ=255)
                                   gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                                   <  2:3936244/101‑1 (MQ=255)
                                   gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                                   >  2:3185466/1‑101 (MQ=255)
                                    cgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcg                                                                  >  1:4461722/1‑101 (MQ=255)
                                    cgcCGGGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcg                                                                  <  1:427694/101‑1 (MQ=255)
                                     gcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgc                                                                 >  1:783063/1‑101 (MQ=255)
                                         tGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTc                                                             <  1:1865930/101‑1 (MQ=255)
                                          gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa                                                            >  2:3474895/1‑101 (MQ=255)
                                           cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGTCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                                           <  1:3977472/101‑1 (MQ=255)
                                           cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                                           >  2:690154/1‑101 (MQ=255)
                                            tCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAGACTGGCCAGCTGCATCGTGCCCAgc                                                          >  2:26012/1‑100 (MQ=0)
                                            tCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAGACTGGCCAGCTGCATCGTGCCCAgc                                                          >  2:1077151/1‑100 (MQ=0)
                                            tCAGCGCCGCCGCCTGCCCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAgg                                                          <  2:1434973/101‑1 (MQ=255)
                                            tCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAgg                                                          >  2:3018374/1‑101 (MQ=255)
                                             cAGCGCCGCCGACCGCACCGAGCAAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGc                                                         <  1:482585/101‑1 (MQ=255)
                                                    gccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAc                                                  >  1:2687503/1‑101 (MQ=255)
                                                     ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                                 <  2:4136256/101‑1 (MQ=255)
                                                      cgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccg                                                >  2:4274378/1‑101 (MQ=255)
                                                        ccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgcc                                              <  1:328537/101‑1 (MQ=255)
                                                         cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGCGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccg                                             >  2:1659954/1‑101 (MQ=255)
                                                          tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                            >  2:2692411/1‑101 (MQ=255)
                                                          tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                            >  2:3813845/1‑101 (MQ=255)
                                                          tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                            >  1:3749042/1‑101 (MQ=255)
                                                            cACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCgcgc                                          >  2:2993956/1‑101 (MQ=255)
                                                             aCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                                         >  1:3209136/1‑101 (MQ=255)
                                                              ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                                        <  1:4421032/101‑1 (MQ=255)
                                                                gAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGcc                                      >  1:963361/1‑101 (MQ=255)
                                                                 aGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGccc                                     <  1:660370/101‑1 (MQ=255)
                                                                  gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                                    <  2:3377573/101‑1 (MQ=255)
                                                                     aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc                                 <  1:3018603/101‑1 (MQ=255)
                                                                     aaGCCAGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc                                 <  1:2318660/101‑1 (MQ=255)
                                                                      aGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGccc                                >  2:249521/1‑101 (MQ=255)
                                                                       gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                               >  1:1590752/1‑101 (MQ=255)
                                                                       gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                               >  1:838888/1‑101 (MQ=255)
                                                                           ggCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGa                           >  2:2917593/1‑101 (MQ=255)
                                                                             ccACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATc                         <  2:2469069/101‑1 (MQ=255)
                                                                              cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt                        >  2:372204/1‑101 (MQ=255)
                                                                                ccGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgc                      <  2:339311/101‑1 (MQ=255)
                                                                                  gCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgcgc                    <  1:2554640/101‑1 (MQ=255)
                                                                                   cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc                   <  1:2136127/101‑1 (MQ=255)
                                                                                   cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc                   <  1:28872/101‑1 (MQ=255)
                                                                                   cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc                   >  1:2153821/1‑101 (MQ=255)
                                                                                    aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCg                  >  2:3998939/1‑101 (MQ=255)
                                                                                     aGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGt                 >  2:870868/1‑101 (MQ=255)
                                                                                     aGGCTGGAGAGCGCGCCGACCTGATCGGCGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGGGCTGCCCAGCCCCGCGATCTGCGCCgg                 >  1:3815859/1‑100 (MQ=255)
                                                                                         tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc             >  2:2877329/1‑101 (MQ=255)
                                                                                         tGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccc             >  1:3133360/1‑101 (MQ=255)
                                                                                           gagaGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAg           >  2:1752501/1‑101 (MQ=255)
                                                                                               gcgcgcCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTc       <  2:4425985/101‑1 (MQ=255)
                                                                                                 gcgcCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGc     >  1:1542045/1‑101 (MQ=255)
                                                                                                  cgcCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCa    <  2:1104519/101‑1 (MQ=255)
                                                                                                    ccGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAg  <  2:1467310/101‑1 (MQ=255)
                                                                                                    ccGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAg  <  2:2092254/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAG  >  NC_014034/681573‑681773

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