Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,120 G→C 40.6% 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→C40.6% 87.0 ‑0.5 ‑0.274V829V (GTC→GTG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (14/16):  ref base (24/20):  total (38/36)
Fisher's exact test strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.62e-01

CCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681021‑681217
                                                                                                   |                                                                                                 
ccGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcga                                                                                                  >  1:1022712/1‑101 (MQ=255)
ccGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcga                                                                                                  <  2:2167112/101‑1 (MQ=255)
 cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 >  1:424716/1‑101 (MQ=255)
 cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 <  2:924064/101‑1 (MQ=255)
 cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 >  2:920855/1‑101 (MQ=255)
 cGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgac                                                                                                 >  2:3840686/1‑101 (MQ=255)
  gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                >  2:3532719/1‑101 (MQ=255)
  gggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacg                                                                                                >  2:335293/1‑101 (MQ=255)
   ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCCCCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               <  2:4445498/101‑1 (MQ=255)
   ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGcgacga                                                                                               <  1:1605449/101‑1 (MQ=255)
   ggTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTAAGCGcgacga                                                                                               >  2:1973808/1‑101 (MQ=255)
    gTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  1:1950002/101‑1 (MQ=255)
    gTCAGATCCCCGACCTCCATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCAGCGCCGCCATCTGATCGCCGGTCAGCGCGACGAt                                                                                              <  2:3648614/101‑1 (MQ=255)
     tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             <  1:3774430/101‑1 (MQ=255)
     tCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATg                                                                                             >  1:922483/1‑101 (MQ=255)
      cAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGGTGt                                                                                            <  2:4240730/101‑1 (MQ=255)
       aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           >  1:3391759/1‑101 (MQ=255)
       aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           >  1:1619183/1‑101 (MQ=255)
       aGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTc                                                                                           <  2:752972/101‑1 (MQ=255)
        gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  2:2714463/1‑100 (MQ=255)
        gATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcg                                                                                          >  2:990562/1‑100 (MQ=255)
         aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         <  2:3847130/101‑3 (MQ=255)
         aTCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgc                                                                                         >  2:1031151/1‑99 (MQ=255)
          tCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgcc                                                                                        <  1:3191195/101‑4 (MQ=255)
           cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccg                                                                                       <  2:3846374/101‑5 (MQ=255)
           cctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGAGGTcgccg                                                                                       >  2:2098926/1‑97 (MQ=255)
            ctcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTcgccgg                                                                                      <  1:2067887/101‑6 (MQ=255)
             tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     >  1:2070060/1‑101 (MQ=255)
             tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  2:2130884/101‑1 (MQ=255)
             tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGt                                                                                     <  1:577544/101‑1 (MQ=255)
             tcgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGCCAGCGCGACGATGTCGCCGGt                                                                                     >  2:3752015/1‑101 (MQ=255)
              cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    <  1:2361815/101‑2 (MQ=255)
              cgacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtc                                                                                    <  2:766218/101‑2 (MQ=255)
               gacctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcg                                                                                   <  1:563625/101‑3 (MQ=255)
                 cctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgac                                                                                 <  1:2013672/101‑5 (MQ=255)
                  ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCGATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                <  2:1512624/101‑6 (MQ=38)
                  ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                >  1:3479584/1‑96 (MQ=255)
                  ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                <  1:668132/101‑6 (MQ=255)
                  ctcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgaca                                                                                >  2:447106/1‑96 (MQ=255)
                   tcgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacag                                                                               <  1:1552349/101‑7 (MQ=38)
                    cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  1:952811/1‑94 (MQ=38)
                    cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              >  1:1000865/1‑94 (MQ=38)
                    cgaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagc                                                                              <  1:4589884/101‑8 (MQ=38)
                     gaTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagcc                                                                             >  1:2311853/1‑93 (MQ=37)
                      aTCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGtcgacagccc                                                                            >  2:4377718/1‑92 (MQ=37)
                                gCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcg                                                                  >  2:1208874/1‑100 (MQ=25)
                                 cTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGTCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgc                                                                 >  2:1793744/1‑99 (MQ=25)
                                 cTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgc                                                                 >  2:4287412/1‑99 (MQ=25)
                                  tGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcc                                                                >  1:2917857/1‑98 (MQ=25)
                                   gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                               >  1:3712510/1‑97 (MQ=25)
                                   gcgcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgcca                                                               >  2:1730515/1‑97 (MQ=25)
                                     gcCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                             >  1:1237129/1‑95 (MQ=25)
                                     cccGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccact                                                             <  1:829920/100‑7 (MQ=21)
                                       cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGGCGGACAGCCCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccaccgg                                                           >  2:3113376/1‑93 (MQ=16)
                                       cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactgg                                                           >  1:4219663/1‑93 (MQ=21)
                                        gTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggc                                                          >  1:1619573/1‑92 (MQ=21)
                                         tGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggcc                                                         >  2:1414237/1‑91 (MQ=255)
                                         tGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGCGcgccactggcc                                                         >  1:4254928/1‑91 (MQ=25)
                                                                                 ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGTCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  1:3921526/11‑101 (MQ=37)
                                                                                 ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  2:1736550/11‑101 (MQ=37)
                                                                                 ctgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcg                 >  2:2361205/11‑101 (MQ=37)
                                                                                  tgcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgt                >  2:2230505/10‑101 (MQ=38)
                                                                                   gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  1:1556247/93‑1 (MQ=255)
                                                                                   gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  2:3294009/93‑1 (MQ=255)
                                                                                   gcgccgtcgTCAGCGCCACGATGTCATCGGTGCCCAGCCCCTTGATCGCATTGGCACCCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtc               <  1:3987090/93‑1 (MQ=37)
                                                                                    cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg              <  1:4575800/94‑1 (MQ=255)
                                                                                    cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg              >  1:1154046/8‑101 (MQ=255)
                                                                                    cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg              >  1:1892995/8‑101 (MQ=255)
                                                                                    cgccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcg              >  2:4181603/8‑101 (MQ=255)
                                                                                     gccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGcgtcgt             <  1:1613432/95‑1 (MQ=255)
                                                                                      ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg            <  1:2105024/96‑1 (MQ=255)
                                                                                      ccgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCAATGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTg            >  2:4037655/6‑101 (MQ=38)
                                                                                       cgtcgttaGCGCCACGATGTCATCGGTGCCGATGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc           <  2:2776385/94‑1 (MQ=38)
                                                                                       cgtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGc           <  1:3260145/97‑1 (MQ=255)
                                                                                        gtcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCt          >  1:768814/4‑101 (MQ=255)
                                                                                         tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         >  1:2245014/3‑101 (MQ=255)
                                                                                         tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  1:2048137/99‑1 (MQ=255)
                                                                                         tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  1:193633/99‑1 (MQ=255)
                                                                                         tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  2:3443084/99‑1 (MQ=255)
                                                                                         tcgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTc         <  2:2264202/99‑1 (MQ=255)
                                                                                          cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  1:2330521/100‑1 (MQ=255)
                                                                                          cgTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCa        <  1:3194457/100‑1 (MQ=255)
                                                                                           gTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       >  1:1785196/1‑101 (MQ=255)
                                                                                           ctCAGCGCCACGATGTCATCGGTGCCGAGGGCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAg       <  2:3315582/100‑1 (MQ=255)
                                                                                            tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc      >  1:2285044/1‑101 (MQ=255)
                                                                                            tCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgc      >  2:245392/1‑101 (MQ=255)
                                                                                             cAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcg     >  1:2296602/1‑101 (MQ=255)
                                                                                              aGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAgcgc    <  2:2748494/101‑1 (MQ=255)
                                                                                               gcgcCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGcgcg   <  2:3760670/101‑1 (MQ=255)
                                                                                                cgccACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCgg  <  1:2384733/101‑1 (MQ=255)
                                                                                                   |                                                                                                 
CCGGGTCAGATCCTCGACCTCGATCGCCGCCAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCATCTGATCGCCGGTCAGCGCGACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTGCGTCGTGCTCAGCGCGG  >  NC_014034/681021‑681217

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