Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,561 C→A 4.3% A682A (GCG→GCT)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5610C→A 4.3% 10.2 ‑0.8 ‑0.2162A682A (GCG→GCT) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (5/2):  ref base (84/71):  total (89/73)
Fisher's exact test strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.29e-01

Only 100 of 179 total aligned reads displayed.

TCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAG  >  NC_014034/681461‑681659
                                                                                                    |                                                                                                  
 cTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                                                                                                   >  2:908970/1‑101 (MQ=255)
 cTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                                                                                                   >  1:2466565/1‑101 (MQ=255)
     gTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                                                                                               >  2:2308907/1‑101 (MQ=255)
       gCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                                                                                             <  1:4405694/101‑1 (MQ=255)
          ggTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                                                                          <  2:3574863/101‑1 (MQ=255)
            tCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc                                                                                        <  1:2101243/101‑1 (MQ=255)
            tCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc                                                                                        <  1:4161094/101‑1 (MQ=255)
              aGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcgc                                                                                      >  2:2347894/1‑101 (MQ=255)
                cgcgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccg                                                                                    <  2:817906/101‑1 (MQ=255)
                 gcgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccgc                                                                                   <  2:546586/101‑1 (MQ=255)
                  cgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCgccccc                                                                                  >  2:636821/1‑101 (MQ=255)
                   gcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                                                                                 <  2:1567954/101‑1 (MQ=255)
                    ccGATCTGCCCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAc                                                                                <  1:826334/101‑1 (MQ=255)
                        tCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGNCCAGCGCCGCCACCGTc                                                                            >  2:416275/1‑101 (MQ=255)
                         cTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGGCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCg                                                                           <  2:1396715/101‑1 (MQ=255)
                          tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCCGCCCCGCCCCCGTCGc                                                                          >  1:1046853/1‑101 (MQ=255)
                          tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGc                                                                          <  2:3125059/101‑1 (MQ=255)
                          tGCGCCGAGCCCAGACGCGCCGCCTGCCCCGTCGTCAGCTCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGc                                                                          <  2:3688540/101‑1 (MQ=255)
                            cgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCg                                                                        >  2:4029198/1‑101 (MQ=255)
                            cgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCg                                                                        >  2:1087591/1‑101 (MQ=255)
                            cgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCg                                                                        >  1:71686/1‑101 (MQ=255)
                             gcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGa                                                                       <  2:3497503/101‑1 (MQ=255)
                             gcCGAGCCCAGACGCCCCGCCTGCGCCGTCGTCAGCGCCACGACGTCATCGGTGCCAAAGCCGCGGATCGCCGCCCTGCCCAGCGCCGCCACCGTCGCCGa                                                                       >  2:315674/1‑101 (MQ=38)
                              ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa                                                                      <  1:3567722/101‑1 (MQ=255)
                               cGAGCCCAGGCGCGCCGCCTGCGCGGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAg                                                                     <  2:941831/101‑1 (MQ=255)
                               cGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAg                                                                     >  2:403839/1‑101 (MQ=255)
                               cGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCCCCGCAGTGCCCAGCGCCGCCACCGCCGCCgcac                                                                     >  1:1727166/1‑98 (MQ=37)
                                gAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGc                                                                    >  1:3406452/1‑101 (MQ=255)
                                 aGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                                   <  2:1664634/101‑1 (MQ=255)
                                 aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCCCCGTCGCCGAAGca                                                                   >  1:3551043/1‑100 (MQ=255)
                                 aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                                   >  2:2144609/1‑101 (MQ=255)
                                 aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                                   <  1:4058288/101‑1 (MQ=255)
                                    ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                                >  1:2971514/1‑101 (MQ=255)
                                     cAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgcg                                                               >  1:4212160/1‑101 (MQ=255)
                                       gACGCGCCCCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc                                                             <  1:1077993/101‑1 (MQ=255)
                                        aCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGGGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccg                                                            >  2:716535/1‑101 (MQ=255)
                                          ccgccgcCTGCGCCGTCGTCAGCGCCCCGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgca                                                          >  2:2096282/2‑100 (MQ=25)
                                            gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgcaag                                                        >  2:1748395/1‑98 (MQ=25)
                                            gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                        >  2:227697/1‑101 (MQ=255)
                                            gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                        >  2:558931/1‑101 (MQ=255)
                                            gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                        <  1:1198918/101‑1 (MQ=255)
                                             ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                                       >  1:3942465/1‑101 (MQ=255)
                                             ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                                       <  1:3820060/101‑1 (MQ=255)
                                             ccgccTGCCCCGTCGTCAGCGCCACGATGGCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                                       <  2:2054681/101‑1 (MQ=255)
                                              cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                                      <  1:3251154/101‑1 (MQ=255)
                                              cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                                      <  1:276880/101‑1 (MQ=255)
                                               gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                                     >  2:1219999/1‑101 (MQ=255)
                                               gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                                     >  1:3640733/1‑101 (MQ=255)
                                                 cTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgcg                                                   <  2:2873621/101‑1 (MQ=255)
                                                 cTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgcg                                                   <  2:1736495/101‑1 (MQ=255)
                                                  tGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgcgc                                                  <  2:1595102/101‑1 (MQ=255)
                                                   gcgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGGCGCCGAAGCCAGCGCCGCGATCTGCGcc                                                 <  2:686285/101‑1 (MQ=255)
                                                   gcgcCGTCGTCAGCGCCACGATGTCATCGGTCCCCAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAACCAGCGCCGCGATCTTCGcc                                                 >  1:2519430/1‑101 (MQ=37)
                                                   gatcggtggtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGcc                                                 <  2:4373723/93‑1 (MQ=255)
                                                    cgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCTCCg                                                >  1:1178403/1‑101 (MQ=255)
                                                      ccGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTg                                              >  2:3880636/1‑101 (MQ=255)
                                                      ccGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTg                                              <  1:1867149/101‑1 (MQ=255)
                                                       cgtcgtCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGc                                             >  2:486234/1‑101 (MQ=255)
                                                       cgtcgtCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGc                                             >  1:3031839/1‑101 (MQ=255)
                                                          cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTcc                                          >  2:691028/1‑100 (MQ=255)
                                                          cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                          >  1:2009159/1‑101 (MQ=255)
                                                          cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                          >  1:3779551/1‑101 (MQ=255)
                                                             cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCaaga                                       >  1:2698278/1‑98 (MQ=255)
                                                             cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                       >  2:2131196/1‑101 (MQ=255)
                                                             cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                       >  1:2690516/1‑101 (MQ=255)
                                                              aGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                      <  2:3332523/101‑1 (MQ=255)
                                                              aGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                      >  1:3332760/1‑101 (MQ=255)
                                                               gcgcCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgcc                                     <  1:1318374/101‑1 (MQ=255)
                                                                  ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                  <  1:4007466/101‑1 (MQ=255)
                                                                  ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                  >  1:1531272/1‑101 (MQ=255)
                                                                   cACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgccg                                 >  1:2216190/1‑101 (MQ=255)
                                                                      gATGTCGTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                              <  2:4190883/101‑1 (MQ=255)
                                                                      gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                              >  2:963163/1‑101 (MQ=255)
                                                                      gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                              >  1:576900/1‑101 (MQ=255)
                                                                      gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                              <  1:2916727/101‑1 (MQ=255)
                                                                      gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                              >  2:2019109/1‑101 (MQ=255)
                                                                      gATGTCATCGGTGCCAAAGCCGCGGAGCCCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTTCGCCGTGCTCAGCGCCGCCGCCt                              <  2:2519267/101‑1 (MQ=255)
                                                                       aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                             >  1:2600220/1‑101 (MQ=255)
                                                                       aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                             >  2:2582160/1‑101 (MQ=255)
                                                                       aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                             <  2:1226531/101‑1 (MQ=255)
                                                                         gTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCa                           <  1:1784879/101‑1 (MQ=255)
                                                                           cATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAcc                         >  2:1966856/1‑101 (MQ=255)
                                                                            aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                        <  2:2422656/101‑1 (MQ=255)
                                                                              cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                      >  2:2520539/1‑101 (MQ=255)
                                                                              cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                      >  2:4356492/1‑101 (MQ=255)
                                                                               ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                     >  2:3590176/1‑101 (MQ=255)
                                                                               ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                     <  2:3406215/101‑1 (MQ=255)
                                                                               ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                     >  2:776302/1‑101 (MQ=255)
                                                                               ggTGCCAAAGCCGCGGACCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCTGTGCTCAGCGCCGCCGCCTGCACCGAGc                     <  2:44738/101‑1 (MQ=255)
                                                                                  gCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCCCCGAGCCaa                  >  1:1978056/1‑101 (MQ=255)
                                                                                   ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                 <  1:2309057/101‑1 (MQ=255)
                                                                                   ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                 >  2:27554/1‑101 (MQ=255)
                                                                                   ccAAAGCCGCGGATCGCCGCCGTCCCCAGCCCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                 <  1:1293138/101‑1 (MQ=255)
                                                                                   ccAAAGCCGCGGATCCCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                 <  2:4211899/101‑1 (MQ=255)
                                                                                     aaaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGcc               >  1:3654916/1‑101 (MQ=255)
                                                                                       aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg             <  1:416315/101‑1 (MQ=255)
                                                                                        gCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCgg            >  1:3543199/1‑101 (MQ=255)
                                                                                          cgcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGcc          >  2:2216707/1‑101 (MQ=255)
                                                                                                tCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCa    <  1:4379167/101‑1 (MQ=255)
                                                                                                  cccgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg  <  1:227724/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
TCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAG  >  NC_014034/681461‑681659

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG