Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,568 C→T 9.6% 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.6% 46.6 ‑2.9 ‑0.1197E347K (GAA→AAA) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (5/14):  ref base (91/87):  total (96/101)
Fisher's exact test strand distribution p-value = 5.27e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.58e-01

Only 100 of 209 total aligned reads displayed.

GCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCG  >  NC_014034/682468‑682668
                                                                                                    |                                                                                                     
gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                                       <  2:2743757/101‑1 (MQ=255)
     ctcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGACAGCGCCGCCAGTTCGATCg                                                                                                  >  1:4723742/1‑101 (MQ=255)
        ctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGc                                                                                              >  1:3628364/1‑101 (MQ=255)
        ctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGc                                                                                              >  2:4171812/1‑101 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             >  1:1913697/1‑101 (MQ=255)
          cGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTc                                                                                            >  2:1987944/1‑101 (MQ=255)
          cGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTc                                                                                            <  1:2760164/101‑1 (MQ=255)
           ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                                           <  1:3240412/101‑1 (MQ=255)
           ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                                           <  2:1869111/101‑1 (MQ=255)
           ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                                           <  1:1619824/101‑1 (MQ=255)
           ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                                           >  2:74562/1‑101 (MQ=255)
            gACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAg                                                                                          <  2:3655904/101‑1 (MQ=255)
              cAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGcc                                                                                        <  1:2364968/101‑1 (MQ=255)
               aGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGccc                                                                                       <  1:1076635/101‑1 (MQ=255)
                gcgcCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCg                                                                                      <  1:2020256/101‑1 (MQ=255)
                  gccgcGATCTGCCGAGTGCCCAGAACCCACGGCTGGGCCGTCCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCCCCAGTTCGATCG.TGCTCAGCCCGGc                                                                                    <  2:3141095/101‑1 (MQ=255)
                   ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                                                   <  2:4520095/101‑1 (MQ=255)
                   ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGCCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                                                   >  2:2868992/1‑101 (MQ=255)
                     gcgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCAt                                                                                 <  1:2139128/101‑1 (MQ=255)
                        aTCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTc                                                                              >  2:2387493/1‑101 (MQ=255)
                           tGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTc                                                                           >  2:1927858/1‑101 (MQ=255)
                              cgcgTGCCCAGACCCCCCAGCGGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                        <  2:2974154/101‑1 (MQ=255)
                                cgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCg                                                                      <  1:2933333/101‑1 (MQ=255)
                                 gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                     >  2:2152582/1‑101 (MQ=255)
                                 gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                     >  1:1827081/1‑101 (MQ=255)
                                 gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                     >  1:650267/1‑101 (MQ=255)
                                  tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                                    >  2:3533789/1‑101 (MQ=255)
                                   gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                   >  1:3476918/1‑101 (MQ=255)
                                     ccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGc                                                                 <  1:1082765/101‑1 (MQ=255)
                                       cgAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                               <  2:905667/100‑1 (MQ=255)
                                       aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                               >  1:2467764/1‑101 (MQ=255)
                                         aaCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGACCCt                                                             >  2:1454725/1‑101 (MQ=255)
                                          gccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCtt                                                            <  1:1420330/97‑1 (MQ=37)
                                          aCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCtt                                                            <  2:4052756/101‑1 (MQ=255)
                                           ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                                           <  2:2600036/98‑2 (MQ=37)
                                              ccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGGTCCCGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgc                                                        <  1:4060728/101‑1 (MQ=255)
                                                aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                                      >  1:431124/1‑101 (MQ=255)
                                                 gCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGcgcg                                                     <  1:3335553/101‑1 (MQ=255)
                                                 gCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGcgcg                                                     >  1:1397443/1‑101 (MQ=255)
                                                 gCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcg                                                     <  2:4218152/101‑1 (MQ=25)
                                                 gCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcg                                                     <  2:426741/101‑1 (MQ=25)
                                                  cTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcgg                                                    <  2:3496815/101‑1 (MQ=255)
                                                   tGGCCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                   <  2:4061720/101‑1 (MQ=255)
                                                   cgggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                   <  2:1130500/100‑1 (MQ=255)
                                                     ggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGcggcgg                                                 <  2:3364864/101‑1 (MQ=255)
                                                     ggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGcggcgg                                                 >  1:4254797/1‑101 (MQ=255)
                                                     cgCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccac                                                 <  2:244022/100‑5 (MQ=25)
                                                      gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                <  2:2833706/101‑6 (MQ=21)
                                                      gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                <  1:4522819/101‑6 (MQ=21)
                                                       ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                               >  1:4587206/1‑101 (MQ=255)
                                                       ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                               <  1:36113/101‑1 (MQ=255)
                                                       ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                               <  1:906199/101‑1 (MQ=255)
                                                       ccGTGCTCAGCGCCGCCACCTGCTGCGAGGCCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccactg                                               <  2:3923360/101‑7 (MQ=21)
                                                        cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccactgg                                              <  1:4692772/101‑8 (MQ=25)
                                                          tGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCg                                            <  2:650203/101‑1 (MQ=255)
                                                          tGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCg                                            <  1:391503/101‑1 (MQ=255)
                                                            cTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGAt                                          >  2:2297673/1‑101 (MQ=255)
                                                              cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                        >  2:4446391/1‑101 (MQ=255)
                                                              cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                        <  2:3008287/101‑1 (MQ=255)
                                                              cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                        <  2:2467621/101‑1 (MQ=255)
                                                                gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc                                      <  2:1706223/101‑1 (MQ=255)
                                                                  gccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccgc                                    >  1:1959417/1‑101 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   <  2:1353024/101‑1 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   >  2:4742427/1‑101 (MQ=255)
                                                                    cgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaa                                  >  2:1091461/1‑101 (MQ=255)
                                                                    cgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGCCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGGGCCGATCGCCGCaa                                  >  2:516674/1‑101 (MQ=255)
                                                                     gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                 <  1:824250/101‑1 (MQ=255)
                                                                     gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                 <  1:2592249/101‑1 (MQ=255)
                                                                       cATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAgg                               <  1:2860437/101‑1 (MQ=255)
                                                                       cATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAgg                               >  2:2909949/1‑101 (MQ=255)
                                                                          cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                            <  1:1666574/101‑1 (MQ=255)
                                                                           tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc                           >  2:2019194/1‑101 (MQ=255)
                                                                           tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAAGCcgc                           >  2:3223827/1‑101 (MQ=255)
                                                                             cgcCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCgcgc                         >  2:3843941/1‑101 (MQ=255)
                                                                                   ggccAGCGCCGCCATTGCGGTCG.TGCTCAGCCCGGCCATGTCCTCGGGCGAAAGCCCCTTGAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaa                   <  1:1273404/98‑1 (MQ=255)
                                                                                    gTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaaa                  <  1:1361577/101‑1 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  1:2503225/1‑101 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  1:3665/1‑101 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  1:3114468/1‑101 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  1:3541182/1‑101 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  1:1039254/1‑101 (MQ=255)
                                                                                        gcgcCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGgcg              <  2:2097179/101‑1 (MQ=255)
                                                                                          gccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGcc            >  1:2794762/1‑101 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  1:2542685/101‑1 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           >  2:3892214/1‑101 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  1:4503962/101‑1 (MQ=255)
                                                                                              ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc        <  2:1384145/101‑1 (MQ=255)
                                                                                              ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc        <  2:4586953/101‑1 (MQ=255)
                                                                                               catttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt       <  2:2176992/98‑1 (MQ=255)
                                                                                               cAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt       <  1:74569/101‑1 (MQ=255)
                                                                                               cAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt       <  1:2047914/101‑1 (MQ=255)
                                                                                                aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg      >  1:3869257/1‑101 (MQ=255)
                                                                                                  ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt    <  2:4657871/101‑1 (MQ=255)
                                                                                                  ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt    >  1:2266128/1‑101 (MQ=255)
                                                                                                  ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt    >  2:262629/1‑101 (MQ=255)
                                                                                                  ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt    >  2:2785400/1‑101 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGTAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   >  1:4086610/1‑101 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   >  1:4408880/1‑101 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   <  2:870231/101‑1 (MQ=255)
                                                                                                    cGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCg  >  2:410076/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCG  >  NC_014034/682468‑682668

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG