Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,120 G→C 34.1% 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→C34.1% 73.4 0.0 ‑0.979V829V (GTC→GTG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (17/10):  ref base (32/20):  total (49/30)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.13e-01

GCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681020‑681217
                                                                                                    |                                                                                                 
gCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGcgcg                                                                                                   <  2:2477040/101‑1 (MQ=255)
gCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGcgcg                                                                                                   <  1:1244851/101‑1 (MQ=255)
 ccGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcga                                                                                                  >  1:1429445/1‑101 (MQ=255)
 ccGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcga                                                                                                  >  1:1393510/1‑101 (MQ=255)
  cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 >  2:534997/1‑101 (MQ=255)
  cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 >  2:4214821/1‑101 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgccg                                                                                                >  2:4209964/1‑101 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                <  1:555343/101‑1 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                <  1:3114526/101‑1 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                <  2:589134/101‑1 (MQ=255)
   gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAACGcgacg                                                                                                >  1:49813/1‑101 (MQ=255)
    ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               >  2:4377847/1‑101 (MQ=255)
    ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               >  2:4710089/1‑101 (MQ=255)
     tgtgctcttccgatctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  1:679817/91‑1 (MQ=38)
     gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              >  2:4356446/1‑101 (MQ=255)
     gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              >  1:1984766/1‑101 (MQ=255)
     gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  1:4107553/101‑1 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             <  1:743167/101‑1 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             >  1:3820219/1‑101 (MQ=255)
      tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             >  2:2475323/1‑101 (MQ=255)
       cAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGt                                                                                            >  1:1223653/1‑101 (MQ=255)
       cAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGt                                                                                            >  2:1536234/1‑101 (MQ=255)
        aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           >  1:3239542/1‑101 (MQ=255)
        aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           <  1:1587771/101‑1 (MQ=255)
        aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           >  2:299230/1‑101 (MQ=255)
        aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGGCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           <  1:4872962/101‑1 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          <  1:3160850/101‑2 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  2:4317455/1‑100 (MQ=255)
         gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  1:829255/1‑100 (MQ=255)
          aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         >  1:1053164/1‑99 (MQ=255)
          aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         >  2:346368/1‑99 (MQ=255)
           tCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgcc                                                                                        <  2:149345/101‑4 (MQ=255)
             ctcgccctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccgg                                                                                      <  2:491389/101‑6 (MQ=255)
              tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  1:2724975/101‑1 (MQ=255)
               cgatctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    <  1:4816163/101‑2 (MQ=255)
               cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    >  1:987758/1‑100 (MQ=255)
               cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    <  2:1212861/101‑2 (MQ=255)
               cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    >  2:3062470/1‑100 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCTCCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATggcgccggtcg                                                                                   <  2:2008867/101‑11 (MQ=38)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:2825681/101‑3 (MQ=255)
                gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   >  2:2565710/1‑99 (MQ=255)
                cgaggagaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   >  2:537041/7‑99 (MQ=37)
                 acctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  <  1:911209/101‑4 (MQ=255)
                 acctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  >  2:1300650/1‑98 (MQ=255)
                 acctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcga                                                                                  <  1:4866358/101‑4 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                >  1:1281880/1‑96 (MQ=255)
                   ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                >  1:1003636/1‑96 (MQ=255)
                    tcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacag                                                                               <  2:1605346/101‑7 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCTGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  1:4500248/1‑94 (MQ=255)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGTCACTTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              <  1:4313600/101‑8 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  2:4616934/1‑94 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              <  1:31857/101‑8 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  1:1075353/1‑94 (MQ=38)
                     cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  2:1454509/1‑94 (MQ=38)
                      gaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcc                                                                             >  1:665238/1‑93 (MQ=37)
                       aTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagccc                                                                            >  1:4940896/1‑92 (MQ=37)
                       aTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagccc                                                                            <  1:3198744/101‑10 (MQ=37)
                        tCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcccc                                                                           <  2:2617247/101‑11 (MQ=37)
                                aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGc                                                                   >  2:2195651/1‑101 (MQ=25)
                                  cTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgc                                                                 >  1:667445/1‑99 (MQ=25)
                                    gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                               >  2:672420/1‑97 (MQ=25)
                                      gcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                             >  1:3001435/1‑95 (MQ=25)
                                        cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactgg                                                           >  2:169848/1‑93 (MQ=21)
                                        cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactgg                                                           >  1:3097594/1‑93 (MQ=21)
                                         gTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggc                                                          >  1:3820076/1‑92 (MQ=21)
                                                                                  ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  1:2438968/11‑101 (MQ=37)
                                                                                  ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 <  2:4796317/91‑1 (MQ=37)
                                                                                  ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  2:2906074/11‑101 (MQ=37)
                                                                                   tgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgt                >  1:3463498/10‑101 (MQ=37)
                                                                                    gcgccgtcgTCAGCGCCCCGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  2:715110/93‑1 (MQ=38)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               >  1:4375626/9‑101 (MQ=38)
                                                                                    gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  2:1906432/93‑1 (MQ=255)
                                                                                      gccgtcgTCAGCGCCACGATTTCATCGGTGCCGAGGCCCTTGATCGCATTGGCCCTCAGCGCCGTCACCCTGCCCGTGTTCAGCGCCGCCAGCGCcgtcgt             >  2:1286727/7‑101 (MQ=25)
                                                                                      gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt             >  2:462632/7‑101 (MQ=255)
                                                                                       ccgtcgTCAGCGCCACGAAGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg            <  2:3375965/96‑1 (MQ=255)
                                                                                        cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc           <  1:1225479/97‑1 (MQ=255)
                                                                                        cgccgTCGGCGCCACGATGTCA.CCGTTCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt          <  1:3041358/97‑1 (MQ=255)
                                                                                          tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  1:4748790/99‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:770652/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  1:1953516/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  2:3352622/100‑1 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  2:1853713/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        >  1:4193540/2‑101 (MQ=255)
                                                                                           cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATACAGCACCGCAAGCAGAGTCCTACTCa        >  2:3613137/2‑101 (MQ=37)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:1369276/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:22574/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  2:2876365/1‑101 (MQ=255)
                                                                                            gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:4521325/1‑101 (MQ=255)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     <  2:4322520/101‑1 (MQ=255)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     >  1:1361216/1‑101 (MQ=255)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     >  1:2567771/1‑101 (MQ=255)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     <  2:1138714/101‑1 (MQ=255)
                                                                                              cAGCGCCACGATGTCATCGGTGCCGAGGCCCGTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     >  1:4474570/1‑101 (MQ=255)
                                                                                                 cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCgg  >  1:1633281/1‑101 (MQ=255)
                                                                                                    |                                                                                                 
GCCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681020‑681217

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