Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,689 T→C 10.4% S640G (AGC→GGC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6890T→C10.4% 41.6 ‑3.1 ‑0.1183S640G (AGC→GGC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (12/7):  ref base (62/102):  total (74/109)
Fisher's exact test strand distribution p-value = 4.68e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.80e-01

Only 100 of 196 total aligned reads displayed.

GAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCGACAGC  >  NC_014034/681589‑681789
                                                                                                    |                                                                                                    
gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCt                                                                                                      <  2:949936/101‑1 (MQ=255)
 aaGCGAGCGCGGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc                                                                                                     <  1:851974/101‑1 (MQ=255)
 aaGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc                                                                                                     <  1:2125079/101‑1 (MQ=255)
  ggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                                                    >  1:2896334/7‑101 (MQ=2)
  ggcaatcgccgcAATTTGCGCCGCGCGCACCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                                                    <  1:1652321/95‑1 (MQ=2)
  ggcaaccgccgcGATCTGCGCCGCGCTCAGCGCTGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                                                    <  1:4306695/95‑1 (MQ=2)
  ggaattcgccgcGACCTGCGCCGCACGCCGCCCCGCCTTTTGCGCCGAGCCGAGCCCGCCGGCCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                                                    <  2:3624355/95‑1 (MQ=255)
  aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                                                    <  2:2340273/101‑1 (MQ=255)
  aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                                                    <  1:1129582/101‑1 (MQ=255)
  aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                                                    <  2:4819944/101‑1 (MQ=255)
   gcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCCCTGGCCAGCGCCGCGACCTGATCGGCGCCCAg                                                                                                   >  1:4915642/6‑101 (MQ=0)
   gcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAg                                                                                                   >  1:581391/6‑101 (MQ=2)
   gCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAg                                                                                                   <  2:143138/101‑1 (MQ=255)
    caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                                                  >  2:2829804/5‑100 (MQ=2)
      atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgatc                                                                                                >  2:4582967/3‑98 (MQ=2)
      atcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgact                                                                                                >  2:4315296/3‑98 (MQ=2)
      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                                                <  2:4112279/101‑1 (MQ=255)
      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                                                >  1:2514625/1‑101 (MQ=255)
      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCtt                                                                                                >  1:4393654/1‑101 (MQ=255)
       gcgcCGCGACCTGCGCCGTGCTCAGCCCCGCCGCCTGCACCGAGCCAAGCCCGTCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTg                                                                                               <  2:2514451/101‑1 (MQ=255)
        cgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactgg                                                                                              >  2:3882991/1‑96 (MQ=2)
         gccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGcc                                                                                             <  2:2942839/101‑1 (MQ=255)
         gccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGcc                                                                                             <  1:101973/101‑1 (MQ=255)
           cgcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGa                                                                                           <  1:4142581/101‑1 (MQ=255)
            gcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGAt                                                                                          >  2:2793756/1‑101 (MQ=255)
            gcgATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGAt                                                                                          <  2:2396388/101‑1 (MQ=255)
            gcgATCGGCGCGGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGAt                                                                                          <  2:4419512/101‑1 (MQ=255)
             cgATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggccagc                                                                                         >  2:4804953/1‑91 (MQ=2)
              gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                                                        >  1:1606620/1‑101 (MQ=255)
              gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                                                        <  1:3338349/101‑1 (MQ=255)
                tCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCTAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGc                                                                                      <  1:4313184/101‑1 (MQ=255)
                 cTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCg                                                                                     >  1:1990359/1‑101 (MQ=255)
                      ccGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTcgcg                                                                                >  2:4694665/1‑101 (MQ=255)
                        gTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCt                                                                              >  2:2656620/1‑101 (MQ=255)
                         tGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTc                                                                             <  2:4486441/101‑1 (MQ=255)
                         tGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTc                                                                             <  1:1817695/101‑1 (MQ=255)
                         tGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTc                                                                             <  2:2413160/101‑1 (MQ=255)
                         tGCTCAGCGCCGCCGCCGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTc                                                                             <  1:3795143/101‑1 (MQ=255)
                          gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa                                                                            <  2:3438006/101‑1 (MQ=255)
                          gCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCa                                                                            >  1:3302988/1‑101 (MQ=255)
                           cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                                                           >  1:1786314/1‑101 (MQ=255)
                           cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                                                           <  2:2064363/101‑1 (MQ=255)
                             cAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGc                                                                         <  1:1323889/101‑1 (MQ=255)
                             cAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGc                                                                         <  2:127941/101‑1 (MQ=255)
                              aGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGcc                                                                        >  1:4355711/1‑101 (MQ=255)
                              aGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGcc                                                                        <  1:1744376/101‑1 (MQ=255)
                                cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                                                      <  2:2719391/101‑1 (MQ=255)
                                    gccgccGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTGGCCGATCTGCGTCGCGCTCAGGCCCGAAAc                                                                  <  2:4010707/101‑1 (MQ=255)
                                     ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                                                 <  1:4726262/101‑1 (MQ=255)
                                       gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                                               >  1:3343477/1‑101 (MQ=255)
                                       gccGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                                               <  1:3455455/101‑1 (MQ=255)
                                        ccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgcc                                                              <  1:1406020/101‑1 (MQ=255)
                                         cTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgccg                                                             >  2:1053566/1‑101 (MQ=255)
                                          tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                                                            >  2:2460641/1‑101 (MQ=255)
                                             aCCGAGCCAACCCCGGCCACCGCAAGGATGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGAGTTGCGCTCAGGCCCGCAACCGCCGCGCt                                                         >  2:3264789/1‑101 (MQ=255)
                                              ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                                                        <  2:1965186/101‑1 (MQ=255)
                                               cGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGc                                                       <  1:3520078/101‑1 (MQ=255)
                                                gAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGcc                                                      >  2:878278/1‑101 (MQ=255)
                                                 aGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGccc                                                     <  2:4064369/101‑1 (MQ=255)
                                                  gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                                                    <  2:4762603/101‑1 (MQ=255)
                                                   ccAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAg                                                   <  1:3247082/101‑1 (MQ=255)
                                                       gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                                               <  1:3192759/101‑1 (MQ=255)
                                                       gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                                               <  1:3395998/101‑1 (MQ=255)
                                                       gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                                               <  2:1862658/101‑1 (MQ=255)
                                                        cccGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCcg                                              <  2:1242222/101‑1 (MQ=255)
                                                           ggCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGa                                           >  1:4648543/1‑101 (MQ=255)
                                                            gCCACCGCAAGGCTCGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGAt                                          <  2:3343214/101‑1 (MQ=255)
                                                             ccACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATc                                         <  2:2908793/101‑1 (MQ=255)
                                                              cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt                                        >  1:1033327/1‑101 (MQ=255)
                                                              cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCCCCGCGCTGCCCAGCCCCGCGATCt                                        <  1:26759/101‑1 (MQ=255)
                                                               aCCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTg                                       >  2:228153/1‑101 (MQ=255)
                                                                   cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc                                   >  2:3042288/1‑101 (MQ=255)
                                                                   cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc                                   >  2:548913/1‑101 (MQ=255)
                                                                   cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc                                   >  2:1597027/1‑101 (MQ=255)
                                                                    aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCg                                  <  1:4985442/101‑1 (MQ=255)
                                                                    aaGGCTGGAGAGCGCGCCGACCTGATCGGGGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACGGACGCGCTGCCCAGCCCCGCGATGTGAGCCg                                  >  1:2909016/1‑101 (MQ=255)
                                                                      ggCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTg                                <  2:1529150/101‑1 (MQ=255)
                                                                      ggCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTg                                <  1:2793958/101‑1 (MQ=255)
                                                                        gtGGAGACCACGCCGCCCTGAAAGGTGCCCAGCTTGCGGCTCTGCGGCGGGCGCAGGCCAGAAACCGCCGCGCTGCCGAGCCCCGCGATCTGAGCCGTGcc                              <  2:2002849/100‑1 (MQ=255)
                                                                         tGGAGAGGGCGACGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCCCTGAGGCACGAAACGGCCGGGCTGGCCACCCCGGCGACATGCGCCGTGccc                             >  2:3359548/1‑101 (MQ=255)
                                                                         ggcagagCCCGACGACCTGTTCGTTGCTCAGCTTGCCGGTCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTTGCCAGCCCCGCGATCTGCGCCGTGccc                             <  2:626952/98‑1 (MQ=37)
                                                                          ggAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCa                            <  2:1786181/101‑1 (MQ=255)
                                                                          ggAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCa                            <  1:3868318/101‑1 (MQ=255)
                                                                          ggAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCa                            >  1:4855343/1‑101 (MQ=255)
                                                                          ggAGAGCGCGCCGACCTGAGGGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCa                            <  1:2996460/101‑1 (MQ=255)
                                                                           gagaGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGccccg                           >  2:2999840/1‑99 (MQ=255)
                                                                            agagCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGc                          >  1:1574237/1‑101 (MQ=255)
                                                                              agCGCGCCGACCTGATCAGTGCTCCGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGt                        <  1:2460815/101‑1 (MQ=255)
                                                                                  cgcCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTTCCCAGCGTCGCa                    <  1:3609268/101‑1 (MQ=255)
                                                                                     cGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGc                 <  1:1597159/101‑1 (MQ=255)
                                                                                          tGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTgcgc            >  2:3846402/1‑101 (MQ=255)
                                                                                           gATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGcc           <  1:514008/101‑1 (MQ=255)
                                                                                             tCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGt         <  2:4294935/101‑1 (MQ=255)
                                                                                             tCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGt         >  2:3132127/1‑101 (MQ=255)
                                                                                               ggTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCg       >  1:2382318/1‑101 (MQ=255)
                                                                                               ggTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCg       <  2:1234143/101‑1 (MQ=255)
                                                                                               ggTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCg       >  1:2241045/1‑101 (MQ=255)
                                                                                                  gCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCGACa    >  2:1666916/1‑101 (MQ=255)
                                                                                                   cTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCGACAg   >  2:3991756/1‑101 (MQ=255)
                                                                                                   cTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCGACAg   >  1:3292379/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGTGCCCAGCGTCGCAAGCTGCGCCGTCGACAGC  >  NC_014034/681589‑681789

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