Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,548 A→G 9.7% F687L (TTT→CTT)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5480A→G 9.7% 37.8 ‑1.5 ‑0.8185F687L (TTT→CTT) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (6/12):  ref base (86/81):  total (92/93)
Fisher's exact test strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.72e-01

Only 100 of 205 total aligned reads displayed.

TCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCG  >  NC_014034/681448‑681648
                                                                                                    |                                                                                                    
tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCTAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                                                                      <  1:3924149/101‑1 (MQ=255)
tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                                                                      <  1:1229753/101‑1 (MQ=255)
tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                                                                                      >  1:4064842/1‑101 (MQ=255)
 cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg                                                                                                     >  2:2554466/1‑101 (MQ=255)
  gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc                                                                                                    <  2:2770451/101‑1 (MQ=255)
   ggcgcCGCGACCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                                                                                                   <  2:215413/100‑1 (MQ=255)
    gcgcCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                                                                                                  <  1:586247/101‑2 (MQ=9)
     cgccgcAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                                                                                 <  2:4453878/101‑3 (MQ=9)
     cgccgcAACCTGCGTGCGCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                                                                                 <  2:867686/101‑3 (MQ=2)
      gccgcGATCTGCGTGCTGGTCCGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                                                                <  2:2013009/101‑1 (MQ=255)
      gccgcGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                                                                >  1:1388195/1‑101 (MQ=255)
      gccgcGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                                                                                <  1:4316790/101‑1 (MQ=255)
      gccgcAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                                                                <  2:3292091/101‑4 (MQ=9)
      gccgcAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                                                                <  1:667447/101‑4 (MQ=2)
      gccgcAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                                                                                <  1:3693229/101‑4 (MQ=2)
       ccgcAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcccttg                                                                                               <  1:1530769/101‑5 (MQ=2)
        cgcgATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa                                                                                              >  1:3146235/1‑101 (MQ=255)
          cgATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                                                                                            >  1:1190688/1‑101 (MQ=255)
          cgATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                                                                                            <  1:1596655/101‑1 (MQ=255)
           gATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                                                                                           >  1:814188/1‑101 (MQ=255)
            aTCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                                                                                          <  1:2143054/101‑1 (MQ=255)
            aTCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                                                                                          >  2:2214444/1‑101 (MQ=255)
             tCTGCTTGCTGTTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                                                                                         <  2:1426763/101‑1 (MQ=255)
             tCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc                                                                                         <  2:2681195/101‑1 (MQ=255)
              cTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg                                                                                        >  1:3959631/1‑101 (MQ=255)
                gCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc                                                                                      <  2:2353201/101‑1 (MQ=255)
                 cGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg                                                                                     >  2:418058/1‑101 (MQ=255)
                    gCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                                                                                  <  1:3611068/101‑1 (MQ=255)
                    gCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                                                                                  >  1:3409667/1‑101 (MQ=255)
                     cTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc                                                                                 >  1:2790705/1‑101 (MQ=255)
                      tGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc                                                                                <  1:281598/101‑1 (MQ=255)
                       ggTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                                                               >  2:1270008/1‑101 (MQ=255)
                       ggTCAGCCCCCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                                                                               <  2:4064556/101‑1 (MQ=255)
                        gTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg                                                                              <  1:4550589/101‑1 (MQ=255)
                              gcgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccgc                                                                        >  1:752865/1‑101 (MQ=255)
                                gcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                                                                      >  1:29892/1‑101 (MQ=255)
                                gcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                                                                      >  2:162564/1‑101 (MQ=255)
                                gcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                                                                      >  1:2191635/1‑101 (MQ=255)
                                gcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                                                                      >  1:1782283/1‑101 (MQ=255)
                                      cTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCg                                                                <  1:2214606/101‑1 (MQ=255)
                                        gcgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGcc                                                              <  2:2669798/101‑1 (MQ=255)
                                        gcgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGcc                                                              <  2:1063557/101‑1 (MQ=255)
                                         cgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCg                                                             >  2:1166015/1‑101 (MQ=255)
                                         cgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACAATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCg                                                             <  2:4022997/101‑1 (MQ=255)
                                          gcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGa                                                            >  2:382028/1‑101 (MQ=255)
                                           ccGAGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa                                                           <  2:2790489/101‑1 (MQ=255)
                                           ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa                                                           <  1:1078106/101‑1 (MQ=255)
                                           ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa                                                           <  1:4393322/101‑1 (MQ=255)
                                           ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAACGTCGCCGaa                                                           >  2:3383419/1‑101 (MQ=255)
                                            cGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCTCGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAg                                                          <  1:3252850/101‑1 (MQ=255)
                                             gAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGc                                                         >  2:1572887/1‑101 (MQ=255)
                                              ggCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                        <  1:345067/100‑1 (MQ=255)
                                              aGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGcc                                                        >  2:2520217/1‑101 (MQ=255)
                                               gCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCa                                                       >  1:4521414/1‑101 (MQ=255)
                                                cccAGCCCCCCCCCCCGCGCCGTCGTCAGCCCCCCGATGTCATCGGTGCCAAACCCGCGGACCCCCCCCGTGCCCAGCCCCCCCCCCGTCGCCGAAGCCaa                                                      <  2:3957075/101‑2 (MQ=255)
                                                 ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                     <  1:263546/101‑1 (MQ=255)
                                                 ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                     <  1:401004/101‑1 (MQ=255)
                                                 ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                     >  1:1435068/1‑101 (MQ=255)
                                                 ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                     <  1:2337569/101‑1 (MQ=255)
                                                    gccgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc                                                  <  2:2623988/99‑1 (MQ=255)
                                                    gccgccCCCCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc                                                  <  2:1388100/99‑1 (MQ=255)
                                                     ccgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccg                                                 <  2:424728/100‑1 (MQ=255)
                                                     aCGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccg                                                 <  1:418095/101‑1 (MQ=255)
                                                      gccgccgcCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgc                                                >  2:4405911/3‑101 (MQ=21)
                                                      cgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccgc                                                <  2:438659/101‑1 (MQ=255)
                                                          ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                            <  1:2978751/101‑1 (MQ=255)
                                                           cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                           <  1:3660460/101‑1 (MQ=255)
                                                           cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCCCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                           <  1:3307417/101‑1 (MQ=255)
                                                            gccTGCGCCGGCGTTAGCGCCGCGATGTCGTCGGTGTGAAAGCCGCGGATCACCACCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                          <  2:2285918/101‑1 (MQ=37)
                                                             ccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTgc                                         >  2:596007/1‑101 (MQ=37)
                                                                gcgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGcc                                      >  2:427657/1‑101 (MQ=255)
                                                                gcgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGcc                                      <  2:3433686/101‑1 (MQ=255)
                                                                 cgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCg                                     >  1:1827140/1‑101 (MQ=255)
                                                                     gtcgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCt                                 <  2:4521119/101‑1 (MQ=255)
                                                                       cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                               <  1:3803105/101‑1 (MQ=255)
                                                                       cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCCCCGCCGTGCCCAGCGCCGCCGCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                               <  2:771592/101‑1 (MQ=255)
                                                                        gtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAg                              >  2:3611504/1‑101 (MQ=255)
                                                                          cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                            <  2:4199873/101‑1 (MQ=255)
                                                                          cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                            >  2:1023701/1‑101 (MQ=255)
                                                                           aGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                           >  1:1005434/1‑101 (MQ=255)
                                                                             cgcCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccg                         >  2:3861546/1‑101 (MQ=255)
                                                                             ccccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccg                         <  1:883841/99‑1 (MQ=255)
                                                                               ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                       <  2:3226469/101‑1 (MQ=255)
                                                                               ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                       >  2:3977187/1‑101 (MQ=255)
                                                                                 aCGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgccgc                     >  1:4598540/1‑101 (MQ=255)
                                                                                  cGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCgccgccgcc                    <  2:458640/101‑1 (MQ=255)
                                                                                  cGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCgccgccgcc                    >  2:3024875/1‑101 (MQ=255)
                                                                                  cGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCgccgccgcc                    >  2:2816404/1‑101 (MQ=255)
                                                                                   gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                   <  1:3492451/101‑1 (MQ=255)
                                                                                   gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                   <  1:2670170/101‑1 (MQ=255)
                                                                                    aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                  >  1:4081772/1‑101 (MQ=255)
                                                                                    aTGTCATCGGTGCAAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                  >  1:2891011/1‑101 (MQ=255)
                                                                                        cATCGGTGCAAAACCGGCGGGTCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAcc              <  2:1769797/101‑1 (MQ=37)
                                                                                         ttCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg             <  2:2890795/100‑1 (MQ=255)
                                                                                           cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg           <  2:2072426/101‑1 (MQ=255)
                                                                                            ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc          >  2:3649712/1‑101 (MQ=255)
                                                                                                ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTCCACCGAGCCAAg      <  1:2520391/101‑1 (MQ=255)
                                                                                                  aaaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGcc    >  1:3377823/1‑101 (MQ=255)
                                                                                                    aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg  >  2:1449248/1‑101 (MQ=255)
                                                                                                    aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg  <  1:3852501/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
TCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCG  >  NC_014034/681448‑681648

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