Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,568 C→T 7.9% 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 7.9% 34.2 0.0 ‑0.1189E347K (GAA→AAA) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (6/9):  ref base (72/102):  total (78/111)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.49e-01

Only 100 of 203 total aligned reads displayed.

GCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCG  >  NC_014034/682468‑682668
                                                                                                    |                                                                                                     
gcggcCTCCTCGGACAGCTCCGCGATCTGCNGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                                       <  1:1386888/101‑1 (MQ=255)
gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                                                       >  1:3358585/1‑101 (MQ=255)
     ctcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCg                                                                                                  <  1:556014/101‑1 (MQ=255)
     ctcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCg                                                                                                  >  1:3291258/1‑101 (MQ=255)
       cctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.Tg                                                                                               <  1:1263303/101‑1 (MQ=255)
        ctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGc                                                                                              >  2:4648367/1‑101 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCGGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             <  2:4032081/101‑1 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             >  1:1344043/1‑101 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             >  1:4642406/1‑101 (MQ=255)
         tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                                             >  2:3166475/1‑101 (MQ=255)
          cGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTc                                                                                            >  2:1110278/1‑101 (MQ=255)
           ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                                           <  2:3817032/101‑1 (MQ=255)
           gcaCAGCGCCGCGATCTACCGCGTGCCCAGAACCCCCAACTGCCCCGTCCTCACGGCCGCCATCTGCGCCCCCGTCTCCCCCGCCATTCCCATCG.TGCTcg                                                                                           >  2:2855768/3‑100 (MQ=255)
            gACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAg                                                                                          <  2:1392703/101‑1 (MQ=255)
               aGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGccc                                                                                       <  1:2995205/101‑1 (MQ=255)
                  gccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGc                                                                                    >  2:536293/1‑101 (MQ=255)
                   ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                                                   <  1:4172715/101‑1 (MQ=255)
                    cgcgCTCTTGCGCGTGCCCAGACCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCa                                                                                  <  2:896729/101‑1 (MQ=38)
                     gcgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCAt                                                                                 >  1:2054857/1‑101 (MQ=255)
                        aTCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTc                                                                              <  1:2473249/101‑1 (MQ=255)
                            gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                                          <  2:1879835/101‑1 (MQ=255)
                            gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                                          <  2:936180/101‑1 (MQ=255)
                              cgcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                        <  2:391157/101‑1 (MQ=255)
                              cgcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                        <  1:1128282/101‑1 (MQ=255)
                               gcgTGCCCAGAACCCCCAGCTGGTCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                                                       <  2:1514111/101‑1 (MQ=255)
                               gcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                                                       <  1:2786758/101‑1 (MQ=255)
                                 gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                     >  2:2163579/1‑101 (MQ=255)
                                 gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                     >  2:3336365/1‑101 (MQ=255)
                                  tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                                    >  1:1021938/1‑101 (MQ=255)
                                   gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                   >  2:1345404/1‑101 (MQ=255)
                                    cccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAg                                                                  <  1:2615838/101‑1 (MQ=255)
                                     ccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGc                                                                 >  2:474960/1‑101 (MQ=255)
                                         aaCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCt                                                             <  2:3633457/101‑1 (MQ=255)
                                          aCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCtt                                                            <  1:2639510/101‑1 (MQ=255)
                                           ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                                           <  1:3083497/98‑2 (MQ=37)
                                           cccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                                           <  2:4531013/101‑1 (MQ=255)
                                           cccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                                           <  1:3309976/101‑1 (MQ=255)
                                           cccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                                           <  1:1834894/101‑1 (MQ=255)
                                            cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                                                          >  1:3597997/3‑99 (MQ=37)
                                            cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                                                          >  2:1468752/3‑99 (MQ=37)
                                            ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                          <  2:1320517/101‑1 (MQ=255)
                                            ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TG..CAGCTCGGTCATGTCATAGATCGGAAgagcacacgtc                                                        >  1:432844/1‑91 (MQ=25)
                                             cccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAg                                                         <  1:4723506/101‑1 (MQ=255)
                                              ccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgc                                                        <  1:1611283/101‑1 (MQ=25)
                                                aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                                      >  1:2389598/1‑101 (MQ=255)
                                                aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                                      >  2:1011348/1‑101 (MQ=255)
                                                aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcgc                                                      <  2:2466930/101‑1 (MQ=25)
                                                 gCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGcgcg                                                     >  2:2495135/1‑101 (MQ=255)
                                                  cTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGcgcgc                                                    >  1:4754103/1‑100 (MQ=255)
                                                   tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                   >  2:4359047/1‑101 (MQ=255)
                                                   tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                   >  2:4998914/1‑101 (MQ=255)
                                                    gggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggcg                                                  <  2:738270/101‑1 (MQ=255)
                                                      gCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGt                                                >  2:1687509/1‑101 (MQ=255)
                                                      gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                <  2:2368605/101‑6 (MQ=21)
                                                       ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                               >  2:2823665/1‑101 (MQ=255)
                                                         gTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGcc                                             <  2:46197/101‑1 (MQ=255)
                                                          tGCTCAGCGCCGCCATCTTCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCg                                            >  2:3738516/1‑101 (MQ=255)
                                                          tGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCg                                            >  1:436853/1‑101 (MQ=255)
                                                            cTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGAt                                          >  2:4218573/1‑101 (MQ=255)
                                                             tCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                                         <  1:1774271/101‑1 (MQ=255)
                                                             tCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                                         <  1:2556462/101‑1 (MQ=255)
                                                              cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                        >  2:4244930/1‑101 (MQ=255)
                                                              cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                        >  2:3056828/1‑101 (MQ=255)
                                                               aGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgc                                       >  2:3188072/1‑101 (MQ=255)
                                                                gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCACCCCGGCCATGTCCTCGGTCGAAAGCCCGTTCAGCGCGGCGGTGCCGATcccc                                      >  2:297821/1‑101 (MQ=255)
                                                                 cgccgcCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccg                                     <  1:2660942/101‑1 (MQ=255)
                                                                   ccggcATCCGCGCCCCGGCCCGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   <  1:3564818/101‑1 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   <  1:3424451/101‑1 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   <  1:1590026/101‑1 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   <  1:716736/101‑1 (MQ=255)
                                                                   ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                   <  1:2163734/101‑1 (MQ=255)
                                                                     gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                 >  1:3684156/1‑101 (MQ=255)
                                                                          cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                            >  1:3456272/1‑101 (MQ=255)
                                                                          cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                            >  1:4852117/1‑101 (MQ=255)
                                                                          cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                            <  1:3166726/101‑1 (MQ=255)
                                                                                cccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCgg                      <  2:2645882/101‑1 (MQ=255)
                                                                                  cGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaa                    <  1:629904/101‑1 (MQ=255)
                                                                                   ggTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaa                   >  2:4468637/1‑101 (MQ=255)
                                                                                    gTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaaa                  <  1:897490/101‑1 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  2:3138130/1‑101 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                <  1:453346/101‑1 (MQ=255)
                                                                                      cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                >  2:1705666/1‑101 (MQ=255)
                                                                                        gcgcCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAATGCCGCGCCGGAAAAGgcg              >  2:4869078/1‑101 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  2:1133609/101‑1 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  1:305085/101‑1 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  1:807958/101‑1 (MQ=255)
                                                                                           ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg           <  1:4420005/101‑1 (MQ=255)
                                                                                             gccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGAt         <  1:3738794/101‑1 (MQ=255)
                                                                                              ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc        <  2:3049152/101‑1 (MQ=255)
                                                                                              ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCCCCCGGGAAAAGGCCCCgttt        >  2:3761368/1‑98 (MQ=255)
                                                                                               cAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt       >  2:4155948/1‑101 (MQ=255)
                                                                                                aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCGGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg      <  1:3991324/101‑1 (MQ=255)
                                                                                                aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg      >  2:3861406/1‑101 (MQ=255)
                                                                                                 gTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGc     >  1:4949378/1‑101 (MQ=255)
                                                                                                  ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt    <  1:2924968/101‑1 (MQ=255)
                                                                                                  ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt    <  2:579389/101‑1 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   <  1:3861691/101‑1 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   >  1:508833/1‑101 (MQ=255)
                                                                                                   tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc   >  1:1587873/1‑101 (MQ=255)
                                                                                                    cGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCg  >  1:1963586/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCG  >  NC_014034/682468‑682668

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