Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,450 G→C 9.1% L719L (CTC→CTG RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4500G→C 9.1% 30.0 ‑4.1 ‑0.0176L719L (CTC→CTGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/13):  ref base (83/77):  total (86/90)
Fisher's exact test strand distribution p-value = 1.63e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

Only 100 of 192 total aligned reads displayed.

TCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGC  >  NC_014034/681350‑681550
                                                                                                    |                                                                                                     
tCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                       <  2:2321493/101‑1 (MQ=255)
tCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                       >  1:1452562/1‑101 (MQ=255)
   ctaGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgc                                                                                                    <  2:3115007/99‑1 (MQ=255)
      gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                                 <  2:986524/101‑1 (MQ=255)
      gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                                 <  2:4423443/101‑1 (MQ=255)
        gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                               <  2:1171606/101‑1 (MQ=255)
             gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                          >  2:3246831/1‑101 (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)
                   cgcgcCTGAAAGCCCAGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGc                                                                                    >  1:2396552/1‑101 (MQ=255)
                    tcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                                   <  2:5041801/100‑1 (MQ=255)
                    gcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                                   <  1:2843335/101‑1 (MQ=255)
                        cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                               <  1:4390299/101‑1 (MQ=255)
                        cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                               <  2:2032147/101‑1 (MQ=255)
                        cTGAAAGCCCGACCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                               >  1:1955285/1‑101 (MQ=255)
                          gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTTCGTGCTGGTCAg                                                                             >  1:2027285/1‑101 (MQ=255)
                          gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                             <  2:1845928/101‑1 (MQ=255)
                            aaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcg                                                                           >  1:4505818/1‑101 (MQ=255)
                             aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                          <  2:1368828/101‑1 (MQ=255)
                             aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                          <  2:3048653/101‑1 (MQ=255)
                             aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCGCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                          <  1:4481860/101‑1 (MQ=255)
                                   gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                                    <  1:690547/101‑1 (MQ=255)
                                    ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                                   <  2:4616803/101‑1 (MQ=255)
                                    ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                                   >  1:3072881/1‑101 (MQ=255)
                                     cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCt                                                                  >  1:3875933/1‑101 (MQ=255)
                                      aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGCGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                                 >  1:4751252/1‑101 (MQ=255)
                                       tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgc                                                                <  2:2505273/101‑1 (MQ=255)
                                       tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgc                                                                >  2:2026532/1‑101 (MQ=255)
                                        cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcg                                                               >  2:3969248/1‑101 (MQ=255)
                                         gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcgc                                                              <  1:2085021/101‑1 (MQ=255)
                                          tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGcc                                                             <  1:4487938/101‑1 (MQ=255)
                                           ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGGGCCGATCTGCGCCg                                                            >  2:2361997/1‑101 (MQ=255)
                                           ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                            >  2:1797641/1‑101 (MQ=255)
                                            gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                           >  1:3513976/1‑101 (MQ=255)
                                             ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAg                                                          <  2:3735592/101‑1 (MQ=255)
                                              cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                         <  2:53838/101‑1 (MQ=255)
                                               gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                        >  2:785457/1‑101 (MQ=255)
                                               gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                        <  2:2828017/101‑1 (MQ=255)
                                               gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                        <  1:547224/101‑1 (MQ=255)
                                                  cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAg                                                     >  1:3356525/1‑101 (MQ=255)
                                                      gcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgc                                                 <  1:3949258/101‑1 (MQ=255)
                                                       cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgc.g                                               >  1:4137795/1‑101 (MQ=255)
                                                         cGGAGAGCTCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                             <  1:4246863/101‑1 (MQ=255)
                                                         cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                             <  2:4839140/101‑1 (MQ=255)
                                                         cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                             <  1:4912142/101‑1 (MQ=255)
                                                         cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gcc                                             <  1:371856/101‑1 (MQ=255)
                                                          ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGc.gccg                                            >  2:330817/1‑101 (MQ=255)
                                                            agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.gccgcc                                          <  1:1019758/101‑1 (MQ=255)
                                                            agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.gccgcc                                          >  2:4523938/1‑101 (MQ=255)
                                                                  ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGcc                                    <  1:4075698/101‑1 (MQ=255)
                                                                   cgccgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATGTGCGTGCTGGTCAGGGCGCCGATCTGCGCCGAGCCAAGACGC.GCCGCCTGCGCcg                                   <  2:619325/98‑1 (MQ=25)
                                                                   cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcg                                   >  2:4075235/1‑101 (MQ=255)
                                                                    gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgt                                  >  2:48974/1‑101 (MQ=255)
                                                                     tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtc                                 <  1:2542191/101‑1 (MQ=255)
                                                                     gcgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtc                                 <  2:1452435/100‑1 (MQ=255)
                                                                      cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcg                                >  2:4795266/1‑101 (MQ=255)
                                                                      cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcg                                >  2:2406248/1‑101 (MQ=255)
                                                                        cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCgtcgtc                              >  1:4803815/1‑101 (MQ=255)
                                                                          cggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAg                            <  1:3216274/91‑1 (MQ=2)
                                                                          cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAg                            >  2:3578378/1‑101 (MQ=255)
                                                                           tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GGCGCCTGCGCCGTCGTCAgc                           >  1:2964804/1‑101 (MQ=255)
                                                                           ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAgc                           <  2:2056722/92‑1 (MQ=2)
                                                                               ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCa                       <  1:5008976/96‑1 (MQ=2)
                                                                               cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCa                       <  1:3995673/101‑1 (MQ=255)
                                                                                cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCAc                      <  2:4606429/97‑1 (MQ=2)
                                                                                 gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACg                     <  2:4261747/101‑1 (MQ=255)
                                                                                 gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACg                     >  2:3791042/1‑101 (MQ=255)
                                                                                 gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCTTGCTCGACAGGGCCAGAAGCTGCGTCGAGCCCAG.CGCGGCCGCCTGCGCCGTCGTCAGCCCCACg                     >  2:1636690/4‑101 (MQ=2)
                                                                                 gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg                     <  2:2598648/98‑1 (MQ=2)
                                                                                 gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg                     <  1:3925425/98‑1 (MQ=2)
                                                                                  caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGa                    <  2:1836725/99‑1 (MQ=9)
                                                                                   gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGAt                   <  1:18655/101‑1 (MQ=255)
                                                                                   accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGAt                   <  1:3689095/100‑1 (MQ=2)
                                                                                    ccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                  <  2:2788658/101‑1 (MQ=9)
                                                                                    ccAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                  >  1:2738498/1‑101 (MQ=255)
                                                                                    ccAGATCCGCCGCCTGGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                  <  1:889882/101‑1 (MQ=255)
                                                                                       gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa               >  2:4499916/1‑101 (MQ=255)
                                                                                        aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt              <  2:3819267/101‑1 (MQ=255)
                                                                                        aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt              >  1:3870981/1‑101 (MQ=255)
                                                                                        aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAAGC.GCAGCCTGCGCCGTCGTCAGCGCCACGATGTCAt              >  2:2497699/1‑101 (MQ=255)
                                                                                          ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg            >  1:4419831/1‑101 (MQ=255)
                                                                                          ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg            <  1:4181143/101‑1 (MQ=255)
                                                                                          ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg            <  1:3280713/101‑1 (MQ=255)
                                                                                            gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt          <  1:1196055/101‑1 (MQ=255)
                                                                                              cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGNCACGATGTCATCGGTGc        >  2:1535291/1‑101 (MQ=255)
                                                                                               gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc       <  2:4154348/101‑1 (MQ=255)
                                                                                                 cTCGAGCGCCGCGATCTGGGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa     <  1:1364855/101‑1 (MQ=255)
                                                                                                 cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa     >  2:1854572/1‑101 (MQ=255)
                                                                                                 cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTAGTCAGCGCCACGATGTCATCGGTGCCaa     <  1:1034975/101‑1 (MQ=255)
                                                                                                  tCGAGCGCCGCGCTCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa    >  2:3519808/1‑101 (MQ=255)
                                                                                                  tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa    <  1:1532426/101‑1 (MQ=255)
                                                                                                   cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg   <  2:4749914/101‑1 (MQ=255)
                                                                                                   cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg   >  2:5019028/1‑101 (MQ=255)
                                                                                                   cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg   >  1:852211/1‑101 (MQ=255)
                                                                                                    gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc  <  1:2579101/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
TCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGC  >  NC_014034/681350‑681550

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