Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,482 C→G 21.5% G709R (GGC→CGA)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4820C→G21.5% 112.2 ‑3.8 ‑0.1198G709R (GGC→CGC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (12/30):  ref base (76/77):  total (91/107)
Fisher's exact test strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.81e-01

Only 100 of 209 total aligned reads displayed.

CCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACC  >  NC_014034/681382‑681582
                                                                                                    |                                                                                                     
   gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.At                                                                                                   <  1:1104107/101‑1 (MQ=255)
     cATCGTGGCCGAGCCGCGCCCGGAGGGGGCCGTCGCGCTGATCTGCGCCAGATCCACCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCACCG.ATCt                                                                                                 <  2:2193460/101‑1 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGGTCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                                                 <  2:2961614/101‑1 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                                                 <  2:874939/101‑1 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                                                 <  1:1578034/101‑1 (MQ=255)
     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                                                 <  2:286813/101‑1 (MQ=255)
      aTCGTGGCCTGGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                                                <  2:2215312/101‑1 (MQ=255)
      aTCGTGGCCGAGGCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                                                >  1:2422813/1‑101 (MQ=255)
       tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgc                                                                                               <  2:3183806/101‑1 (MQ=255)
            gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGa                                                                                          <  2:1157123/101‑1 (MQ=255)
            gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGa                                                                                          <  1:1806921/101‑1 (MQ=255)
             ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAg                                                                                         <  1:225089/101‑1 (MQ=255)
             ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAg                                                                                         <  1:1484222/101‑1 (MQ=255)
               gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGcc                                                                                       >  1:206418/1‑101 (MQ=255)
                agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGccc                                                                                      >  2:4395246/1‑101 (MQ=255)
                 gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCa                                                                                     >  2:2238423/1‑101 (MQ=255)
                  cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAg                                                                                    <  1:971956/101‑1 (MQ=255)
                      gcccAGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcgc                                                                                >  2:160311/1‑101 (MQ=255)
                        ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACgcgc                                                                              >  1:2029343/1‑101 (MQ=255)
                         cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgcc                                                                             >  1:3054402/1‑101 (MQ=255)
                            agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCATCGCGCCG.ATCTGCGCCGAGCCCAGACGCgccgcc                                                                          <  1:2657903/101‑1 (MQ=255)
                            agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCgccgcc                                                                          <  1:1329364/101‑1 (MQ=255)
                             gcggtccGTCCCGCTGATCTGGGCAGGCGCCGCCGCATCGATCGCCGCGATATGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                                         <  2:1995220/96‑1 (MQ=21)
                             gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                                         <  1:3746190/101‑1 (MQ=255)
                             gagCGCCGGCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                                         <  1:197067/101‑1 (MQ=255)
                                  ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                                                                    >  2:3057984/1‑101 (MQ=255)
                                  ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                                                                    >  1:1376502/1‑101 (MQ=255)
                                  ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                                                                    <  2:1107803/101‑1 (MQ=255)
                                   cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                   >  1:1664745/1‑101 (MQ=255)
                                      cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                                >  2:2523639/1‑101 (MQ=255)
                                      cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                                >  1:1057054/1‑101 (MQ=255)
                                      cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                                                >  1:1595208/1‑101 (MQ=255)
                                        cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                              <  1:2182652/101‑1 (MQ=255)
                                        cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                              <  2:2726811/101‑1 (MQ=255)
                                         gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                             <  1:1386350/101‑1 (MQ=255)
                                           ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                           <  1:738253/92‑1 (MQ=2)
                                           ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                           <  2:4043686/92‑1 (MQ=2)
                                            gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                          >  1:68898/1‑101 (MQ=255)
                                            gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                          <  2:2636672/101‑1 (MQ=255)
                                            gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                          <  2:3263145/101‑1 (MQ=255)
                                              tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                        <  2:153344/101‑1 (MQ=255)
                                              tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                        <  1:2002385/101‑1 (MQ=255)
                                               ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGccc                                                       >  2:3550824/6‑100 (MQ=2)
                                               ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                                                       <  2:229577/96‑1 (MQ=2)
                                                cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                      <  1:4292345/97‑1 (MQ=2)
                                                  cgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                                    >  2:826252/1‑101 (MQ=255)
                                                   gcCAGATCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                   <  2:3447877/101‑1 (MQ=255)
                                                   cccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                   <  2:1298262/100‑1 (MQ=2)
                                                   cccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                   <  2:248389/100‑1 (MQ=9)
                                                   accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                   <  2:704030/100‑1 (MQ=9)
                                                    gcAGTTCCGCCGCCCCGAGCGCCGCGATCTGCGCGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGCCGCGCCGCCTGCGCCGCCGTCAGCGCCACGATg                                                  <  2:1471628/100‑1 (MQ=21)
                                                    ccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                  <  1:2497133/101‑1 (MQ=2)
                                                    ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                  <  1:570036/101‑1 (MQ=255)
                                                    ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                  >  1:2843771/1‑101 (MQ=255)
                                                     cAGATCCGCCGCCTAGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                 <  1:205131/101‑1 (MQ=255)
                                                          ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                            <  1:1628757/101‑1 (MQ=255)
                                                           ggcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                           >  1:3009061/11‑101 (MQ=9)
                                                            gcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                                          <  1:3699636/92‑1 (MQ=9)
                                                            gcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                                          <  1:2636600/92‑1 (MQ=9)
                                                            gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                                          <  1:3058205/101‑1 (MQ=255)
                                                            gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                                          >  2:4058028/1‑101 (MQ=255)
                                                              ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                        <  1:3794663/94‑1 (MQ=9)
                                                              cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGGGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                        >  2:1867025/1‑101 (MQ=255)
                                                               gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                       >  2:1377146/7‑101 (MQ=9)
                                                               gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                       <  2:2042354/101‑1 (MQ=255)
                                                               gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                       <  1:1151264/101‑1 (MQ=255)
                                                                  tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                    <  2:1263787/101‑1 (MQ=255)
                                                                  tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                    <  2:4280718/101‑1 (MQ=255)
                                                                  tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                    >  1:940213/1‑101 (MQ=255)
                                                                   cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                   <  1:780467/101‑1 (MQ=255)
                                                                   cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                   >  2:175137/1‑101 (MQ=255)
                                                                    gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                  >  2:2287087/1‑101 (MQ=255)
                                                                    caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                  <  2:3630541/100‑6 (MQ=9)
                                                                     aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                 >  1:3504570/1‑101 (MQ=2)
                                                                      gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                                <  1:2208014/101‑2 (MQ=9)
                                                                       cgccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgc                               >  2:3915570/1‑101 (MQ=255)
                                                                       cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                               >  1:2597334/1‑99 (MQ=2)
                                                                        gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                              <  1:3082810/101‑4 (MQ=9)
                                                                        gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                              <  2:798942/101‑4 (MQ=9)
                                                                          cgcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa                            >  2:4303367/1‑101 (MQ=255)
                                                                           gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                           >  1:4170241/1‑101 (MQ=255)
                                                                           gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                           <  1:1917335/101‑1 (MQ=255)
                                                                             gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                         <  2:1585465/101‑1 (MQ=255)
                                                                             gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                         <  1:3220519/101‑1 (MQ=255)
                                                                               tCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                       >  1:3190304/1‑101 (MQ=255)
                                                                                    gTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccgg                  >  2:2128132/1‑100 (MQ=255)
                                                                                     tGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg                 >  2:4161584/1‑101 (MQ=255)
                                                                                      gCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                <  1:4405694/101‑1 (MQ=255)
                                                                                      gCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                >  2:3647900/1‑101 (MQ=255)
                                                                                         ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCATAGCCGCGGATCGCCGCCGTGCCCa             <  2:2563403/101‑1 (MQ=255)
                                                                                         ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa             <  2:3574863/101‑1 (MQ=255)
                                                                                         ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa             >  1:4339994/1‑101 (MQ=255)
                                                                                          atctcggtgcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg            >  2:2030951/9‑101 (MQ=255)
                                                                                           tCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc           <  1:2101243/101‑1 (MQ=255)
                                                                                           tCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc           <  1:685771/101‑1 (MQ=255)
                                                                                             aGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcgc         >  2:2347894/1‑101 (MQ=255)
                                                                                               cgcgcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccg       <  1:4016881/101‑1 (MQ=255)
                                                                                                 cgcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCgccccc     >  2:636821/1‑101 (MQ=255)
                                                                                                   ccG.ATCTGCCCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAc   <  1:826334/101‑1 (MQ=255)
                                                                                                    cG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc  >  1:1226624/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
CCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACC  >  NC_014034/681382‑681582

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