Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,120 G→C 31.4% 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→C31.4% 57.6 ‑0.2 ‑0.070V829V (GTC→GTG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (10/12):  ref base (25/23):  total (35/35)
Fisher's exact test strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-01

GCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681020‑681217
                                                                                                    |                                                                                                 
gCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGcgcg                                                                                                   <  2:2754439/101‑1 (MQ=255)
 ccGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcga                                                                                                  >  2:2871165/1‑101 (MQ=255)
  ggggTGAGATGCTCGACCTCGATAGCCGCCAGCTGCGCCGTGCTCAGGGCCGCCACCTGCCGCGAGGTCATTGCCGCCTTCTAATCGCAGGTCAGCGctac                                                                                                 >  2:719188/2‑101 (MQ=255)
  cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 <  2:666620/101‑1 (MQ=255)
  cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 <  1:1611090/101‑1 (MQ=255)
  cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 <  2:1240284/101‑1 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                <  1:3827783/101‑1 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                >  2:2691507/1‑101 (MQ=255)
    ggTCAGATCCTCGACCTCGATCGCCGCCAGCTTCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTTAGAGcgacga                                                                                               >  2:2050643/1‑101 (MQ=255)
     gTCAGATGCTCGACCTTGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  1:719264/101‑1 (MQ=255)
     gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              >  1:2069690/1‑101 (MQ=255)
     gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              >  1:417726/1‑101 (MQ=255)
     gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  2:3689593/101‑1 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             <  1:4033477/101‑1 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             >  1:1981502/1‑101 (MQ=255)
       cAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGt                                                                                            <  2:162988/101‑1 (MQ=255)
        aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           >  2:3048576/1‑101 (MQ=255)
          aTCCTCGCCCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCCCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         <  1:2691713/101‑3 (MQ=255)
           tCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgcc                                                                                        >  2:4029812/1‑98 (MQ=255)
           tCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgcc                                                                                        <  1:3977231/101‑4 (MQ=255)
            cctcgatctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccg                                                                                       <  1:678183/101‑5 (MQ=255)
            cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccg                                                                                       <  1:2643631/101‑5 (MQ=255)
            cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccg                                                                                       >  1:220971/1‑97 (MQ=255)
            cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccg                                                                                       <  2:4201164/101‑5 (MQ=255)
             ctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccgg                                                                                      >  2:4070713/1‑96 (MQ=255)
              tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     >  1:845373/1‑101 (MQ=255)
              tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  2:118290/101‑1 (MQ=255)
              tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  2:429193/101‑1 (MQ=255)
               cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    <  2:1975304/101‑2 (MQ=255)
               cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    <  2:495185/101‑2 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  2:217306/101‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   >  1:3421741/1‑99 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  2:3024289/101‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   >  2:3622186/1‑99 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   >  2:526922/1‑99 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:115973/101‑3 (MQ=255)
                 ccctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  <  2:1436745/100‑4 (MQ=255)
                 acttcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  <  2:2498079/98‑4 (MQ=38)
                 acctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  <  1:3788679/101‑4 (MQ=255)
                  cctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgac                                                                                 <  2:4065647/101‑5 (MQ=255)
                  cctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgac                                                                                 <  1:980097/101‑5 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                >  2:2110448/1‑96 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                >  2:8417/1‑96 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                <  1:3437513/101‑6 (MQ=255)
                   ctcgaTCCCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                <  2:2392939/101‑6 (MQ=38)
                       aTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagccc                                                                            >  1:23078/1‑92 (MQ=37)
                        tCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcccc                                                                           >  1:1460837/1‑91 (MQ=38)
                                aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGc                                                                   >  1:3916708/1‑101 (MQ=25)
                                aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGc                                                                   >  1:1506451/1‑101 (MQ=25)
                                   tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                >  1:2103216/1‑98 (MQ=25)
                                   tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                >  1:3539066/1‑98 (MQ=25)
                                    gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                               >  2:3074367/1‑97 (MQ=25)
                                     cgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccac                                                              >  1:2113462/1‑96 (MQ=25)
                                     cgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccaa                                                              >  2:2617187/1‑96 (MQ=25)
                                      cccttGCTCAGCGCCGCCACCTGCTGCGAGGTTATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                             <  1:1625932/97‑7 (MQ=21)
                                       ccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGctccactg                                                            >  1:3933321/1‑94 (MQ=25)
                                       ccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactg                                                            >  1:3101068/1‑94 (MQ=25)
                                         gTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggc                                                          >  1:2777898/1‑92 (MQ=21)
                                                                                   tgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgt                <  1:503909/92‑1 (MQ=38)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  2:1253766/93‑1 (MQ=255)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  1:520833/93‑1 (MQ=255)
                                                                                     cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg              >  1:2463190/8‑101 (MQ=255)
                                                                                      gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt             >  1:2199396/7‑101 (MQ=255)
                                                                                       ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg            <  1:2791572/96‑1 (MQ=255)
                                                                                       ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg            <  2:1750863/96‑1 (MQ=255)
                                                                                        cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc           >  2:570188/5‑101 (MQ=255)
                                                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt          <  1:2684716/98‑1 (MQ=255)
                                                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt          >  2:4012875/4‑101 (MQ=255)
                                                                                         gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGCCGTGCt          >  1:3482095/4‑101 (MQ=255)
                                                                                         gtcgTCAGCGCCACGATGTCAGCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt          >  1:3118963/4‑101 (MQ=255)
                                                                                          tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  1:3109555/99‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:3173425/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:3768971/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  1:1727738/100‑1 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  2:386776/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  2:724654/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       <  2:2268405/101‑1 (MQ=255)
                                                                                             tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc      <  1:188795/101‑1 (MQ=255)
                                                                                             gcAGCGCCACGATGTCATCGGTGCCGAGCCCCTTGGTCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCGGCGTCGTGCTCAgc      <  2:2476925/100‑1 (MQ=38)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     >  1:2965703/1‑101 (MQ=255)
                                                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc    <  1:116740/101‑1 (MQ=255)
                                                                                               aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc    >  2:3045110/1‑101 (MQ=255)
                                                                                                gcgcGACGGTGTCATGGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGcgcg   <  2:1124803/101‑1 (MQ=255)
                                                                                                 cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCgg  <  2:1116459/101‑1 (MQ=255)
                                                                                                    |                                                                                                 
GCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681020‑681217

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