Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,444 C→G 9.6% A721A (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,4440C→G 9.6% 31.0 ‑4.1 ‑0.1163A721A (GCG→GCC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/13):  ref base (77/70):  total (80/83)
Fisher's exact test strand distribution p-value = 1.58e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.75e-01

Only 100 of 187 total aligned reads displayed.

GCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGC  >  NC_014034/681344‑681544
                                                                                                    |                                                                                                    
gCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgcc                                                                                                      >  2:4710085/1‑101 (MQ=255)
 ccGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgccg                                                                                                     <  1:4937373/101‑1 (MQ=255)
 ccGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgccg                                                                                                     <  1:613770/101‑1 (MQ=255)
   gcgTCGACCGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGTCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                   <  2:2451869/101‑1 (MQ=255)
     gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                 <  1:1954177/101‑1 (MQ=255)
     gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                 >  2:2845/1‑101 (MQ=255)
     gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                 >  2:2888385/1‑101 (MQ=255)
     gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                 <  1:2888598/101‑1 (MQ=255)
      tCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                <  2:2321493/101‑1 (MQ=255)
      acggcAGCGCCGCAATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                <  1:1519002/97‑1 (MQ=255)
       cGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGa                                                                                               <  2:355524/101‑1 (MQ=255)
        gACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAg                                                                                              >  2:2352019/1‑101 (MQ=255)
         ctaGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgc                                                                                             <  2:3115007/99‑1 (MQ=255)
           aGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgcag                                                                                           >  1:3115245/1‑99 (MQ=255)
            gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                          <  1:4393216/101‑1 (MQ=255)
            gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                          <  2:986524/101‑1 (MQ=255)
            gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                          <  2:4423443/101‑1 (MQ=255)
              gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                        <  2:1171606/101‑1 (MQ=255)
                cgcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGa                                                                                      >  2:4309452/1‑101 (MQ=255)
                   gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                   >  2:3246831/1‑101 (MQ=255)
                   gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                   <  2:2097249/101‑1 (MQ=255)
                    aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                  >  2:1129468/1‑101 (MQ=255)
                     tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                 <  1:4519515/101‑1 (MQ=255)
                      cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                >  2:3995371/1‑101 (MQ=255)
                      cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                >  1:514676/1‑101 (MQ=255)
                       tGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGt                                                                               <  1:1656132/101‑1 (MQ=255)
                        gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                              >  2:2928914/1‑101 (MQ=255)
                        gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                              >  2:889777/1‑101 (MQ=255)
                        gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                              <  2:2013700/101‑1 (MQ=255)
                         cgcgcCTGAAGGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCCGATCCGCCGCCTCGAGCGGCGCGAGCTGCGTGc                                                                             >  1:4363226/1‑101 (MQ=255)
                          tcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                            <  2:5041801/100‑1 (MQ=255)
                              cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                        <  2:2056033/101‑1 (MQ=255)
                                gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTTCGTGCTGGTCAg                                                                      >  1:2027285/1‑101 (MQ=255)
                                   aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGGCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                   >  2:272043/1‑101 (MQ=255)
                                   aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                   <  2:1368828/101‑1 (MQ=255)
                                   aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                   <  2:3048653/101‑1 (MQ=255)
                                      ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGNGTGCTGGTCAGCGCGcc                                                                >  2:1838227/1‑101 (MQ=255)
                                       cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                               >  2:1532296/1‑101 (MQ=255)
                                         gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                             <  2:4790044/101‑1 (MQ=255)
                                          ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                            >  1:3072881/1‑101 (MQ=255)
                                          ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                            <  1:1387547/101‑1 (MQ=255)
                                           cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCt                                                           >  1:3875933/1‑101 (MQ=255)
                                            aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGCGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                          >  1:4751252/1‑101 (MQ=255)
                                            aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                          >  2:5014923/1‑101 (MQ=255)
                                              cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcg                                                        >  2:3969248/1‑101 (MQ=255)
                                              cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGACCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcg                                                        <  2:2595383/101‑1 (MQ=255)
                                               gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcgc                                                       <  1:2085021/101‑1 (MQ=255)
                                                tggcttagcgcagcccGGAGAGCGCCGTCGCGCGGATCTACGTCAGATCTTCCGAATAGAGCGCCGCGTTCTGCGTGGTGGTCAGCGCGCCGATCTGCGcc                                                      <  2:2027141/95‑1 (MQ=19)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                     >  2:1797641/1‑101 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                     >  1:4486761/1‑101 (MQ=255)
                                                  gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                    >  1:3513976/1‑101 (MQ=255)
                                                    cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                  <  2:53838/101‑1 (MQ=255)
                                                     gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                 >  2:785457/1‑101 (MQ=255)
                                                     gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                 <  1:547224/101‑1 (MQ=255)
                                                        cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                              >  1:3356525/1‑101 (MQ=255)
                                                        cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                              <  2:4143183/101‑1 (MQ=255)
                                                            gcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgc                                          <  1:3949258/101‑1 (MQ=255)
                                                             cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgcg                                         >  1:4137795/1‑101 (MQ=255)
                                                               cGGAGAGCTCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgcc                                       <  1:4246863/101‑1 (MQ=255)
                                                               cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgcc                                       <  2:4839140/101‑1 (MQ=255)
                                                                ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgccg                                      >  2:330817/1‑101 (MQ=255)
                                                                 gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgccgc                                     >  1:4873147/1‑101 (MQ=255)
                                                                  agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCgccgcc                                    >  2:3818583/1‑101 (MQ=255)
                                                                  agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCgccgcc                                    >  2:4523938/1‑101 (MQ=255)
                                                                  agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCgccgcc                                    <  1:1019758/101‑1 (MQ=255)
                                                                    agCGCCGTCGCGCTGATCTGCGCTAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTg                                  <  2:1829619/101‑1 (MQ=255)
                                                                      cgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgcg                                >  1:1941853/1‑101 (MQ=255)
                                                                       gcCGTCGCGCTGATCTGCGCCAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgcgc                               >  1:4490819/1‑101 (MQ=255)
                                                                        ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                              <  1:4075698/101‑1 (MQ=255)
                                                                         cgccgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATGTGCGTGCTGGTCAGGGCGCCGATCTGCGCCGAGCCAAGACGCGCCGCCTGCGCcg                             <  2:619325/98‑1 (MQ=25)
                                                                         cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                             >  2:4075235/1‑101 (MQ=255)
                                                                         cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                             <  2:1850689/101‑1 (MQ=255)
                                                                          gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgt                            >  2:48974/1‑101 (MQ=255)
                                                                            cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                          >  2:2406248/1‑101 (MQ=255)
                                                                              cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                        >  1:4803815/1‑101 (MQ=255)
                                                                               gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                       <  1:3279146/101‑1 (MQ=255)
                                                                                cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAg                      >  2:3578378/1‑101 (MQ=255)
                                                                                cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGACTGCGCCGTCGCCAg                      >  1:4445062/1‑101 (MQ=255)
                                                                                 tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGGCGCCTGCGCCGTCGTCAgc                     >  1:2964804/1‑101 (MQ=255)
                                                                                 ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                     <  1:1093843/92‑1 (MQ=2)
                                                                                     ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                 <  1:5008976/96‑1 (MQ=2)
                                                                                     cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                 <  1:3995673/101‑1 (MQ=255)
                                                                                      cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                <  1:901341/97‑1 (MQ=2)
                                                                                      cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                <  2:4606429/97‑1 (MQ=2)
                                                                                       gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg               <  1:3925425/98‑1 (MQ=2)
                                                                                        cgcCAGATCCGCCGCCTCGAGGGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGa              >  1:3717894/1‑101 (MQ=255)
                                                                                        cgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGa              >  1:794114/1‑101 (MQ=255)
                                                                                        cgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGa              >  2:1629356/1‑101 (MQ=255)
                                                                                         gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt             <  1:18655/101‑1 (MQ=255)
                                                                                         accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCCCGAt             >  2:2502085/2‑101 (MQ=2)
                                                                                         accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt             <  1:3689095/100‑1 (MQ=2)
                                                                                          ccAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg            >  1:2738498/1‑101 (MQ=255)
                                                                                          ccAGATCCGCCGCCTGGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg            <  1:889882/101‑1 (MQ=255)
                                                                                             gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa         >  2:4499916/1‑101 (MQ=255)
                                                                                              aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt        >  1:2413151/1‑101 (MQ=255)
                                                                                              aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt        >  1:3870981/1‑101 (MQ=255)
                                                                                               tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGATCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc       >  2:2825426/1‑101 (MQ=255)
                                                                                                ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg      <  1:4181143/101‑1 (MQ=255)
                                                                                                  gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCggg    >  2:3741002/1‑100 (MQ=255)
                                                                                                  gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt    <  1:1196055/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGC  >  NC_014034/681344‑681544

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