Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,590 A→T 2.8% S673T (TCG→ATG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5900A→T 2.8% 2.6 ‑1.9 ‑0.0141S673T (TCG→ACG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (79/58):  total (83/58)
Fisher's exact test strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.74e-01

Only 100 of 155 total aligned reads displayed.

CGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTC  >  NC_014034/681492‑681690
                                                                                                  |                                                                                                    
cGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAg                                                                                                    >  1:4274446/1‑101 (MQ=255)
     ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                                                               >  2:1778851/1‑101 (MQ=255)
     ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                                                               >  1:458776/1‑101 (MQ=255)
          gccgccgcCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgc                                                                                          >  2:4623526/3‑101 (MQ=21)
          cgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCACCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccgc                                                                                          >  2:1109758/1‑101 (MQ=255)
           gcgcCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCcgcg                                                                                         <  1:1003607/101‑1 (MQ=255)
           ccgccgcCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgca                                                                                         >  1:762990/2‑100 (MQ=25)
            cgccgcCTGCGCCGTCGTCACCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCCTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCCCGa                                                                                        >  2:611011/1‑101 (MQ=255)
             gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGACGAAGCCAGCGCCGCGAt                                                                                       <  1:3477587/101‑1 (MQ=255)
              ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgcaagc                                                                                      >  1:2156272/1‑97 (MQ=25)
              ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                                                                      <  2:2451976/101‑1 (MQ=255)
                gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGACGCCAGCGCCGCGATCTg                                                                                    >  2:4266711/1‑101 (MQ=255)
                gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                                                                    >  1:4439962/1‑101 (MQ=255)
                  cTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgcg                                                                                  >  2:1574047/1‑101 (MQ=255)
                   tGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgcgc                                                                                 >  2:2932050/1‑101 (MQ=255)
                    gcgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGcc                                                                                <  2:3342928/101‑1 (MQ=255)
                     cgcTGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCg                                                                               >  2:4336254/1‑101 (MQ=255)
                       ccGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTg                                                                             <  1:1039683/101‑1 (MQ=255)
                        cgtcgtCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGc                                                                            >  2:685097/1‑101 (MQ=255)
                          gcgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTc                                                                          <  1:3100794/100‑1 (MQ=255)
                           cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                         >  1:1542141/1‑101 (MQ=255)
                           cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                         <  2:3340780/101‑1 (MQ=255)
                              cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                      >  2:1409236/1‑101 (MQ=255)
                              cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                      >  2:173651/1‑101 (MQ=255)
                              cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                      >  1:3229965/1‑101 (MQ=255)
                               aGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                                                     <  2:2616239/101‑1 (MQ=255)
                               aGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                                                     <  2:2815000/101‑1 (MQ=255)
                                  gcCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                  >  1:2268997/1‑101 (MQ=255)
                                   ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                                                 <  1:999382/101‑1 (MQ=255)
                                   ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                                                 <  2:4282560/101‑1 (MQ=255)
                                   ccACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgcc                                                                 >  2:2930200/1‑101 (MQ=255)
                                    cACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgccg                                                                >  2:1021113/1‑101 (MQ=255)
                                    cACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgccg                                                                >  2:3781638/1‑101 (MQ=255)
                                      ggATGTCATCGGTGCCAAAGCCGCGGACCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCgccgccgcc                                                              <  2:114324/100‑1 (MQ=255)
                                       gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                                             >  1:48430/1‑101 (MQ=255)
                                       gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                                             >  2:3962861/1‑101 (MQ=255)
                                       gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                                             >  2:1671833/1‑101 (MQ=255)
                                           tCATCGGTGCCAAGGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                                         <  2:3843179/101‑1 (MQ=255)
                                           tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                                         <  1:4266958/101‑1 (MQ=255)
                                           tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                                         >  2:4130361/1‑101 (MQ=255)
                                           tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                                         >  1:1796806/1‑101 (MQ=255)
                                           tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                                                         <  2:2628280/101‑1 (MQ=255)
                                            cATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAcc                                                        >  1:3005167/1‑101 (MQ=255)
                                             aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                                                       >  2:2070493/1‑101 (MQ=255)
                                             aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                                                       <  2:2381198/101‑1 (MQ=255)
                                             aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                                                       <  2:2980091/101‑1 (MQ=255)
                                             aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                                                       <  1:4529952/101‑1 (MQ=255)
                                              tCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGa                                                      <  1:4632524/101‑1 (MQ=255)
                                               cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                                                     >  1:1621610/1‑101 (MQ=255)
                                               cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                                                     >  2:3472764/1‑101 (MQ=255)
                                                ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                                    <  2:34450/101‑1 (MQ=255)
                                                ggTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                                                    >  2:2673371/1‑101 (MQ=255)
                                                 gTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGcc                                                   <  1:1778960/101‑1 (MQ=255)
                                                 gTGCCAAAGCCCCCGACCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCTTGCTCAGCGCCGCCGCCTGCACCGAGcc                                                   <  1:611075/101‑1 (MQ=38)
                                                   gCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCaa                                                 >  1:2800090/1‑101 (MQ=255)
                                                    ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                                                <  1:4664792/101‑1 (MQ=255)
                                                    ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                                                >  2:472406/1‑101 (MQ=255)
                                                    ccAAAGCCCCGGGCCGCCACCGTCCCCCGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                                                <  2:3118988/101‑1 (MQ=37)
                                                     cAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGc                                               >  1:1249426/1‑101 (MQ=255)
                                                      aaaGCCGCGGATCCCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGcc                                              <  1:1672675/101‑1 (MQ=255)
                                                       aaGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGccc                                             <  1:3195716/101‑1 (MQ=255)
                                                        aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg                                            >  1:3525011/1‑101 (MQ=255)
                                                        aGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCg                                            >  1:916148/1‑101 (MQ=255)
                                                         gCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCgg                                           >  2:4754802/1‑101 (MQ=255)
                                                          cgggggacccccgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGc                                          <  1:2673525/91‑1 (MQ=255)
                                                           cgcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGcc                                         <  1:1895396/101‑1 (MQ=255)
                                                           cccgGTTCCCCGCCGGCCCCAGCGCGGCCACCGGCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGcc                                         <  2:2986754/99‑1 (MQ=37)
                                                             cgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAc                                       <  1:2810972/101‑1 (MQ=255)
                                                               gATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCg                                     >  2:2220922/1‑101 (MQ=255)
                                                                 tCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAACCCCGCCCATCGCa                                   >  2:1216975/1‑101 (MQ=255)
                                                                  cgccgcCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCaa                                  >  1:1434534/1‑101 (MQ=255)
                                                                   gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                                 >  1:2451836/1‑101 (MQ=255)
                                                                   gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg                                 >  1:1307996/1‑101 (MQ=255)
                                                                    ccgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAgg                                <  1:3811402/101‑1 (MQ=255)
                                                                        cgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTgg                            <  1:4087081/101‑1 (MQ=255)
                                                                         gTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGga                           >  1:240117/1‑101 (MQ=255)
                                                                          tGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                          <  1:4128431/101‑1 (MQ=255)
                                                                          tGCCCAGCCCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCGGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGgag                          <  2:3789085/101‑1 (MQ=255)
                                                                            cccCGCCCCCCCACCGTCGCGGACGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGagag                        <  2:4600928/101‑1 (MQ=37)
                                                                            cccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGagag                        <  2:3812106/101‑1 (MQ=255)
                                                                            cccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGAACGCCTGGagag                        >  1:1591494/1‑101 (MQ=255)
                                                                             ccAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgc                       >  2:4050138/1‑101 (MQ=255)
                                                                              cAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcg                      >  1:1044029/1‑101 (MQ=255)
                                                                                gcgcCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcg                    <  2:1621510/101‑1 (MQ=255)
                                                                                 cgccgcCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAgcgcgc                   >  2:4148526/1‑101 (MQ=255)
                                                                                   ccgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCg                 <  1:1924801/101‑1 (MQ=255)
                                                                                    cgccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGa                >  1:405037/1‑101 (MQ=255)
                                                                                     gccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAc               <  1:1021192/101‑1 (MQ=255)
                                                                                     cccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAc               <  2:1968122/100‑1 (MQ=255)
                                                                                      ccACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGAcc              <  2:2942022/101‑1 (MQ=255)
                                                                                       cACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCt             >  1:4141652/1‑101 (MQ=255)
                                                                                        aCCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTg            >  2:4375079/1‑101 (MQ=255)
                                                                                           gTCGCCGAAGCCAGCGCCGAGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCACGGCCACAGCAAGGCTGGAGAGCGCGCAGCCCTGAtg         >  2:3141982/1‑100 (MQ=37)
                                                                                             cGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCgg       >  1:549350/1‑101 (MQ=255)
                                                                                              gCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAACGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACAGCAAGGCTGGAGAGCGCCCCGACCTGATCggg      >  2:723540/1‑100 (MQ=255)
                                                                                                cGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGc    >  2:3573016/1‑101 (MQ=255)
                                                                                                 gAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCt   >  1:2358759/1‑101 (MQ=255)
                                                                                                  aaGCCAGCGCCGCGATCTGCGCCGTGCTCATCGCCGCCGCCTGAACCGAGCCAAGCCCGCCCACCGCAAGGCTGGAAAGCGCGCCGACCTGATCGGTGCTc  >  2:4440771/1‑101 (MQ=38)
                                                                                                  aaGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc  <  1:2221046/101‑1 (MQ=255)
                                                                                                  aaGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTc  <  1:701067/101‑1 (MQ=255)
                                                                                                  |                                                                                                    
CGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCTGGAGAGCGCGCCGACCTGATCGGTGCTC  >  NC_014034/681492‑681690

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