Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,568 C→T 9.5% E347K (GAA→AAA)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034682,5680C→T 9.5% 59.2 ‑3.5 ‑0.0252E347K (GAA→AAA) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (5/19):  ref base (103/125):  total (108/144)
Fisher's exact test strand distribution p-value = 2.88e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.27e-01

Only 100 of 260 total aligned reads displayed.

GCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCG  >  NC_014034/682468‑682668
                                                                                                    |                                                                                                     
gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                                       <  1:2619489/101‑1 (MQ=255)
gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                                       >  1:1407869/1‑101 (MQ=255)
   gcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGAt                                                                                                    <  1:2956879/101‑1 (MQ=255)
    cctcctCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCCGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATc                                                                                                   <  1:2745970/101‑1 (MQ=255)
      tcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.t                                                                                                >  2:2859175/1‑101 (MQ=255)
       cctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.Tg                                                                                               >  1:3417981/1‑101 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             >  2:2675656/1‑101 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             >  2:1399925/1‑101 (MQ=255)
          cGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTc                                                                                            >  1:2886689/1‑101 (MQ=255)
             aCAGCGCCGCGATCTCCCGCGTGCCCAGACCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGc                                                                                         <  1:4563548/101‑1 (MQ=255)
              cggcgcCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGcc                                                                                        <  1:1873691/99‑1 (MQ=255)
                 cgccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCgg                                                                                     >  1:944043/1‑101 (MQ=255)
                   ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                                                   <  2:1245628/101‑1 (MQ=255)
                   ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                                                   <  1:361754/101‑1 (MQ=255)
                      cgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATg                                                                                <  2:4447415/101‑1 (MQ=255)
                       gATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGt                                                                               >  2:4323289/1‑101 (MQ=255)
                       gATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGt                                                                               >  2:3214659/1‑101 (MQ=255)
                       gATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGt                                                                               >  1:2032479/1‑101 (MQ=255)
                           tGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTc                                                                           <  2:71554/101‑1 (MQ=255)
                            gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                                          <  2:2141195/101‑1 (MQ=255)
                            gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                                          <  1:88971/101‑1 (MQ=255)
                            gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                                          >  1:1787036/1‑101 (MQ=255)
                             ccGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCgg                                                                         >  1:1222253/1‑101 (MQ=255)
                              cgcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                        <  1:1086165/101‑1 (MQ=255)
                               gcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                                                       <  1:1445328/101‑1 (MQ=255)
                                cgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCg                                                                      <  1:748541/101‑1 (MQ=255)
                                 gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                     >  2:419898/1‑101 (MQ=255)
                                   gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                   <  1:2943586/101‑1 (MQ=255)
                                   gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                   <  2:4225654/101‑1 (MQ=255)
                                     ccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGc                                                                 >  2:1359922/1‑101 (MQ=255)
                                     ccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGc                                                                 <  1:493666/101‑1 (MQ=255)
                                       aGAGCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                               <  1:1493073/101‑1 (MQ=255)
                                       aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                               >  2:2339384/1‑101 (MQ=255)
                                       aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                               >  2:652004/1‑101 (MQ=255)
                                        gAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGcccc                                                              <  2:3737773/101‑1 (MQ=255)
                                         cgccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCt                                                             <  1:1003370/96‑1 (MQ=37)
                                         cgccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCt                                                             >  2:2441010/6‑101 (MQ=37)
                                          gccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCtt                                                            <  2:1343206/97‑1 (MQ=37)
                                          gccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCtt                                                            >  1:135916/5‑101 (MQ=37)
                                           ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                                           <  1:4103193/98‑2 (MQ=37)
                                           ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                                           <  1:1718649/98‑2 (MQ=37)
                                           ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                                           <  2:147721/98‑2 (MQ=37)
                                           cccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                                           <  2:2050342/101‑1 (MQ=255)
                                           cccccAGCTCGGCCTCCCTCGTCCCCCCCACTTTCGCCCCGGTCAGCTCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                                           <  2:2314312/101‑1 (MQ=255)
                                            cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                                                          <  2:166127/99‑3 (MQ=37)
                                            ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                          <  2:1114567/101‑1 (MQ=255)
                                             cccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAg                                                         <  1:2711059/101‑1 (MQ=255)
                                             cccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAg                                                         >  2:3913621/1‑101 (MQ=255)
                                              ccAGCTGGGCCGTCCTCAGCGCCGCCATCTGCGTCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgc                                                        <  1:4009140/101‑1 (MQ=255)
                                              ccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgc                                                        <  1:2909646/101‑1 (MQ=25)
                                               cAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcg                                                       >  2:3834459/1‑101 (MQ=255)
                                               cAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcg                                                       >  2:2270660/1‑101 (MQ=25)
                                                aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                                      <  2:85535/101‑1 (MQ=255)
                                                aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                                      >  2:309080/1‑101 (MQ=255)
                                                  cTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcgg                                                    <  1:2275677/101‑1 (MQ=255)
                                                   tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                   >  2:1776078/1‑101 (MQ=255)
                                                   tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                   >  1:533009/1‑101 (MQ=255)
                                                     ggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGcggcgg                                                 >  1:3947150/1‑101 (MQ=255)
                                                     ggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGcggcgg                                                 >  1:657399/1‑101 (MQ=255)
                                                      gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                <  2:4524971/101‑6 (MQ=21)
                                                      gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                <  2:2175809/101‑6 (MQ=21)
                                                        cGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGc                                              <  1:2859391/101‑1 (MQ=255)
                                                        cGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGc                                              >  2:3633730/1‑101 (MQ=255)
                                                        cGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGc                                              <  2:775294/101‑1 (MQ=255)
                                                        cGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGc                                              <  2:2551579/101‑1 (MQ=255)
                                                         gTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGcc                                             <  1:2339546/101‑1 (MQ=255)
                                                           gCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGa                                           <  2:3720110/101‑1 (MQ=255)
                                                            cTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGAt                                          <  2:2889449/101‑1 (MQ=255)
                                                               aGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgc                                       <  1:3485614/101‑1 (MQ=255)
                                                               aGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgc                                       >  1:2561047/1‑101 (MQ=255)
                                                                gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc                                      >  1:3152702/1‑101 (MQ=255)
                                                                  gccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccgc                                    <  2:306475/101‑1 (MQ=255)
                                                                  gccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccgc                                    >  2:4627122/1‑101 (MQ=255)
                                                                  gccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccgc                                    <  2:1808149/101‑1 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   >  1:3869317/1‑101 (MQ=255)
                                                                     gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                 >  2:690988/1‑101 (MQ=255)
                                                                     gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                 <  1:1604196/101‑1 (MQ=255)
                                                                      ccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcagatcgg                                >  2:3954129/1‑94 (MQ=255)
                                                                       cATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAgg                               >  2:2421056/1‑101 (MQ=255)
                                                                        aTCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGc                              <  2:943963/101‑1 (MQ=255)
                                                                         tCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGcc                             <  2:2949676/101‑1 (MQ=255)
                                                                          cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAACGCcg                            >  2:1960084/1‑101 (MQ=255)
                                                                           tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc                           <  1:3291456/101‑1 (MQ=255)
                                                                           tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc                           >  1:4066069/1‑101 (MQ=255)
                                                                           tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc                           >  2:2648351/1‑101 (MQ=255)
                                                                            gcgcCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgcg                          <  1:419935/101‑1 (MQ=255)
                                                                                cccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCgg                      >  1:4277981/1‑101 (MQ=255)
                                                                                 ccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGa                     >  2:3464493/1‑101 (MQ=255)
                                                                                   ggTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaa                   <  1:1288526/101‑1 (MQ=255)
                                                                                   ggTCAGCGCCGCCAGTTCGATCG.GGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaa                   >  1:3751467/1‑101 (MQ=255)
                                                                                    gTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGGTCGCCGCAAAGGCCGCGCCGTaaaa                  >  1:4455142/1‑101 (MQ=255)
                                                                                     tCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAg                 <  2:1524011/101‑1 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  2:594194/101‑1 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCCGAAAAGGCGCCg           <  1:2530244/101‑1 (MQ=255)
                                                                                              ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc        >  2:1198627/1‑101 (MQ=255)
                                                                                              ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc        >  2:2899883/1‑101 (MQ=255)
                                                                                                aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg      >  2:2820512/1‑101 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   <  1:691050/101‑1 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   >  1:2116593/1‑101 (MQ=255)
                                                                                                    cGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCg  <  1:1776197/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
GCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCG  >  NC_014034/682468‑682668

Base quality scores:  ATCG  < 3 ≤  ATCG  < 22 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 40 ≤  ATCG