Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,459 G→A 20.1% I716I (ATC→ATT RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4590G→A20.1% 99.5 ‑2.6 ‑0.0189I716I (ATC→ATTRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (13/25):  ref base (77/74):  total (90/99)
Fisher's exact test strand distribution p-value = 7.13e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.34e-01

Only 100 of 217 total aligned reads displayed.

CCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCG  >  NC_014034/681359‑681559
                                                                                                    |                                                                                                     
ccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCcgcg                                                                                                       <  1:841572/101‑1 (MQ=255)
ccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCcgcg                                                                                                       >  1:44751/1‑101 (MQ=255)
  gcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGAt                                                                                                     <  2:2338672/101‑1 (MQ=255)
   cgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATc                                                                                                    >  2:1077753/1‑101 (MQ=255)
    gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                                   <  1:2171344/101‑1 (MQ=255)
     aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                                  >  1:73075/1‑101 (MQ=255)
      tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                                 >  2:2598380/1‑101 (MQ=255)
       cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                                >  1:4219819/1‑101 (MQ=255)
       cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                                <  2:4227612/101‑1 (MQ=255)
        tGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGt                                                                                               <  1:976342/101‑1 (MQ=255)
        tGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGt                                                                                               >  2:2182513/1‑101 (MQ=255)
        tGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGt                                                                                               <  1:3794412/101‑1 (MQ=255)
           gcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                                            <  2:1305946/101‑1 (MQ=255)
               cTGAAAGCCCGGCCAGCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                                        <  2:4114890/101‑1 (MQ=255)
                tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                       >  1:112510/1‑101 (MQ=255)
                tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                       >  1:1736606/1‑101 (MQ=255)
                 gAAGGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                                      <  2:161796/101‑1 (MQ=255)
                     gCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgcg                                                                                  >  2:1644431/1‑101 (MQ=255)
                      cccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgcgc                                                                                 <  2:4406255/101‑1 (MQ=255)
                       ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGcc                                                                                >  1:219222/1‑101 (MQ=255)
                        cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                                               >  1:758789/1‑101 (MQ=255)
                         ggCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.a                                                                             >  1:2091647/1‑101 (MQ=255)
                         ggCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.a                                                                             >  2:1849205/1‑101 (MQ=255)
                          gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.At                                                                            <  1:1104107/101‑1 (MQ=255)
                           cactCGTGCCCGAGCCCAGCCCGGAGCGCGCCCTCGCCCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCCCGCCG.ATc                                                                           <  2:2912958/98‑1 (MQ=37)
                            cATCGTGGCCGAGCCGCGCCCGGAGGGGGCCGTCGCGCTGATCTGCGCCAGATCCACCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCACCG.ATCt                                                                          <  2:2193460/101‑1 (MQ=255)
                            cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                          <  2:874939/101‑1 (MQ=255)
                            cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                          <  2:73062/101‑1 (MQ=255)
                             aTCGTGGCCTGGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                         <  2:2215312/101‑1 (MQ=255)
                              tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgc                                                                        <  2:3183806/101‑1 (MQ=255)
                                 tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGcc                                                                     >  1:2796337/1‑101 (MQ=255)
                                    ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAg                                                                  <  1:225089/101‑1 (MQ=255)
                                    ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAg                                                                  >  2:1198827/1‑101 (MQ=255)
                                      gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AGCTGCGCCGAGcc                                                                >  2:2024083/1‑101 (MQ=255)
                                      gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCCTGCTGGTCCGCGCGCCG.ATCTGCGCCGAGcc                                                                >  1:3056405/1‑101 (MQ=255)
                                       agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGccc                                                               >  2:1151173/1‑101 (MQ=255)
                                          ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGa                                                            <  2:2645941/101‑1 (MQ=255)
                                          ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGa                                                            >  2:3122573/1‑101 (MQ=255)
                                          ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGa                                                            >  2:43273/1‑101 (MQ=255)
                                           cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAc                                                           <  1:1376004/101‑1 (MQ=255)
                                              cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcgcg                                                        <  2:2830636/101‑1 (MQ=255)
                                               ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACgcgc                                                       <  1:2353246/101‑1 (MQ=255)
                                               ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACgcgc                                                       >  1:2029343/1‑101 (MQ=255)
                                                  gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgccgc                                                    >  1:1385939/1‑101 (MQ=255)
                                                    gagCGCCGGCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                  <  1:197067/101‑1 (MQ=255)
                                                     agCGCCGTCGCCCTGATCTGCGCCAGATCCACCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTg                                                 <  2:2125952/101‑1 (MQ=255)
                                                       cgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTgcg                                               <  2:3588081/101‑1 (MQ=255)
                                                        gcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTgcgc                                              >  2:227220/1‑101 (MQ=255)
                                                          cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                            >  1:2873829/1‑101 (MQ=255)
                                                             cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                         >  2:2523639/1‑101 (MQ=255)
                                                             cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                                         >  1:1057054/1‑101 (MQ=255)
                                                              gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCACAGACGCGCCGCCTGCGCcgtcgt                                        >  1:1750279/1‑101 (MQ=255)
                                                               cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                       <  2:2726811/101‑1 (MQ=255)
                                                                 gggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAg                                     <  1:97186/91‑1 (MQ=2)
                                                                  ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                    <  2:4043686/92‑1 (MQ=2)
                                                                   gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                                   <  1:123826/93‑1 (MQ=2)
                                                                   gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                                   <  2:1592623/93‑1 (MQ=2)
                                                                   gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                   <  1:158849/101‑1 (MQ=255)
                                                                   gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                   <  2:3263145/101‑1 (MQ=255)
                                                                   gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                   >  1:68898/1‑101 (MQ=255)
                                                                   gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                   <  1:2598551/101‑1 (MQ=255)
                                                                   gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                   >  2:1361671/1‑101 (MQ=255)
                                                                     tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                 <  1:2002385/101‑1 (MQ=255)
                                                                     tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                 <  2:153344/101‑1 (MQ=255)
                                                                      ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGccc                                >  2:3550824/6‑100 (MQ=2)
                                                                       tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                               <  2:14506/101‑1 (MQ=255)
                                                                        gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                              <  1:3196885/98‑1 (MQ=2)
                                                                        gcaccAGATCGGCGGCTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCCTCGACCCCAGCGCCGCCGCCTGCGCCGTCGTCACGGCCACg                              >  1:3295658/4‑101 (MQ=255)
                                                                          cccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                            <  2:248389/100‑1 (MQ=9)
                                                                          accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                            <  2:704030/100‑1 (MQ=9)
                                                                           ccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                           <  1:2497133/101‑1 (MQ=2)
                                                                           ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                           <  1:570036/101‑1 (MQ=255)
                                                                            cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                          >  1:2054818/1‑101 (MQ=255)
                                                                               aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                       <  1:160332/101‑1 (MQ=255)
                                                                                  ccccgccCCGGGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                    >  2:3673217/3‑101 (MQ=255)
                                                                                   gcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                   <  1:3699636/92‑1 (MQ=9)
                                                                                   gcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                   >  1:3979609/10‑101 (MQ=9)
                                                                                   gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                   >  2:4058028/1‑101 (MQ=255)
                                                                                    cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                  <  1:1512639/93‑1 (MQ=9)
                                                                                     ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                 <  1:3033885/94‑1 (MQ=9)
                                                                                      gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                <  2:2042354/101‑1 (MQ=255)
                                                                                       tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg               >  1:263843/6‑100 (MQ=9)
                                                                                       tttccaGCCCCGCAACCTGCGTGCGCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccg               <  2:3043888/96‑2 (MQ=9)
                                                                                       ccgcGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa               <  1:1648844/98‑1 (MQ=255)
                                                                                         tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa             >  1:940213/1‑101 (MQ=255)
                                                                                           gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc           >  2:1917211/1‑101 (MQ=255)
                                                                                           caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGGCATCGGTGccgaggc           <  1:3750074/100‑6 (MQ=9)
                                                                                            aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc          >  2:2289594/1‑101 (MQ=255)
                                                                                            aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc          >  1:3504570/1‑101 (MQ=2)
                                                                                             gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc         <  1:2208014/101‑2 (MQ=9)
                                                                                              cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct        >  2:4282110/1‑99 (MQ=2)
                                                                                               gccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg       >  2:1293045/1‑101 (MQ=255)
                                                                                               gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt       <  2:3937978/101‑4 (MQ=2)
                                                                                               gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt       <  1:3082810/101‑4 (MQ=9)
                                                                                               gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt       <  2:798942/101‑4 (MQ=9)
                                                                                                  gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt    <  1:4334650/101‑1 (MQ=255)
                                                                                                  gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt    <  1:1917335/101‑1 (MQ=255)
                                                                                                  gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt    >  1:4170241/1‑101 (MQ=255)
                                                                                                    gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg  <  1:3220519/101‑1 (MQ=255)
                                                                                                    gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg  >  1:654593/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
CCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCG  >  NC_014034/681359‑681559

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