Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,501 A→C 18.5% 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→C18.5% 189.8 ‑3.9 ‑1.4403R702R (CGT→CGGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (23/52):  ref base (150/178):  total (173/230)
Fisher's exact test strand distribution p-value = 1.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.88e-02

Only 100 of 430 total aligned reads displayed.

CCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCG  >  NC_014034/681401‑681601
                                                                                                     |                                                                                                    
ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGa                                                                                                      <  1:7957038/101‑1 (MQ=255)
  agccaGGGGAGCGCCGTGGCGCTGATCTGCGCCCGATCCGCCGCTTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcg                                                                                                    <  2:6989282/101‑1 (MQ=37)
      cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGTCCAGACGcgcc                                                                                                <  1:5735789/101‑1 (MQ=255)
      cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgcc                                                                                                >  1:4181556/1‑101 (MQ=255)
       ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgccg                                                                                               <  2:7650308/101‑1 (MQ=255)
        gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgccgc                                                                                              >  1:433832/1‑101 (MQ=255)
         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGGgccccc                                                                                             >  1:7022306/1‑101 (MQ=255)
         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCgccgcc                                                                                             <  2:6083718/101‑1 (MQ=255)
          gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                                                            >  1:8500082/1‑101 (MQ=255)
          gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                                                            >  2:8078913/1‑101 (MQ=255)
          gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                                                            <  2:7119944/101‑1 (MQ=255)
            gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTgc                                                                                          <  2:5484785/101‑1 (MQ=255)
               ccGTCGCGCTTATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCGGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                                                                                       <  2:2078072/101‑1 (MQ=255)
                cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                                      >  1:1254342/1‑101 (MQ=255)
                cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                                                                                      <  2:1879837/101‑1 (MQ=255)
                  tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                                    >  1:2967200/1‑101 (MQ=255)
                  tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                                    <  1:387027/101‑1 (MQ=255)
                    gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcgt                                                                                  <  2:8547386/101‑1 (MQ=255)
                    gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcgc                                                                                  >  2:9027387/1‑100 (MQ=255)
                      gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCGAGACGCGCCGCCTGCGCCGTCGTCa                                                                                >  1:3781378/1‑101 (MQ=255)
                      gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                                                >  1:1413430/1‑101 (MQ=255)
                       cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAg                                                                               >  1:4902202/1‑101 (MQ=255)
                        tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcggggcattt                                                                              <  2:3788813/101‑9 (MQ=255)
                        ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                                              <  1:5501695/92‑1 (MQ=2)
                         gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                                             >  2:1026708/1‑101 (MQ=255)
                          agccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGAGCCGCCGCCTGCGCCGTCGTCAgcgc                                                                            >  1:8956321/8‑101 (MQ=9)
                          aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                                            >  2:6534962/1‑101 (MQ=255)
                           tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                                           <  1:204451/101‑1 (MQ=255)
                            ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                          <  2:4591692/96‑1 (MQ=2)
                            ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                                                                          <  2:5586329/96‑1 (MQ=2)
                             cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                         <  1:2408226/97‑1 (MQ=2)
                             cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                                         <  2:6207677/97‑1 (MQ=2)
                              gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCATACGCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                                        <  2:9175464/101‑1 (MQ=255)
                              gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                                        <  1:871275/101‑1 (MQ=255)
                              gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                                        <  1:5068900/98‑1 (MQ=2)
                                 ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg                                                                     <  2:7280903/101‑1 (MQ=255)
                                  cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                                    >  2:3976649/1‑101 (MQ=255)
                                    gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                  <  1:7369233/101‑1 (MQ=255)
                                    gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                  >  1:294042/1‑101 (MQ=255)
                                    gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                  >  2:8433566/1‑101 (MQ=255)
                                    gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                                  >  2:7908291/1‑101 (MQ=255)
                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                                 <  1:6381518/101‑1 (MQ=255)
                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                                 <  1:2453450/101‑1 (MQ=255)
                                      tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                                <  2:485663/101‑1 (MQ=255)
                                      tCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                                <  1:1011637/101‑1 (MQ=255)
                                          ccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                                                            >  2:2207637/1‑101 (MQ=255)
                                           ggtttgcaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           >  2:4529928/8‑101 (MQ=9)
                                           ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           <  1:3341139/94‑1 (MQ=9)
                                           ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                                           >  1:584162/8‑101 (MQ=9)
                                            gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                                          <  1:7099441/95‑1 (MQ=9)
                                            gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                                          <  2:139398/95‑1 (MQ=9)
                                            gccGCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                                          <  1:1074174/101‑1 (MQ=255)
                                             tttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCTTCGGTGccg                                                         <  1:8802082/96‑2 (MQ=9)
                                             ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                                                         <  1:5411617/101‑1 (MQ=255)
                                               tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                       >  1:3537636/1‑101 (MQ=255)
                                               tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCCCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                       >  2:3681176/1‑101 (MQ=255)
                                                ccaGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                                      >  1:1770656/3‑97 (MQ=2)
                                                 gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                                     >  1:8281636/1‑101 (MQ=255)
                                                  aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                                                    <  2:7750037/101‑1 (MQ=255)
                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                                    >  1:4450472/1‑101 (MQ=9)
                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                                    >  2:8913765/1‑101 (MQ=9)
                                                   gccgcgtgATCTGCGTGCCGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcg                                                   <  2:2923858/94‑1 (MQ=255)
                                                    cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                                  <  2:5074496/101‑3 (MQ=9)
                                                     gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                 >  1:6052557/1‑98 (MQ=2)
                                                      ccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                                                >  2:5708540/1‑97 (MQ=2)
                                                          gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                            <  2:2345395/101‑1 (MQ=255)
                                                          gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                            >  2:1229520/1‑101 (MQ=255)
                                                           aTCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                                           >  2:8336724/1‑101 (MQ=255)
                                                           aTCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                                           >  2:1696301/1‑101 (MQ=255)
                                                           aTCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGGGCCACAGCCGCGGCTccc                                           >  1:4285097/1‑99 (MQ=255)
                                                               gCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                                       <  1:7195830/101‑1 (MQ=255)
                                                               gCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                                       <  1:7716925/101‑1 (MQ=255)
                                                               gCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                                       <  1:1877876/101‑1 (MQ=255)
                                                                cGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg                                      <  2:5211326/101‑1 (MQ=255)
                                                                 gTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                                     >  2:7370179/1‑101 (MQ=255)
                                                                 gTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCAGCCTGAGCCGTCGTCAGCGCAACGATGTAATCGGTGCCAAAGACGCAGATAGCCGCCGt                                     >  1:9098958/1‑101 (MQ=25)
                                                                 gTGCGGGTCAGCCCCCCG.ATCGGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                                     <  2:2989397/101‑1 (MQ=37)
                                                                    cTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCCCCGCCTGCGCCGTCGTCATCGCCACGATATCATCGGTGCCAAACCCACCGAGcaccgcgccgcc                                  >  2:4257220/1‑95 (MQ=18)
                                                                      ggccAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                <  2:3574253/98‑1 (MQ=255)
                                                                      ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                <  1:6231834/101‑1 (MQ=255)
                                                                      ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                <  2:1986255/101‑1 (MQ=255)
                                                                       gTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg                               >  2:7877737/1‑101 (MQ=255)
                                                                         cAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg                             >  1:2527603/1‑101 (MQ=255)
                                                                         cAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg                             >  2:642622/1‑101 (MQ=255)
                                                                             gcgcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccgc                         <  1:7566144/101‑1 (MQ=255)
                                                                               gcCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                       <  1:7858681/101‑1 (MQ=255)
                                                                                 cgaaTCTGCGCCGGGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAc                      <  2:2684192/98‑1 (MQ=255)
                                                                                 cG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc                     <  2:2169105/101‑1 (MQ=255)
                                                                                 cG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc                     <  2:8187661/101‑1 (MQ=255)
                                                                                      cTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCg                 <  1:6718068/101‑1 (MQ=255)
                                                                                       tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGc                <  2:590886/101‑1 (MQ=255)
                                                                                          gcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGa             >  1:3151372/1‑101 (MQ=255)
                                                                                           ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa            <  1:524894/101‑1 (MQ=255)
                                                                                              aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc         >  1:2754991/1‑101 (MQ=255)
                                                                                               gCCCAGCCGCGCCCCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCa        <  1:3681285/101‑1 (MQ=255)
                                                                                               gCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCa        <  1:223469/101‑1 (MQ=255)
                                                                                                 ccAGGCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGCATCGCCGCCGTGCCCAGAGCCGCCACCGTCGCCGAAGCCAgc      <  2:4371235/101‑1 (MQ=255)
                                                                                                  cAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgcg     >  2:2601470/1‑101 (MQ=255)
                                                                                                    gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc   >  2:7470134/1‑101 (MQ=255)
                                                                                                    gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc   <  1:4837606/101‑1 (MQ=255)
                                                                                                     |                                                                                                    
CCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCG  >  NC_014034/681401‑681601

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