Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,594 T→C 9.5% E338G (GAG→GGG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034682,5940T→C 9.5% 43.3 ‑3.6 ‑1.9200E338G (GAG→GGG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (5/14):  ref base (99/82):  total (104/96)
Fisher's exact test strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.39e-02

Only 100 of 208 total aligned reads displayed.

CTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGCCGA  >  NC_014034/682494‑682694
                                                                                                     |                                                                                                    
cTGTCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCt                                                                                                      <  1:3909966/101‑1 (MQ=255)
cTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCt                                                                                                      >  1:1353116/1‑101 (MQ=255)
    cgcgTGCCCAGACCCCCCAGCGGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                                                  <  2:2974154/101‑1 (MQ=255)
        tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                                                              >  2:3533789/1‑101 (MQ=255)
         gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                                             <  2:718807/101‑1 (MQ=255)
          cccAGACCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAg                                                                                            <  1:3144448/101‑1 (MQ=255)
             aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                                                         >  2:2824666/1‑101 (MQ=255)
              gAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGcccc                                                                                        >  1:1384239/1‑101 (MQ=255)
               cgccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCt                                                                                       <  1:4392785/96‑1 (MQ=38)
               cgccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCt                                                                                       >  1:1885415/6‑101 (MQ=37)
                gccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCtt                                                                                      <  2:2782014/97‑1 (MQ=37)
                gccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCtt                                                                                      <  1:3296719/97‑1 (MQ=37)
                gccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCtt                                                                                      <  2:424278/97‑1 (MQ=37)
                aCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCtt                                                                                      <  2:4052756/101‑1 (MQ=255)
                aCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCtt                                                                                      >  2:2837347/1‑101 (MQ=255)
                 ccgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttg                                                                                     <  2:2600036/98‑2 (MQ=37)
                  ccccAGCGGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                                                    <  2:2165414/101‑1 (MQ=255)
                   gccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAg                                                                                   >  1:233123/2‑101 (MQ=37)
                    ccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGGTCCCGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgc                                                                                  <  1:4060728/101‑1 (MQ=255)
                        cTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcgg                                                                              >  1:512498/1‑101 (MQ=255)
                        cTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCACCTTCAGcgcgc                                                                              >  2:1273313/1‑100 (MQ=255)
                         tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                                             >  2:466941/1‑101 (MQ=255)
                         tGGCCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                                             <  2:4061720/101‑1 (MQ=255)
                          gcgCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgcca                                                                            <  1:3391529/99‑4 (MQ=25)
                           ggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGcggcgg                                                                           >  1:4254797/1‑101 (MQ=255)
                           cgCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccac                                                                           <  2:244022/100‑5 (MQ=25)
                            gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                                          <  1:4522819/101‑6 (MQ=21)
                             ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                                                         <  1:36113/101‑1 (MQ=255)
                             ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                                                         >  1:4587206/1‑101 (MQ=255)
                                tGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCg                                                                      <  1:391503/101‑1 (MQ=255)
                                  cTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGAt                                                                    >  2:2297673/1‑101 (MQ=255)
                                  cTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGAt                                                                    >  1:1049712/1‑101 (MQ=255)
                                   tCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                                                                   >  2:3146290/1‑101 (MQ=255)
                                    cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                                                  >  2:4446391/1‑101 (MQ=255)
                                    cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                                                  <  2:3008287/101‑1 (MQ=255)
                                    cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                                                  <  2:2467621/101‑1 (MQ=255)
                                     aGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgc                                                                 >  2:641612/1‑101 (MQ=255)
                                        gccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccgc                                                              >  1:1959417/1‑101 (MQ=255)
                                           gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                                           <  1:824250/101‑1 (MQ=255)
                                           gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                                           <  1:2592249/101‑1 (MQ=255)
                                            ccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAg                                                          <  2:3472299/101‑1 (MQ=255)
                                             cATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAgg                                                         <  1:2860437/101‑1 (MQ=255)
                                                cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                                                      <  1:1666574/101‑1 (MQ=255)
                                                 tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc                                                     >  2:2019194/1‑101 (MQ=255)
                                                 tGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcgc                                                     <  1:3100073/101‑1 (MQ=255)
                                                   cgcCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCgcgc                                                   >  2:3843941/1‑101 (MQ=255)
                                                    gcCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGcc                                                  <  1:2387627/101‑1 (MQ=255)
                                                    gcCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGcc                                                  <  2:512447/101‑1 (MQ=255)
                                                      cccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCgg                                                <  2:3477940/101‑1 (MQ=255)
                                                        cGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaa                                              <  1:1591513/101‑1 (MQ=255)
                                                         ggccAGCGCCGCCATTGCGGTCG.TGCTCAGCCCGGCCATGTCCTCGGGCGAAAGCCCCTTGAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaa                                             <  1:1273404/98‑1 (MQ=255)
                                                          gTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaaa                                            >  2:2884699/1‑101 (MQ=255)
                                                          gTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaaa                                            <  1:2548057/101‑1 (MQ=255)
                                                           tCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAg                                           <  1:2469584/101‑1 (MQ=255)
                                                            cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                                          >  1:2503225/1‑101 (MQ=255)
                                                            cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                                          >  1:3541182/1‑101 (MQ=255)
                                                            cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                                          >  1:3665/1‑101 (MQ=255)
                                                            cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                                          <  2:1742455/101‑1 (MQ=255)
                                                             aGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGgc                                         <  2:1644784/101‑1 (MQ=255)
                                                                gccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGcc                                      >  1:2794762/1‑101 (MQ=255)
                                                                 ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                     <  2:1913585/101‑1 (MQ=255)
                                                                 ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                     <  1:4503962/101‑1 (MQ=255)
                                                                 ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                     >  2:3892214/1‑101 (MQ=255)
                                                                  cgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGa                                    >  2:4188249/1‑101 (MQ=255)
                                                                    ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc                                  >  2:3225121/1‑101 (MQ=255)
                                                                    ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc                                  <  2:1384145/101‑1 (MQ=255)
                                                                    ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc                                  <  1:2152706/101‑1 (MQ=255)
                                                                    ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc                                  <  2:4586953/101‑1 (MQ=255)
                                                                     cAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt                                 <  1:74569/101‑1 (MQ=255)
                                                                      aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg                                >  1:3869257/1‑101 (MQ=255)
                                                                      aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATATg                                >  1:4705161/1‑101 (MQ=255)
                                                                       gTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGc                               >  1:1415264/1‑101 (MQ=255)
                                                                        ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt                              >  2:262629/1‑101 (MQ=255)
                                                                         tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc                             <  2:870231/101‑1 (MQ=255)
                                                                         tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc                             >  2:2119011/1‑101 (MQ=255)
                                                                         tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc                             >  1:4408880/1‑101 (MQ=255)
                                                                          cGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCg                            >  1:4580085/1‑101 (MQ=255)
                                                                            aTCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGc                          <  1:2014557/101‑1 (MQ=255)
                                                                              cG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCg                        <  1:2521087/101‑1 (MQ=255)
                                                                              cG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCg                        <  1:2928332/101‑1 (MQ=255)
                                                                                 tGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCgag                      <  1:641676/101‑1 (MQ=255)
                                                                                  gCTCAGCCCGGCCGTGTCCTCGGTCGAAAGCCCCGTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCgaga                     <  1:4078340/101‑1 (MQ=255)
                                                                                  gCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCgaga                     >  2:104523/1‑101 (MQ=255)
                                                                                    tCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAgc                   >  2:4115270/1‑101 (MQ=255)
                                                                                    tCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAgc                   >  2:1165860/1‑101 (MQ=255)
                                                                                    tCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAgc                   <  1:1000084/101‑1 (MQ=255)
                                                                                     cAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAgcg                  >  1:1080587/1‑101 (MQ=255)
                                                                                         ccGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGa              >  2:1449831/1‑101 (MQ=255)
                                                                                          cGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGaa             <  2:2310216/101‑1 (MQ=255)
                                                                                          cGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGaa             <  1:3820869/101‑1 (MQ=255)
                                                                                            gCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAg           <  2:3607042/101‑1 (MQ=255)
                                                                                            gCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAg           <  2:4254550/101‑1 (MQ=255)
                                                                                            gCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAg           >  1:918143/1‑101 (MQ=255)
                                                                                                tGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTgc       >  2:1949534/1‑101 (MQ=255)
                                                                                                tGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTgc       >  2:149266/1‑101 (MQ=255)
                                                                                                 gTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTgcg      <  2:4598147/101‑1 (MQ=255)
                                                                                                  tCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTgcgc     >  1:259481/1‑101 (MQ=255)
                                                                                                   ccTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGcc    >  1:1098221/1‑101 (MQ=255)
                                                                                                    cTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGCCg   <  2:3660/101‑1 (MQ=255)
                                                                                                     tCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGCCGa  <  2:1098135/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
CTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGCCGA  >  NC_014034/682494‑682694

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