Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,485 T→G 19.9% M708L (ATC→CTC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4850T→G19.9% 116.3 ‑2.8 ‑0.0210M708L (ATC→CTC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (15/27):  ref base (88/80):  total (103/107)
Fisher's exact test strand distribution p-value = 5.92e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.18e-01

Only 100 of 225 total aligned reads displayed.

GCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTC  >  NC_014034/681385‑681585
                                                                                                     |                                                                                                     
gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.At                                                                                                       <  2:4790044/101‑1 (MQ=255)
 ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATc                                                                                                      >  1:3072881/1‑101 (MQ=255)
 ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATc                                                                                                      <  1:1387547/101‑1 (MQ=255)
  cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                                                     >  1:3875933/1‑101 (MQ=255)
   aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGCGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                                                    >  1:4751252/1‑101 (MQ=255)
   aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                                                    >  1:3997497/1‑101 (MQ=255)
    tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgc                                                                                                   <  2:2505273/101‑1 (MQ=255)
    tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgc                                                                                                   >  2:2026532/1‑101 (MQ=255)
        ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCg                                                                                               >  1:4486761/1‑101 (MQ=255)
         gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGa                                                                                              >  1:3513976/1‑101 (MQ=255)
           cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                                            >  1:573455/1‑101 (MQ=255)
           cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                                            <  2:53838/101‑1 (MQ=255)
           cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                                            >  2:866258/1‑101 (MQ=255)
              gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCa                                                                                         <  2:3982639/101‑1 (MQ=255)
               cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAg                                                                                        <  1:519500/101‑1 (MQ=255)
               cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAg                                                                                        >  1:3356525/1‑101 (MQ=255)
                ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGa                                                                                       >  1:1621578/1‑101 (MQ=255)
                    cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcgc.g                                                                                  >  1:4137795/1‑101 (MQ=255)
                     ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACgc.gc                                                                                 >  2:1169326/1‑101 (MQ=255)
                      cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGc.gcc                                                                                <  1:4912142/101‑1 (MQ=255)
                       ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGc.gccg                                                                               >  1:127957/1‑101 (MQ=255)
                         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.gccgcc                                                                             >  2:3818583/1‑101 (MQ=255)
                         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.gccgcc                                                                             >  2:4523938/1‑101 (MQ=255)
                         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.gccgcc                                                                             <  1:1019758/101‑1 (MQ=255)
                              gcCGTCGCGCTGATCTGCGCCAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTgcgc                                                                        >  1:4490819/1‑101 (MQ=255)
                                cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcg                                                                      <  2:1850689/101‑1 (MQ=255)
                                 gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCATACGCcgt                                                                     >  2:4007452/1‑101 (MQ=255)
                                  tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtc                                                                    <  1:2542191/101‑1 (MQ=255)
                                   cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcg                                                                   >  2:2406248/1‑101 (MQ=255)
                                    gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcgt                                                                  <  2:2917503/101‑1 (MQ=255)
                                     cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCgtcgtc                                                                 >  1:4803815/1‑101 (MQ=255)
                                      gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCa                                                                <  1:3279146/101‑1 (MQ=255)
                                       cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAg                                                               >  2:3578378/1‑101 (MQ=255)
                                       cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGACTGCGCCGTCGCCAg                                                               >  1:4445062/1‑101 (MQ=255)
                                        tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GGCGCCTGCGCCGTCGTCAgc                                                              >  1:2964804/1‑101 (MQ=255)
                                        ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAgc                                                              <  1:1093843/92‑1 (MQ=2)
                                             cgcgcCAGATCCCCCGCCCCGAGCGCCGCCATCTTCCTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCAc                                                         <  2:5045812/100‑1 (MQ=25)
                                             cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCAc                                                         <  1:901341/97‑1 (MQ=2)
                                             cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCAc                                                         <  2:4606429/97‑1 (MQ=2)
                                              gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACg                                                        >  2:3791042/1‑101 (MQ=255)
                                              gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCTTGCTCGACAGGGCCAGA.AGCTGCGTCGAGCCCAG.CGCGGCCGCCTGCGCCGTCGTCAGCCCCACg                                                        >  2:1636690/4‑101 (MQ=2)
                                              gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg                                                        <  2:2598648/98‑1 (MQ=2)
                                               caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGa                                                       <  2:1836725/99‑1 (MQ=9)
                                               caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGACGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGGGCCACGa                                                       >  2:5021035/3‑101 (MQ=9)
                                                gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                      <  1:18655/101‑1 (MQ=255)
                                                accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                      <  1:3689095/100‑1 (MQ=2)
                                                 ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                     >  1:3937935/1‑101 (MQ=255)
                                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                 >  1:3870981/1‑101 (MQ=255)
                                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAAGC.GCAGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                 >  2:2497699/1‑101 (MQ=255)
                                                      tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                <  1:4607355/101‑1 (MQ=255)
                                                      tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGATCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                >  2:2825426/1‑101 (MQ=255)
                                                        cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                              >  2:599939/1‑101 (MQ=255)
                                                            gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                          <  2:2089590/95‑1 (MQ=12)
                                                            gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                          <  2:4154348/101‑1 (MQ=255)
                                                            gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                          >  2:173597/1‑101 (MQ=255)
                                                             ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATCTCATCGGTGCCa                                         <  1:2742357/101‑1 (MQ=255)
                                                             ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATCTCATCGGTGCCa                                         <  1:2732474/101‑1 (MQ=255)
                                                              ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCCCGATGTCATCGGTGccga                                        >  2:4967358/5‑99 (MQ=9)
                                                              ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                                        >  2:932608/5‑99 (MQ=9)
                                                              ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                                        <  2:1078842/97‑3 (MQ=9)
                                                              cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                                        >  1:2907733/1‑101 (MQ=255)
                                                               tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                                       <  1:3852182/98‑4 (MQ=9)
                                                                cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                      >  2:4819878/1‑101 (MQ=255)
                                                                cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                      >  1:852211/1‑101 (MQ=255)
                                                                 caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                     <  2:1892595/100‑6 (MQ=9)
                                                                  aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGCCATCGGTGCCAAAGcc                                    >  1:1315237/1‑101 (MQ=255)
                                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                    <  1:201605/101‑1 (MQ=9)
                                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                    <  2:4721859/101‑1 (MQ=9)
                                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                    <  2:3834733/101‑1 (MQ=9)
                                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                    >  1:2113145/1‑101 (MQ=2)
                                                                    cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                  >  1:4160371/1‑99 (MQ=2)
                                                                     gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                 <  2:1268038/101‑4 (MQ=9)
                                                                      ccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                                <  2:3644085/101‑5 (MQ=2)
                                                                      ccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGGCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                                <  2:394344/101‑5 (MQ=9)
                                                                       cgcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa                               <  1:229757/101‑1 (MQ=255)
                                                                        gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                              <  1:3971375/101‑1 (MQ=255)
                                                                        gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                              <  2:1467376/101‑1 (MQ=255)
                                                                         cgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                             >  2:1597303/1‑101 (MQ=255)
                                                                          gATCTGCGTGCTGGTCGGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                            <  2:794012/101‑1 (MQ=255)
                                                                          gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                            <  2:3072643/101‑1 (MQ=255)
                                                                          gATCGGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                            <  2:3513706/101‑1 (MQ=255)
                                                                             cTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                         >  2:762985/1‑101 (MQ=255)
                                                                             cTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                         <  2:4872820/101‑1 (MQ=255)
                                                                              tGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccgc                        <  2:3677560/101‑1 (MQ=255)
                                                                               gCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                       >  2:3277184/1‑101 (MQ=255)
                                                                                 gTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCCTCGGTGCAAAAGCCGCGGATCGCCGCCGt                     >  2:2361677/1‑101 (MQ=255)
                                                                                   gCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                   >  1:1947071/1‑101 (MQ=255)
                                                                                   gCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                   >  1:2064511/1‑101 (MQ=255)
                                                                                     tGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc                 >  1:1529280/1‑101 (MQ=255)
                                                                                     tGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc                 <  1:4384760/101‑1 (MQ=255)
                                                                                      tgTCAGCGCGCCG.ATCTGCGCCGACCCCAGACGC.GCCGCCTGGGCCGTCGGCAGCGCCACGATGTCATCGGGGCCAAAGCCGCGGATCGCCGCCGTGCCCa                <  2:3304852/100‑1 (MQ=37)
                                                                                      ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                <  2:4137478/101‑1 (MQ=255)
                                                                                        tCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc              <  1:1918060/101‑1 (MQ=255)
                                                                                        tCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc              >  1:2828010/1‑101 (MQ=255)
                                                                                         cAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg             <  2:4422640/101‑1 (MQ=255)
                                                                                          cgcgcgcCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcgc            <  1:2362161/100‑1 (MQ=255)
                                                                                           gcgcgcCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgcc           <  1:785559/101‑1 (MQ=255)
                                                                                               gcCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa       >  2:3324296/1‑101 (MQ=255)
                                                                                                ccG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAc      >  1:365400/1‑101 (MQ=255)
                                                                                                    aTCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGt   >  1:1812013/1‑101 (MQ=255)
                                                                                                     |                                                                                                     
GCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTC  >  NC_014034/681385‑681585

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