Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,541 T→C 7.8% M356M (ATG→GTG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034682,5410T→C 7.8% 29.4 ‑0.5 ‑0.8191M356M (ATG→GTG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (6/9):  ref base (86/90):  total (92/99)
Fisher's exact test strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.48e-01

Only 100 of 205 total aligned reads displayed.

TCGATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCC  >  NC_014034/682441‑682641
                                                                                                    |                                                                                                     
tCGATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCAt                                                                                                       <  2:4262919/101‑1 (MQ=255)
 cGATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATc                                                                                                      >  1:2136681/1‑101 (MQ=255)
 cGATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATc                                                                                                      >  1:1792604/1‑101 (MQ=255)
  gATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCt                                                                                                     >  2:3632300/1‑101 (MQ=255)
  gATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCt                                                                                                     >  1:4711629/1‑101 (MQ=255)
  gATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCt                                                                                                     >  2:1402508/1‑101 (MQ=255)
   aTCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTg                                                                                                    >  1:4142656/1‑101 (MQ=255)
    tCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTgc                                                                                                   <  1:5034165/101‑1 (MQ=255)
           cGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCgg                                                                                            <  2:1689232/101‑1 (MQ=255)
           cGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCgg                                                                                            >  2:897385/1‑101 (MQ=255)
            gAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGTCCGTGCTCAGCGCCGCCATCTGCGCCCCGGt                                                                                           >  1:1514241/1‑101 (MQ=255)
            gAACTCAGCGCCTCGGCGGCCTCCTCGGACAGAGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGt                                                                                           >  1:1414022/1‑101 (MQ=255)
             aaCTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTc                                                                                          <  2:4125258/101‑1 (MQ=255)
             aaCTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTc                                                                                          <  1:2179143/101‑1 (MQ=255)
             aaCTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTc                                                                                          <  1:4204822/101‑1 (MQ=255)
                tCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAgc                                                                                       <  1:951539/101‑1 (MQ=255)
                   gcgcCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGcgcc                                                                                    >  2:4889524/1‑101 (MQ=255)
                    cgcCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGcgccg                                                                                   <  1:1206750/101‑1 (MQ=255)
                      ccTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCgccgcc                                                                                 >  1:1320638/1‑101 (MQ=255)
                      ccTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCgccgcc                                                                                 >  2:3921250/1‑101 (MQ=255)
                         cggcggCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGt                                                                              >  2:629815/1‑101 (MQ=255)
                         cggcggCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGt                                                                              >  2:2765746/1‑101 (MQ=255)
                          ggcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGtt                                                                             <  2:3536939/101‑1 (MQ=255)
                           gcggcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTc                                                                            >  1:3358585/1‑101 (MQ=255)
                              gcCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGAt                                                                         <  1:3842543/101‑1 (MQ=255)
                              gcCTCCGCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGAt                                                                         <  1:798382/101‑1 (MQ=255)
                               cctcctCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATc                                                                        <  2:3686922/101‑1 (MQ=255)
                                ctcctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCg                                                                       >  1:3291258/1‑101 (MQ=255)
                                  cctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.Tg                                                                    <  1:1263303/101‑1 (MQ=255)
                                   ctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGc                                                                   >  2:1357301/1‑101 (MQ=255)
                                   ctcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGc                                                                   <  2:3743538/101‑1 (MQ=255)
                                    tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                  >  2:3166475/1‑101 (MQ=255)
                                    tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                  >  1:1344043/1‑101 (MQ=255)
                                    tcGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCt                                                                  >  1:4642406/1‑101 (MQ=255)
                                      ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                >  1:3147225/1‑101 (MQ=255)
                                      ggACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCa                                                                <  2:3817032/101‑1 (MQ=255)
                                      gcaCAGCGCCGCGATCTACCGCGTGCCCAGAACCCCCAACTGCCCCGTCCTCACGGCCGCCATCTGCGCCCCCGTCTCCCCCGCCATTCCCATCG.TGCTcg                                                                >  2:2855768/3‑100 (MQ=255)
                                         cAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGcc                                                             <  2:2375135/101‑1 (MQ=255)
                                          aGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGccc                                                            <  1:2995205/101‑1 (MQ=255)
                                             gccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGc                                                         >  2:536293/1‑101 (MQ=255)
                                             gccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGc                                                         >  2:2639325/1‑101 (MQ=255)
                                              ccgcGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGcc                                                        <  1:4172715/101‑1 (MQ=255)
                                                gcgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCAt                                                      >  1:2054857/1‑101 (MQ=255)
                                                 cgATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATg                                                     <  2:4025890/101‑1 (MQ=255)
                                                   aTCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTc                                                   <  1:2473249/101‑1 (MQ=255)
                                                   aTCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTc                                                   <  1:2564807/101‑1 (MQ=255)
                                                    tCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTcc                                                  >  2:2849372/1‑101 (MQ=255)
                                                     cTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCt                                                 >  2:305032/1‑101 (MQ=255)
                                                         cgcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                             <  1:1128282/101‑1 (MQ=255)
                                                          gcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                            >  1:4530633/1‑101 (MQ=255)
                                                          gcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                            <  2:1432816/101‑1 (MQ=255)
                                                            gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                          >  1:579468/1‑101 (MQ=255)
                                                            gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                          <  1:2659120/101‑1 (MQ=255)
                                                             tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                         <  2:184440/101‑1 (MQ=255)
                                                             tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                         >  1:1021938/1‑101 (MQ=255)
                                                              gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                        >  2:473556/1‑101 (MQ=255)
                                                              gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                        <  1:2765948/101‑1 (MQ=255)
                                                              gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                        >  1:713118/1‑101 (MQ=255)
                                                               cccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAg                                       <  1:2615838/101‑1 (MQ=255)
                                                                ccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGc                                      >  2:474960/1‑101 (MQ=255)
                                                                  aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                    >  1:2839268/1‑101 (MQ=255)
                                                                  aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGGCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                    >  1:1940018/1‑101 (MQ=255)
                                                                   gAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGcccc                                   <  2:2071888/101‑1 (MQ=255)
                                                                      ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                <  1:3083497/98‑2 (MQ=37)
                                                                      cccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                <  1:3309976/101‑1 (MQ=255)
                                                                      cccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTc                                <  1:1834894/101‑1 (MQ=255)
                                                                       cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                               <  1:2314078/99‑3 (MQ=37)
                                                                       cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                               >  2:1132245/3‑99 (MQ=37)
                                                                       cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                               >  1:3597997/3‑99 (MQ=37)
                                                                       ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                               <  2:4642083/101‑1 (MQ=255)
                                                                       ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                               <  2:1320517/101‑1 (MQ=255)
                                                                       ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TG..CAGCTCGGTCATGTCATAGATCGGAAgagcacacgtc                             >  1:432844/1‑91 (MQ=25)
                                                                        cccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAg                              <  2:2485230/101‑1 (MQ=255)
                                                                         ccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgc                             <  1:1611283/101‑1 (MQ=25)
                                                                          cAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcg                            >  1:853404/1‑101 (MQ=25)
                                                                           aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                           >  1:2389598/1‑101 (MQ=255)
                                                                           aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcgc                           <  2:2466930/101‑1 (MQ=25)
                                                                             cTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGcgcgc                         >  1:4754103/1‑100 (MQ=255)
                                                                              tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                        >  2:4221088/1‑101 (MQ=255)
                                                                              tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                        >  2:4998914/1‑101 (MQ=255)
                                                                               gggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggcg                       <  2:738270/101‑1 (MQ=255)
                                                                                 gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                     <  2:4358505/101‑6 (MQ=21)
                                                                                  ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                    <  2:2592166/101‑1 (MQ=255)
                                                                                   cGTGCTCGGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGc                   <  1:3445179/101‑1 (MQ=255)
                                                                                   cGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccactgg                   <  2:1053054/101‑8 (MQ=25)
                                                                                    gTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGcc                  >  2:4419686/1‑101 (MQ=255)
                                                                                      gCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGa                >  2:2326604/1‑101 (MQ=255)
                                                                                        tCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc              <  1:1774271/101‑1 (MQ=255)
                                                                                        tCAGCGCCGCCATCTGCCCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc              <  2:1588908/101‑1 (MQ=255)
                                                                                         cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg             >  2:3056828/1‑101 (MQ=255)
                                                                                          aGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgc            >  2:3188072/1‑101 (MQ=255)
                                                                                           gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc           <  2:3870579/101‑1 (MQ=255)
                                                                                           gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc           <  2:4579640/101‑1 (MQ=255)
                                                                                           gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCACCCCGGCCATGTCCTCGGTCGAAAGCCCGTTCAGCGCGGCGGTGCCGATcccc           >  2:297821/1‑101 (MQ=255)
                                                                                              ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa        <  1:1590026/101‑1 (MQ=255)
                                                                                               cgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaa       >  2:1482470/1‑101 (MQ=255)
                                                                                               cgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaa       <  1:3753910/101‑1 (MQ=255)
                                                                                                gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa      <  2:4210620/101‑1 (MQ=255)
                                                                                                gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa      >  1:3684156/1‑101 (MQ=255)
                                                                                                    tCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGcc  <  2:4470436/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
TCGATCTGCGCCGAACTCAGCGCCTCGGCGGCCTCCTCGGACAGCGCCGCGATCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCC  >  NC_014034/682441‑682641

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