Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,658 A→G 12.6% L650P (CTT→CCA)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,6580A→G12.6% 52.3 ‑3.1 ‑0.0188L650P (CTT→CCT) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (15/9):  ref base (64/100):  total (79/109)
Fisher's exact test strand distribution p-value = 4.47e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.90e-01

Only 100 of 201 total aligned reads displayed.

CGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGT  >  NC_014034/681558‑681758
                                                                                                    |                                                                                                    
cgccgcCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCaa                                                                                                      <  2:3527060/101‑1 (MQ=255)
cgccgcCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCaa                                                                                                      >  1:2700687/1‑101 (MQ=255)
 gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                                                                                                     >  2:2773247/1‑101 (MQ=255)
 gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                                                                                                     <  2:220353/101‑1 (MQ=255)
   cgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGc                                                                                                   >  2:2697249/1‑101 (MQ=255)
   cgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGc                                                                                                   <  1:2900303/101‑1 (MQ=255)
    gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCACGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCt                                                                                                  <  2:3329643/101‑1 (MQ=255)
      cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                                                                                                >  1:1721406/1‑101 (MQ=255)
        tGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                                                                                              <  1:638798/101‑1 (MQ=255)
        tGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                                                                                              <  1:2441310/101‑1 (MQ=255)
         gCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgaga                                                                                             <  2:1367173/101‑1 (MQ=255)
         gCCCAGCGCCCCCACCGTGGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgaga                                                                                             <  2:1047108/101‑1 (MQ=255)
             aGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCGCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgc                                                                                         >  1:2947552/1‑101 (MQ=255)
              gcgcCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcg                                                                                        >  2:2910587/1‑101 (MQ=255)
               cgccgcCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcgc                                                                                       <  2:834271/101‑1 (MQ=255)
                   gccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAc                                                                                   >  1:1006707/1‑101 (MQ=255)
                    ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCCCCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc                                                                                  <  1:4205901/101‑1 (MQ=255)
                    ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCCCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGCcc                                                                                  >  1:747687/1‑101 (MQ=255)
                     cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt                                                                                 <  2:3675123/101‑1 (MQ=255)
                       ccGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCAGa                                                                               >  2:190420/1‑101 (MQ=255)
                        cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                              <  1:500548/101‑1 (MQ=255)
                        cGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGAt                                                                              <  1:1089326/101‑1 (MQ=255)
                             ccGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCCCAAAGCTGGAGAGCGAGCCGACCTGATCGGTg                                                                         >  1:705985/1‑101 (MQ=255)
                              cGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGc                                                                        <  1:2889695/101‑1 (MQ=255)
                               gtggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGcc                                                                       >  1:3398693/9‑100 (MQ=2)
                               gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCt                                                                       >  2:3212533/1‑101 (MQ=255)
                                tggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGccc                                                                      >  1:1375244/8‑99 (MQ=2)
                                tggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGccc                                                                      >  2:2847620/8‑99 (MQ=2)
                                aaGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc                                                                      >  2:1568506/1‑101 (MQ=255)
                                 ggcaatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCa                                                                     <  1:2510584/95‑1 (MQ=2)
                                 aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                     <  2:357191/101‑1 (MQ=255)
                                 aGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCa                                                                     <  1:3006872/101‑1 (MQ=255)
                                  gCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAg                                                                    >  1:1321171/1‑101 (MQ=255)
                                  gCCAGCGCCGCGATCTGCGCCGTGCGCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAg                                                                    <  2:1987696/101‑1 (MQ=255)
                                   caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                   >  1:1901590/5‑100 (MQ=2)
                                   caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                   >  1:858780/5‑100 (MQ=2)
                                   caatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAga                                                                   >  1:2582136/5‑100 (MQ=2)
                                    cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCt                                                                  <  2:3189532/101‑1 (MQ=255)
                                    aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                  >  1:3433027/4‑99 (MQ=2)
                                    aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                  >  2:4120183/4‑99 (MQ=2)
                                    aatcgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgac                                                                  >  1:4047628/4‑99 (MQ=2)
                                    aaGCGCAGAGATCTGAGCCGGGCAGCGCGCGGCCGCGTACACTAAACCAAGCCCGCCCACCGCAAGGCTGGAGAGCTCGCCGACCTGATCGGTGCTCAGCt                                                                  <  2:2956129/100‑1 (MQ=255)
                                       cgccgcGATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGACCAGCGCCGCGACCTGATCGGTGCCCAgactgg                                                               >  1:3610598/1‑96 (MQ=2)
                                       ccccgcGATCTGCGCCGCGCTCGGCGCCGCGATCGGCGCGGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactgg                                                               <  2:2267179/99‑6 (MQ=2)
                                        gccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGcc                                                              >  1:92791/1‑101 (MQ=255)
                                        gccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGcc                                                              >  2:2528702/1‑101 (MQ=255)
                                        gccgcGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGcc                                                              <  1:2827358/101‑1 (MQ=255)
                                          cgcgATCTGCGCCGCGCTCAGCGCCGCGATCTGCGCCGAGCCGAGCCCGCCGACCGCGGCACTGGCCAGCGCCGCGACCTGATCGGTGCCCAgactggcca                                                            <  2:3905978/101‑9 (MQ=2)
                                             gATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCt                                                         <  2:92767/101‑1 (MQ=255)
                                              aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                        <  2:4090517/101‑1 (MQ=255)
                                              aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                        >  1:2358092/1‑101 (MQ=255)
                                              aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                        >  2:1806225/1‑101 (MQ=255)
                                              aTCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGACCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTg                                                        <  1:2626221/101‑1 (MQ=255)
                                                  gcgcCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTc                                                    <  1:3212772/101‑1 (MQ=255)
                                                       gTGCTCAGCGCCGCCGCCTGCACCGAACCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCt                                               <  2:2090585/101‑1 (MQ=255)
                                                        tGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTc                                              <  2:495916/101‑1 (MQ=255)
                                                          cTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAg                                            <  1:182599/101‑1 (MQ=255)
                                                            cAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCCCGCTCAGGc                                          >  1:3894783/1‑101 (MQ=255)
                                                               cgccgccgcCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCg                                       <  1:3332894/101‑1 (MQ=255)
                                                                gccgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGa                                      <  2:2998060/101‑1 (MQ=255)
                                                                gccgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGa                                      <  2:2129847/101‑1 (MQ=255)
                                                                  cgccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGaaa                                    <  1:1414177/101‑1 (MQ=255)
                                                                   gccgccTCCACCGAGCCAAGCCCGCCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAc                                   <  2:552308/101‑1 (MQ=255)
                                                                    ccgccTGCCCCGAGCCAAGCCCGGCCACCGCAAGTCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                  <  1:2472735/101‑1 (MQ=255)
                                                                    ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                  <  2:1321071/101‑1 (MQ=255)
                                                                    ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                  <  1:2149336/101‑1 (MQ=255)
                                                                    ccgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                  >  2:2351609/1‑101 (MQ=255)
                                                                    ccgccGGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAcc                                  <  1:1258183/101‑1 (MQ=255)
                                                                     cgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccg                                 <  1:2773466/101‑1 (MQ=255)
                                                                     cgccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccg                                 >  2:735042/1‑101 (MQ=255)
                                                                      gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                >  1:2003087/1‑101 (MQ=255)
                                                                      gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                >  1:81650/1‑101 (MQ=255)
                                                                      gccTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgc                                >  2:1596598/1‑101 (MQ=255)
                                                                       ccTGCCCCGAGCCAATCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAAccgcc                               <  2:4069730/101‑1 (MQ=255)
                                                                         tGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACcgccgc                             <  2:2143761/101‑1 (MQ=255)
                                                                          gCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCcgcg                            >  1:4141520/1‑101 (MQ=255)
                                                                           cACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCgcgc                           >  2:468369/1‑101 (MQ=255)
                                                                           cACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTCATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCgcgc                           >  2:333747/1‑101 (MQ=255)
                                                                            cccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCt                          <  1:241396/100‑1 (MQ=255)
                                                                             ccGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTg                         <  1:3789870/101‑1 (MQ=255)
                                                                              cGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGc                        >  1:235959/1‑101 (MQ=255)
                                                                                 gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGcccc                     >  2:3062738/1‑100 (MQ=255)
                                                                                 gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                     <  2:2885003/101‑1 (MQ=255)
                                                                                 gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                     <  1:1568615/101‑1 (MQ=255)
                                                                                 gCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCa                     <  1:895307/101‑1 (MQ=255)
                                                                                   cAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGc                   >  1:2464731/1‑101 (MQ=255)
                                                                                   cAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCCCGCAGCCCCGc                   >  1:4153894/1‑101 (MQ=255)
                                                                                    aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcc                  <  2:2473775/101‑1 (MQ=255)
                                                                                    aaGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACAGCCGCGCTGCCCAGcc                  >  2:1066421/1‑101 (MQ=255)
                                                                                      gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                <  1:2200595/101‑1 (MQ=255)
                                                                                      gCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGcccc                <  2:2367566/101‑1 (MQ=255)
                                                                                            ccACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATc          >  2:3470938/1‑101 (MQ=255)
                                                                                            ccACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATc          <  2:264786/101‑1 (MQ=255)
                                                                                             cACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCt         <  2:2600474/101‑1 (MQ=255)
                                                                                               ccGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTgc       <  2:4174359/101‑1 (MQ=255)
                                                                                                  cAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGcc    >  1:760534/1‑101 (MQ=255)
                                                                                                   aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCg   >  2:1099317/1‑101 (MQ=255)
                                                                                                   aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCg   >  2:1401504/1‑101 (MQ=255)
                                                                                                   aaGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCg   >  1:2608452/1‑101 (MQ=255)
                                                                                                    aGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGt  >  2:2687837/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTCAGCTTGCCGATCTGCGTCGCGCTCAGGCCCGAAACCGCCGCGCTGCCCAGCCCCGCGATCTGCGCCGT  >  NC_014034/681558‑681758

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