Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,447 C→T 8.5% E720E (GAG→ACT)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4470C→T 8.5% 31.9 ‑2.9 ‑0.2179E720E (GAG→GAA) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/12):  ref base (77/86):  total (81/98)
Fisher's exact test strand distribution p-value = 5.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.95e-01

Only 100 of 196 total aligned reads displayed.

GCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCC  >  NC_014034/681347‑681545
                                                                                                    |                                                                                                  
gcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                    <  2:1613421/101‑1 (MQ=255)
gcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                    >  1:2961835/1‑101 (MQ=255)
 cgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCt                                                                                                   <  2:3413588/101‑1 (MQ=255)
  gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                  <  2:637024/101‑1 (MQ=255)
  gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                  >  1:3688793/1‑101 (MQ=255)
   tCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                 <  2:863123/101‑1 (MQ=255)
   gcgccAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                 <  1:586702/97‑1 (MQ=255)
     tcagatctcgcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGGCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCACCGCCTCGAg                                                                                               <  2:727727/93‑1 (MQ=255)
     gACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAg                                                                                               >  2:3536806/1‑101 (MQ=255)
      aCAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgc                                                                                              <  1:221463/101‑1 (MQ=255)
      aCAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgc                                                                                              >  2:1318277/1‑101 (MQ=255)
      aCAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCACTGATCTGCGCCAGATCCGCCGCCTCGAgc                                                                                              >  2:1116833/1‑101 (MQ=255)
       cAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgcg                                                                                             <  1:3746817/101‑1 (MQ=255)
        aGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgcgc                                                                                            >  1:3497743/1‑101 (MQ=255)
          cgccgcGATCCGCGCGCCTAAAAGCCCGCCCATCGGGGCCGAGCCCAGCCCGGAGAGCGCCGGCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccg                                                                                          <  2:205135/101‑1 (MQ=255)
           gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                         <  1:3900821/101‑1 (MQ=255)
           gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                         <  2:4049730/101‑1 (MQ=255)
             cgcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgagatcgga                                                                                       >  1:727807/1‑93 (MQ=255)
              gcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGAt                                                                                      <  2:2338672/101‑1 (MQ=255)
              gcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGAt                                                                                      <  1:377909/101‑1 (MQ=255)
               cgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATc                                                                                     >  2:1077753/1‑101 (MQ=255)
                gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccagatcgg                                                                                    >  1:429427/1‑94 (MQ=255)
                gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                    <  1:2171344/101‑1 (MQ=255)
                gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                    <  2:1919085/101‑1 (MQ=255)
                gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                    >  1:4115150/1‑101 (MQ=255)
                 aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                   >  1:2126091/1‑101 (MQ=255)
                   cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                 >  1:4219819/1‑101 (MQ=255)
                   cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                 <  2:4227612/101‑1 (MQ=255)
                    tGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGt                                                                                >  2:2182513/1‑101 (MQ=255)
                        cgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCTGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGcca                                                                            <  2:3133296/101‑3 (MQ=255)
                        cgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTg                                                                            >  1:941916/1‑101 (MQ=255)
                        cgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTg                                                                            >  2:789553/1‑101 (MQ=255)
                           cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                         <  1:789635/101‑1 (MQ=255)
                           cTGAAAGCCCGGCCAGCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                         <  2:4114890/101‑1 (MQ=255)
                            tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                        >  1:1736606/1‑101 (MQ=255)
                            tGAAAGCCCGGCCATCGGGGCCGAGCCCAGCCCGGAGAGCGCCGGCGCGCTGATCTGCGCCAGATCCGCCGCGTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                        >  2:2657731/1‑101 (MQ=255)
                                 gCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgcg                                                                   >  2:1644431/1‑101 (MQ=255)
                                   ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGcc                                                                 >  1:219222/1‑101 (MQ=255)
                                     ggCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGa                                                               >  2:1849205/1‑101 (MQ=255)
                                      gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                              <  1:1104107/101‑1 (MQ=255)
                                       ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                             >  2:2433275/1‑101 (MQ=255)
                                       cactCGTGCCCGAGCCCAGCCCGGAGCGCGCCCTCGCCCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCCCGCCGATc                                                             <  2:2912958/98‑1 (MQ=37)
                                        cATCGTGGCCGAGCCGCGCCCGGAGGGGGCCGTCGCGCTGATCTGCGCCAGATCCACCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCACCGATCt                                                            <  2:2193460/101‑1 (MQ=255)
                                        cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGGTCTGCGTGCTGGTCAGCGCGCCGATCt                                                            <  2:2961614/101‑1 (MQ=255)
                                        cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCt                                                            <  1:1578034/101‑1 (MQ=255)
                                         aTCGTGGCCTGGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                           <  2:2215312/101‑1 (MQ=255)
                                         aTCGTGGCCGAGGCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                           >  1:2422813/1‑101 (MQ=255)
                                             tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGcc                                                       >  1:2796337/1‑101 (MQ=255)
                                              ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcaacg                                                      >  1:1995349/1‑97 (MQ=255)
                                              ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                      <  1:2016558/101‑1 (MQ=255)
                                               gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                     <  2:1157123/101‑1 (MQ=255)
                                               gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                     <  1:1806921/101‑1 (MQ=255)
                                                ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAg                                                    <  1:225089/101‑1 (MQ=255)
                                                  gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                  >  1:206418/1‑101 (MQ=255)
                                                  gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAGCTGCGCCGAGcc                                                  >  2:2024083/1‑101 (MQ=255)
                                                  gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCCTGCTGGTCCGCGCGCCGATCTGCGCCGAGcc                                                  >  1:3056405/1‑101 (MQ=255)
                                                   agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGccc                                                 >  2:4395246/1‑101 (MQ=255)
                                                     cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                               <  1:971956/101‑1 (MQ=255)
                                                      ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGa                                              >  2:43273/1‑101 (MQ=255)
                                                      ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGa                                              >  2:3122573/1‑101 (MQ=255)
                                                        agcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcg                                            <  1:4397344/101‑1 (MQ=255)
                                                          cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgcg                                          <  2:2830636/101‑1 (MQ=255)
                                                           ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACgcgc                                         <  1:2353246/101‑1 (MQ=255)
                                                              gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTAGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgccgc                                      >  2:3220290/1‑101 (MQ=255)
                                                               agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCgccgcc                                     <  1:1329364/101‑1 (MQ=255)
                                                                 agCGCCGTCGCCCTGATCTGCGCCAGATCCACCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTg                                   <  2:2125952/101‑1 (MQ=255)
                                                                  gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgc                                  <  1:477497/101‑1 (MQ=255)
                                                                    gcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgcgc                                >  2:227220/1‑101 (MQ=255)
                                                                     ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                               >  1:1376502/1‑101 (MQ=255)
                                                                     ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                               <  2:1107803/101‑1 (MQ=255)
                                                                         cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                           >  2:2523639/1‑101 (MQ=255)
                                                                         cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                           >  1:1057054/1‑101 (MQ=255)
                                                                         cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                           >  1:1595208/1‑101 (MQ=255)
                                                                          gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCACAGACGCGCCGCCTGCGCcgtcgt                          >  1:1750279/1‑101 (MQ=255)
                                                                           cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                         <  2:2726811/101‑1 (MQ=255)
                                                                            gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                        >  2:2484861/1‑101 (MQ=255)
                                                                            gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                        <  1:1386350/101‑1 (MQ=255)
                                                                             cggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAg                       <  2:4092945/91‑1 (MQ=2)
                                                                             cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAg                       >  2:3567474/1‑101 (MQ=255)
                                                                              ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                      <  2:4043686/92‑1 (MQ=2)
                                                                              ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                      <  1:738253/92‑1 (MQ=2)
                                                                               gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                     <  1:123826/93‑1 (MQ=2)
                                                                               gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                     <  2:1592623/93‑1 (MQ=2)
                                                                               gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                     <  2:2636672/101‑1 (MQ=255)
                                                                               gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                     >  2:1361671/1‑101 (MQ=255)
                                                                               gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                     <  1:2598551/101‑1 (MQ=255)
                                                                                   tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                 <  2:14506/101‑1 (MQ=255)
                                                                                    gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACg                >  2:917669/1‑101 (MQ=255)
                                                                                    gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                <  1:3196885/98‑1 (MQ=2)
                                                                                    gcaccAGATCGGCGGCTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCCTCGACCCCAGCGCCGCCGCCTGCGCCGTCGTCACGGCCACg                >  1:3295658/4‑101 (MQ=255)
                                                                                     caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa               >  2:1485320/3‑101 (MQ=2)
                                                                                      cccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt              <  2:1298262/100‑1 (MQ=2)
                                                                                      accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt              <  2:704030/100‑1 (MQ=9)
                                                                                       ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg             <  1:570036/101‑1 (MQ=255)
                                                                                       ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg             >  1:2843771/1‑101 (MQ=255)
                                                                                        cAGATCCGCCGCCTAGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt            <  1:205131/101‑1 (MQ=255)
                                                                                            tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCACCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc        <  2:1056966/101‑1 (MQ=255)
                                                                                             ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg       <  1:1628757/101‑1 (MQ=255)
                                                                                              ccccgccCCGGGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg      >  2:3673217/3‑101 (MQ=255)
                                                                                                  gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc  <  1:1151264/101‑1 (MQ=255)
                                                                                                    |                                                                                                  
GCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCC  >  NC_014034/681347‑681545

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG