Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 678,785 T→A 7.9% S1608C (AGC→TGC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034678,7850T→A 7.9% 46.2 ‑0.7 ‑0.0239S1608C (AGC→TGC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (12/7):  ref base (115/105):  total (127/112)
Fisher's exact test strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.43e-01

Only 100 of 248 total aligned reads displayed.

GTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGC  >  NC_014034/678685‑678885
                                                                                                    |                                                                                                    
  gCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                                    <  1:1094182/101‑1 (MQ=255)
   cTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAg                                                                                                   <  1:2332254/101‑1 (MQ=255)
     cAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                                                 >  2:1482960/1‑101 (MQ=255)
       gtcTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgcc                                                                                               <  1:3068016/99‑1 (MQ=255)
       ggCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgcc                                                                                               <  1:787657/101‑1 (MQ=255)
        gCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccg                                                                                              <  2:3385354/101‑1 (MQ=255)
        gCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccg                                                                                              <  1:1887758/101‑1 (MQ=255)
          ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgag                                                                                            >  1:1039127/1‑99 (MQ=255)
          ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCcgcg                                                                                            >  1:3158048/1‑101 (MQ=255)
            ggCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGAt                                                                                          >  1:3419878/1‑101 (MQ=255)
             gCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATc                                                                                         >  1:1119536/1‑101 (MQ=255)
                tCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGc                                                                                      <  1:332145/101‑1 (MQ=255)
                 cAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCg                                                                                     <  1:3692535/101‑1 (MQ=255)
                   gcgcCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTc                                                                                   <  2:2823892/101‑1 (MQ=255)
                        gTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCt                                                                              >  2:1766575/1‑101 (MQ=255)
                        gTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCt                                                                              >  2:3077954/1‑101 (MQ=255)
                        gTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCt                                                                              >  1:2550491/1‑101 (MQ=255)
                          cccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                            <  1:1668365/101‑1 (MQ=255)
                           ccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAg                                                                           >  2:1226598/1‑101 (MQ=255)
                           ccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAg                                                                           <  1:2408117/101‑1 (MQ=255)
                            cGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGc                                                                          <  1:2071993/101‑1 (MQ=255)
                             gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                                         >  2:3772205/1‑101 (MQ=255)
                             gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                                         <  2:5328117/101‑1 (MQ=255)
                              aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                        >  2:1263822/1‑101 (MQ=255)
                              aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                        >  1:1601275/1‑101 (MQ=255)
                               tCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCgg                                                                       >  2:4350879/1‑100 (MQ=255)
                                ccccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTg                                                                      <  2:133631/101‑1 (MQ=255)
                                 cccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgc                                                                     <  2:54262/101‑1 (MQ=255)
                                  ccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgcc                                                                    <  2:3966706/101‑1 (MQ=255)
                                    gCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccgc                                                                  >  1:4395710/1‑101 (MQ=255)
                                       aTGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTg                                                               <  2:2701461/101‑1 (MQ=255)
                                       aTGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTg                                                               <  2:5612786/101‑1 (MQ=255)
                                       aTGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTg                                                               >  1:4035/1‑101 (MQ=255)
                                        tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                              >  2:5216658/1‑101 (MQ=255)
                                         ggCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgcg                                                             >  1:5270140/1‑101 (MQ=255)
                                         ggCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgcg                                                             <  2:2578948/101‑1 (MQ=255)
                                            cGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCg                                                          >  1:3984379/1‑101 (MQ=255)
                                             gagaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGt                                                         <  2:1640616/101‑1 (MQ=255)
                                              agaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTg                                                        >  1:3393043/1‑101 (MQ=255)
                                              agaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGAGCTCAGCCGTGCCGCCTGCGCCGTg                                                        >  2:3249800/1‑101 (MQ=255)
                                               gaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTgg                                                       >  1:301650/1‑101 (MQ=255)
                                                aCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGt                                                      >  2:3684590/1‑101 (MQ=255)
                                                 ccGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTc                                                     <  1:714316/101‑1 (MQ=255)
                                                  cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGGCGCCTGCGCCGTGGTCa                                                    >  2:5197856/1‑101 (MQ=255)
                                                   gAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAg                                                   >  2:5773114/1‑101 (MQ=255)
                                                   gAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAg                                                   >  1:3206146/1‑101 (MQ=255)
                                                    aaaGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgc                                                  >  2:4436652/1‑101 (MQ=255)
                                                      aGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgcgc                                                >  1:3776776/1‑101 (MQ=255)
                                                      aGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgcgc                                                >  1:5380993/1‑101 (MQ=255)
                                                        ccgcaccagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCa                                              >  1:2120853/10‑101 (MQ=2)
                                                        cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCa                                              >  1:1711652/1‑101 (MQ=255)
                                                         tGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGcccc                                             >  1:1361434/1‑99 (MQ=255)
                                                          gcaccagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                            >  1:1257668/8‑101 (MQ=2)
                                                           caccagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                           >  1:5514965/7‑101 (MQ=2)
                                                           caccagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                           >  1:3549038/7‑101 (MQ=2)
                                                             cGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATg                                         >  1:4754360/1‑101 (MQ=255)
                                                              gATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                                        >  2:4342354/1‑101 (MQ=255)
                                                              gATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                                        <  2:3556655/101‑1 (MQ=255)
                                                              cagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                        >  2:2604188/4‑101 (MQ=2)
                                                               atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                       >  2:4818359/1‑101 (MQ=255)
                                                               atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                       >  1:2872800/1‑101 (MQ=255)
                                                               agatCGGCGGTTTCCAGCGCCGCAACCTGCGTGTTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                       <  1:4934867/99‑1 (MQ=2)
                                                                gatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                      <  1:3954294/100‑1 (MQ=2)
                                                                  tcggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                                    >  2:5356616/4‑100 (MQ=2)
                                                                  tCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAtt                                    <  2:4853/101‑1 (MQ=255)
                                                                   cggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtcg                                   <  2:5539381/99‑3 (MQ=9)
                                                                   cggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtcg                                   >  1:1460551/3‑99 (MQ=9)
                                                                   ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTc                                   <  1:1760110/101‑1 (MQ=255)
                                                                    cgcgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCg                                  <  2:4030/101‑1 (MQ=255)
                                                                      cgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTg                                <  1:4180159/101‑1 (MQ=255)
                                                                      cgGTTTCCACCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTg                                >  2:4737789/1‑101 (MQ=255)
                                                                       ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                               <  1:5615790/101‑1 (MQ=255)
                                                                       ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                               <  1:1451377/101‑1 (MQ=255)
                                                                         tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                             >  2:4623711/1‑101 (MQ=255)
                                                                         tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                             >  1:5785543/1‑101 (MQ=255)
                                                                         tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                             >  2:1043965/1‑101 (MQ=255)
                                                                         tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                             <  2:3419673/101‑1 (MQ=255)
                                                                          ttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCa                            <  2:1119455/101‑1 (MQ=255)
                                                                          ttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCa                            <  2:4400097/101‑1 (MQ=255)
                                                                           tCCAGCGCCGCGATCTGAGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAg                           >  2:3879472/1‑101 (MQ=255)
                                                                                cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                      <  1:2639733/101‑1 (MQ=255)
                                                                                cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                      <  1:1949818/101‑1 (MQ=255)
                                                                                cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGCTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                      >  1:4980961/1‑101 (MQ=255)
                                                                                cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATATGCGTAGCGCTCAGCCGCGCCGCCTGCGCCGTGGTCAGCGCTACGATGTCATTCGTGCCCAGCCCCt                      >  2:5135487/1‑101 (MQ=25)
                                                                                 gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt                     <  2:631006/101‑1 (MQ=255)
                                                                                   cgcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCa                   >  2:3837890/1‑101 (MQ=255)
                                                                                   cgcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCa                   <  2:2444853/101‑1 (MQ=255)
                                                                                    gcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAg                  <  2:5074816/101‑1 (MQ=255)
                                                                                        tCTGCGCCGTGCTCAGCGCCGCGATCTGCGTAGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGcgcc              >  1:5238047/1‑101 (MQ=255)
                                                                                          tGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGcgccgc            <  2:1256772/101‑1 (MQ=255)
                                                                                            cgcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCgg          >  2:2379524/1‑101 (MQ=255)
                                                                                             gcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGt         >  2:672955/1‑101 (MQ=255)
                                                                                                gTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCt      >  2:1968942/1‑101 (MQ=255)
                                                                                                gTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCt      >  1:3878626/1‑101 (MQ=255)
                                                                                                 tGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTc     <  2:3403880/101‑1 (MQ=255)
                                                                                                 tGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTc     <  1:4818745/101‑1 (MQ=255)
                                                                                                  gCTCAGCGCCGCGATCTGCTTCGCGCTCAGCCGTGCTGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCa    <  2:2051311/101‑1 (MQ=255)
                                                                                                  gCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCa    <  2:2295105/101‑1 (MQ=255)
                                                                                                   cTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAg   <  2:3135554/101‑1 (MQ=255)
                                                                                                    tCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgc  >  2:300841/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGC  >  NC_014034/678685‑678885

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG