Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,120 G→C 47.7% V829V (GTC→GTG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,1200G→C47.7% 127.2 ‑0.7 ‑0.986V829V (GTC→GTG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (25/16):  ref base (24/21):  total (49/37)
Fisher's exact test strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.33e-01

GCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681020‑681217
                                                                                                    |                                                                                                 
gCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGtgcg                                                                                                   <  1:3440930/101‑1 (MQ=255)
gCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGcgcg                                                                                                   <  1:3612383/101‑1 (MQ=255)
 ccGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcga                                                                                                  >  2:1233737/1‑101 (MQ=255)
  cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 <  1:3581714/101‑1 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                >  2:2698765/1‑101 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                >  1:98785/1‑101 (MQ=255)
    ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               >  2:424094/1‑101 (MQ=255)
    ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               >  1:124984/1‑101 (MQ=255)
    ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               <  1:1778732/101‑1 (MQ=255)
     atctgATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  2:1406085/97‑1 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             <  1:1475955/101‑1 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             >  2:423705/1‑101 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             <  2:3376670/101‑1 (MQ=255)
       cAGCTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGt                                                                                            <  2:503742/101‑1 (MQ=255)
        aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           <  2:3841679/101‑1 (MQ=255)
        aGATCCTCGACCTCCAACGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGCTCGCCGGTCAGCGCGACGATGTc                                                                                           >  2:2465269/1‑101 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAtagat                                                                                          >  1:1406192/1‑97 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  1:3668575/1‑100 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  1:2221485/1‑100 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  1:2564391/1‑100 (MQ=255)
          aTCCTCGCCCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         <  2:1715254/101‑3 (MQ=255)
          aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         >  2:348246/1‑99 (MQ=255)
          aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         >  2:981732/1‑99 (MQ=255)
          aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         >  1:2643372/1‑99 (MQ=255)
           tCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgcc                                                                                        <  1:1090896/101‑4 (MQ=255)
            cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGCCGATGTcgccg                                                                                       >  2:2150415/1‑97 (MQ=255)
            cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccg                                                                                       <  2:1948939/101‑5 (MQ=255)
              tcggcctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  1:3999496/101‑1 (MQ=255)
              tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  2:2688422/101‑1 (MQ=255)
              tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  2:3898223/101‑1 (MQ=255)
                gccctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:3110741/99‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:918091/101‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:3507612/101‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:1712535/101‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCAGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:3745959/101‑3 (MQ=255)
                 acctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  >  2:3875821/1‑98 (MQ=255)
                  cctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgac                                                                                 <  2:3932646/101‑5 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                <  1:1987424/101‑6 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                <  2:3525930/101‑6 (MQ=255)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  2:2687679/1‑94 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  1:4100770/1‑94 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGACATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              <  1:2535682/101‑8 (MQ=38)
                      gaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcc                                                                             >  2:2415352/1‑93 (MQ=37)
                      gaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcc                                                                             <  2:880062/101‑9 (MQ=38)
                      gaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcc                                                                             >  1:2581855/1‑93 (MQ=37)
                       aTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagccc                                                                            >  2:1339444/1‑92 (MQ=37)
                                  cTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgc                                                                 >  1:2107865/1‑99 (MQ=25)
                                   tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                >  2:199625/1‑98 (MQ=25)
                                   tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                >  2:721911/1‑98 (MQ=25)
                                    gcgcCGTGCTAAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                               >  2:1045940/1‑97 (MQ=255)
                                     cgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccac                                                              >  1:2687819/1‑96 (MQ=25)
                                       ccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactg                                                            >  1:2297696/1‑94 (MQ=25)
                                         gTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggc                                                          >  2:4264993/1‑92 (MQ=21)
                                          tGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggcc                                                         >  1:3112244/1‑91 (MQ=255)
                                          tGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggcc                                                         >  2:1258836/1‑91 (MQ=25)
                                                                                  ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGGGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  1:3731015/11‑101 (MQ=37)
                                                                                  ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 <  1:823420/91‑1 (MQ=37)
                                                                                  ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  2:3459068/11‑101 (MQ=37)
                                                                                   tgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgc                >  2:3884138/10‑100 (MQ=37)
                                                                                   gtcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgt                >  2:905471/10‑101 (MQ=38)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               >  2:3080092/9‑101 (MQ=38)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               >  1:1656611/9‑101 (MQ=38)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  1:2216150/93‑1 (MQ=255)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  1:2477200/93‑1 (MQ=255)
                                                                                     cgccgtcgTCAGCGCCACGTTGACATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg              >  1:3499414/8‑101 (MQ=38)
                                                                                      gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt             <  2:4178248/95‑1 (MQ=255)
                                                                                      gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt             >  2:3136335/7‑101 (MQ=255)
                                                                                       ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg            >  1:2299009/6‑101 (MQ=255)
                                                                                        cgtcgTCCGCGCCGCGATGTCATCGGTGGCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCCAGCTGCGTCGTGc           >  1:2979692/5‑101 (MQ=255)
                                                                                        cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc           >  2:2391040/5‑101 (MQ=255)
                                                                                        cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc           >  1:906473/5‑101 (MQ=255)
                                                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt          <  1:4306569/98‑1 (MQ=255)
                                                                                          tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  2:2815130/99‑1 (MQ=255)
                                                                                          tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  2:2916626/99‑1 (MQ=255)
                                                                                           cgtgaGCGCCACGATGTCATCGGTGCCGAGGCCATTGATCGCATTGGCTCTAAGCGCTGTCAAGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  1:1041198/97‑1 (MQ=37)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:987706/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  1:1068242/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  1:2604043/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  2:466137/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  1:3563682/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  1:80315/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:1349224/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:3021743/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  2:2691636/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  2:2363738/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  2:3252325/2‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:2851909/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:3257642/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:3591720/1‑101 (MQ=255)
                                                                                             tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc      <  1:1659381/101‑1 (MQ=255)
                                                                                             tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc      <  2:1827950/101‑1 (MQ=255)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     >  2:2535232/1‑101 (MQ=255)
                                                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc    >  2:2456589/1‑101 (MQ=255)
                                                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc    >  2:2523015/1‑101 (MQ=255)
                                                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc    <  1:1485425/101‑1 (MQ=255)
                                                                                                gcgcCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGcgcg   <  2:4431867/101‑1 (MQ=255)
                                                                                                gcgcCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGcgcg   <  2:3664115/101‑1 (MQ=255)
                                                                                                 cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCgg  <  1:668269/101‑1 (MQ=255)
                                                                                                 cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCACgg  >  1:3997554/1‑101 (MQ=255)
                                                                                                    |                                                                                                 
GCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681020‑681217

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