Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,501 A→C 26.1% R702R (CGT→CGG RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5010A→C26.1% 160.1 ‑3.3 ‑1.2226R702R (CGT→CGGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (22/38):  ref base (88/78):  total (110/116)
Fisher's exact test strand distribution p-value = 3.51e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.84e-02

Only 100 of 238 total aligned reads displayed.

CAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCG  >  NC_014034/681402‑681601
                                                                                                    |                                                                                                    
cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAc                                                                                                     <  2:23854/101‑1 (MQ=255)
cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAc                                                                                                     >  2:1492412/1‑101 (MQ=255)
 agcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcg                                                                                                    <  2:3460454/101‑1 (MQ=255)
     cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgcc                                                                                                >  1:424765/1‑101 (MQ=255)
        agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCgccgcc                                                                                             >  1:2995339/1‑101 (MQ=255)
        agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCgccgcc                                                                                             >  2:28866/1‑101 (MQ=255)
            cgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTgcg                                                                                         <  2:2278728/101‑1 (MQ=255)
               cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                                      >  1:109252/1‑101 (MQ=255)
                 tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                                    <  2:2791967/101‑1 (MQ=255)
                    cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                                                 >  1:3846367/1‑101 (MQ=255)
                     gcTGATCTGCGCCAGCTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                                                <  1:4392501/101‑1 (MQ=255)
                     gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                                                >  2:177063/1‑101 (MQ=255)
                      gtGACCTGCGCCAGCTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAg                                                                               <  2:4567994/100‑1 (MQ=255)
                       tGATCTGCGCCAGATCCGCGGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                                              >  2:2594700/1‑101 (MQ=255)
                       tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                                              >  1:4382044/1‑101 (MQ=255)
                        gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAggg                                                                             >  1:2286080/1‑99 (MQ=255)
                         agccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcgc                                                                            <  2:501424/94‑1 (MQ=2)
                         aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGAGCCGAGCCCAGACGCGCCGGCTGCCACGTCGTCAgcgc                                                                            >  2:2780064/1‑101 (MQ=25)
                          gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                                                                           <  1:3896678/95‑1 (MQ=2)
                           cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                          <  2:3555257/101‑1 (MQ=255)
                           cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                          <  1:2866079/101‑1 (MQ=255)
                            tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                         <  2:4500227/101‑1 (MQ=255)
                            cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                         <  1:1572588/97‑1 (MQ=2)
                             gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                                        <  1:513803/98‑1 (MQ=2)
                               cccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                                      <  1:870736/100‑1 (MQ=255)
                                 cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                                    >  1:1467419/1‑101 (MQ=255)
                                  aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                                                   <  1:4219847/101‑1 (MQ=255)
                                   gacccgccgCCCCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                  <  1:1641566/98‑1 (MQ=255)
                                    aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGACAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTcct                                                                 >  2:2318498/1‑99 (MQ=255)
                                     tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                                >  1:2728406/1‑101 (MQ=255)
                                         cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCCCGATGTCATCGGTg                                                            >  2:2957314/9‑101 (MQ=9)
                                         cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                                                            <  2:1107852/93‑1 (MQ=9)
                                          ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           >  2:2166681/8‑101 (MQ=9)
                                          ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           >  2:3986805/8‑101 (MQ=9)
                                          cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           >  1:2114349/1‑101 (MQ=255)
                                          cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           <  2:336626/101‑1 (MQ=255)
                                           gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                                          <  1:262082/95‑1 (MQ=9)
                                           gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                                          >  2:773499/7‑101 (MQ=9)
                                             ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCTCCACGATGTCATCGGTGccga                                                        <  2:219624/97‑3 (MQ=9)
                                             cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                                                        >  1:2082330/1‑101 (MQ=255)
                                              tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                                                       <  1:2825950/98‑4 (MQ=9)
                                              tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                       >  2:4027787/1‑101 (MQ=255)
                                               ccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                                      <  2:332547/99‑5 (MQ=9)
                                               cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                                      >  2:2554466/1‑101 (MQ=255)
                                                caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                                                     >  1:1743236/2‑96 (MQ=9)
                                                 aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                                                    <  2:3846082/101‑1 (MQ=255)
                                                 aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                                    >  2:270665/1‑101 (MQ=9)
                                                 aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                                    >  1:1259292/1‑101 (MQ=2)
                                                  gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                                                   <  2:4083528/101‑2 (MQ=9)
                                                  gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                                                   <  1:586247/101‑2 (MQ=9)
                                                   cgccgcAACCTGCGTGCGCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                                  <  2:867686/101‑3 (MQ=2)
                                                    gccgcGATCTGCGTGCTGGTCCGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                 <  2:2013009/101‑1 (MQ=255)
                                                    gccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                 >  1:1388195/1‑101 (MQ=255)
                                                    gccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                 >  2:2669985/1‑101 (MQ=255)
                                                    gccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                 >  2:2504900/1‑101 (MQ=255)
                                                    gccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                 >  2:4420740/1‑101 (MQ=255)
                                                    gccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                 <  1:4316790/101‑1 (MQ=255)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                 <  1:3693229/101‑4 (MQ=2)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                 <  2:3292091/101‑4 (MQ=9)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                 <  1:667447/101‑4 (MQ=2)
                                                    gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                 <  2:3957563/101‑4 (MQ=9)
                                                    cccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                 <  2:2828514/100‑4 (MQ=2)
                                                      cgcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa                                               <  1:1554733/101‑1 (MQ=255)
                                                        cgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                                             <  1:3662809/101‑1 (MQ=255)
                                                         gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                            <  2:3503218/101‑1 (MQ=255)
                                                         gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                            >  1:814188/1‑101 (MQ=255)
                                                          aTCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                                           >  1:4453877/1‑101 (MQ=255)
                                                          aTCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                                           <  1:2143054/101‑1 (MQ=255)
                                                            cTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                                         <  1:1115076/101‑1 (MQ=255)
                                                                gTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                                     <  1:4079524/101‑1 (MQ=255)
                                                                  gCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                                   >  1:3409667/1‑101 (MQ=255)
                                                                   cTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc                                  >  1:2790705/1‑101 (MQ=255)
                                                                    tGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc                                 >  1:914778/1‑101 (MQ=255)
                                                                     ggTCGGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                <  1:243664/101‑1 (MQ=255)
                                                                     ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                >  2:1270008/1‑101 (MQ=255)
                                                                     ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                <  2:2995110/101‑1 (MQ=255)
                                                                           cgcgcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATTGCCGCCGTGCCCAGcgccg                          >  1:2643352/1‑101 (MQ=255)
                                                                           cccgcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccg                          <  2:1032194/99‑1 (MQ=255)
                                                                              gcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                       >  1:2191635/1‑101 (MQ=255)
                                                                              gcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                       >  2:1409090/1‑101 (MQ=255)
                                                                              gcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                       >  1:1782283/1‑101 (MQ=255)
                                                                              gcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                       >  1:4553438/1‑101 (MQ=255)
                                                                                 g.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg                    >  1:3936529/1‑101 (MQ=255)
                                                                                   acacgggctgAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGt                   >  2:3768488/10‑101 (MQ=255)
                                                                                     cggcgcGGAGCCCAGCCGCCCGTCCTGGGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCTCCGCCGTGCCCAGCGCCGCCACCGTCg                 <  2:4254640/99‑1 (MQ=25)
                                                                                     cTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCg                 <  1:2214606/101‑1 (MQ=255)
                                                                                      tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGc                >  1:4193086/1‑101 (MQ=255)
                                                                                       gcgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGcc               <  1:2903853/101‑1 (MQ=255)
                                                                                       gcgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGcc               <  2:3205687/101‑1 (MQ=255)
                                                                                        cgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCg              >  2:1166015/1‑101 (MQ=255)
                                                                                         gcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGa             >  2:418317/1‑101 (MQ=255)
                                                                                          ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa            <  1:1078106/101‑1 (MQ=255)
                                                                                           cGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCTCGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAg           <  1:3252850/101‑1 (MQ=255)
                                                                                            gAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGc          >  2:1572887/1‑101 (MQ=255)
                                                                                             aGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGTCACCGTCGCCGAAGcc         <  2:1190596/101‑1 (MQ=255)
                                                                                              gCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCa        >  1:4521414/1‑101 (MQ=255)
                                                                                                ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc      <  1:2337569/101‑1 (MQ=255)
                                                                                                ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc      <  1:401004/101‑1 (MQ=255)
                                                                                                   gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc   >  2:4383874/1‑101 (MQ=255)
                                                                                                   gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc   >  1:3471400/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCG  >  NC_014034/681402‑681601

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