Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,444 C→G 7.3% A721A (GCG→GCC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4440C→G 7.3% 24.3 ‑0.6 ‑0.0177A721A (GCG→GCC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (4/9):  ref base (69/95):  total (73/104)
Fisher's exact test strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.52e-01

Only 100 of 205 total aligned reads displayed.

GCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGC  >  NC_014034/681344‑681544
                                                                                                    |                                                                                                    
gCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgcc                                                                                                      <  2:1785077/101‑1 (MQ=255)
   gcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                   >  2:4550294/1‑101 (MQ=255)
   gcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                   >  2:2124225/1‑101 (MQ=255)
    cgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCt                                                                                                  <  1:3417035/101‑1 (MQ=255)
     gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                 <  2:4421479/101‑1 (MQ=255)
     gTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTc                                                                                                 <  1:4506278/101‑1 (MQ=255)
      tCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCg                                                                                                <  1:1752747/101‑1 (MQ=255)
            gcgcCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgcc                                                                                          <  1:270697/101‑1 (MQ=255)
             cgccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTTATCTGCGCCAGATCCGCCGCCTCGAGcgccg                                                                                         <  1:4222297/101‑1 (MQ=255)
             cgccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccg                                                                                         <  1:2290425/101‑1 (MQ=255)
                 gcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGAt                                                                                     <  1:2478933/101‑1 (MQ=255)
                 gcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGAt                                                                                     >  2:2242969/1‑101 (MQ=255)
                  cgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATc                                                                                    <  1:2026569/101‑1 (MQ=255)
                  cgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATc                                                                                    >  2:1030188/1‑101 (MQ=255)
                   gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                   <  2:1097922/101‑1 (MQ=255)
                    aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                  >  1:1803693/1‑101 (MQ=255)
                    aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                  >  2:3852216/1‑101 (MQ=255)
                     tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                 >  1:3226712/1‑101 (MQ=255)
                     gcggcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                 <  2:3369022/98‑1 (MQ=255)
                          gcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                            <  2:1989613/101‑1 (MQ=255)
                           cgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTg                                                                           >  2:883765/1‑101 (MQ=255)
                            gcCTGAAGGCCCGGCCACCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTgg                                                                          <  2:4184701/101‑1 (MQ=255)
                              cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCCCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                        <  1:768589/101‑1 (MQ=255)
                               tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                       <  1:4167438/101‑1 (MQ=255)
                               tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                       <  2:45794/101‑1 (MQ=255)
                               tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                       >  1:1063646/1‑101 (MQ=255)
                                gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                      <  2:3312655/101‑1 (MQ=255)
                                gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                      <  1:3183882/101‑1 (MQ=255)
                                 aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                     >  2:2433479/1‑101 (MQ=255)
                                  aaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcg                                                                    <  2:1263023/101‑1 (MQ=255)
                                  aaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcg                                                                    <  1:4505691/101‑1 (MQ=255)
                                   aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                   >  2:2583400/1‑101 (MQ=255)
                                     cccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgcgc                                                                 <  1:1894531/101‑1 (MQ=255)
                                      ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGcc                                                                <  2:1083541/101‑1 (MQ=255)
                                       cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                               >  2:678213/1‑101 (MQ=255)
                                         gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                             >  2:863654/1‑101 (MQ=255)
                                         gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                             <  2:4220612/101‑1 (MQ=255)
                                          ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTCCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                            <  2:487485/101‑1 (MQ=255)
                                           cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCt                                                           <  2:376051/101‑1 (MQ=255)
                                           cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCt                                                           >  2:4324546/1‑101 (MQ=255)
                                              cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcg                                                        <  1:2585586/101‑1 (MQ=255)
                                               gTGGCGGAGCCCAGCCCGGAGAGCGCCGTCGCGGTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcgc                                                       <  1:3998359/101‑1 (MQ=255)
                                               gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcgc                                                       >  2:1675517/1‑101 (MQ=255)
                                                tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGcc                                                      <  1:130423/101‑1 (MQ=255)
                                                tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGcc                                                      >  2:288156/1‑101 (MQ=255)
                                                tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGcc                                                      >  1:2239848/1‑101 (MQ=255)
                                                tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGcc                                                      <  2:3627433/101‑1 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                     >  1:3503476/1‑101 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                     <  1:4434042/101‑1 (MQ=255)
                                                  gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                    <  2:543501/101‑1 (MQ=255)
                                                    cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                  >  1:3873790/1‑101 (MQ=255)
                                                     gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                 <  2:2081011/101‑1 (MQ=255)
                                                     gAGCCCAGACCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGNGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                 >  1:1426858/1‑101 (MQ=255)
                                                      agcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGccc                                                >  2:3802820/1‑101 (MQ=255)
                                                       gcccagcccGTAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCa                                               <  2:2613461/101‑1 (MQ=255)
                                                          cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAc                                            <  2:23854/101‑1 (MQ=255)
                                                          cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAc                                            <  2:3572486/101‑1 (MQ=255)
                                                          cagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAc                                            >  2:1492412/1‑101 (MQ=255)
                                                           agcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcg                                           <  2:3460454/101‑1 (MQ=255)
                                                               cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgcc                                       <  2:3065080/101‑1 (MQ=255)
                                                               cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgcc                                       >  2:473905/1‑101 (MQ=255)
                                                               cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgcc                                       >  1:424765/1‑101 (MQ=255)
                                                                ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgccg                                      >  1:4200161/1‑101 (MQ=255)
                                                                ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgccg                                      <  1:2775615/101‑1 (MQ=255)
                                                                 gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGcgccgc                                     <  2:2897994/101‑1 (MQ=255)
                                                                  agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCgccgcc                                    >  1:2995339/1‑101 (MQ=255)
                                                                  agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCgccgcc                                    <  1:83661/101‑1 (MQ=255)
                                                                   gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCt                                   <  2:2739261/101‑1 (MQ=255)
                                                                    agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTg                                  >  2:281566/1‑101 (MQ=255)
                                                                     gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgc                                 <  2:1198459/101‑1 (MQ=255)
                                                                      cgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgcg                                <  2:2278728/101‑1 (MQ=255)
                                                                         cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcg                             <  1:2253307/101‑1 (MQ=255)
                                                                         cGTCCCGCTGAGCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCCCGCCGATCTTCCCCGAGCCCAGACGCGCCGCCTGCGCcg                             <  2:1159285/101‑1 (MQ=25)
                                                                          gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgt                            >  2:124797/1‑101 (MQ=255)
                                                                              cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                        <  1:3797802/101‑1 (MQ=255)
                                                                               gcTGATCTGCGCCAGCTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                       <  1:4392501/101‑1 (MQ=255)
                                                                               gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                       >  1:2072566/1‑101 (MQ=255)
                                                                                cggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAg                      <  1:1383044/91‑1 (MQ=2)
                                                                                cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAg                      <  2:3355794/101‑1 (MQ=255)
                                                                                 tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                     <  2:3964048/101‑1 (MQ=255)
                                                                                   aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTTCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                   >  1:1087011/1‑101 (MQ=255)
                                                                                   aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGAGCCGAGCCCAGACGCGCCGGCTGCCACGTCGTCAgcgc                   >  2:2780064/1‑101 (MQ=25)
                                                                                    tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                  <  1:3668210/101‑1 (MQ=255)
                                                                                    gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                  <  2:2861264/95‑1 (MQ=2)
                                                                                    gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                  <  1:3896678/95‑1 (MQ=2)
                                                                                     cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                 <  2:3555257/101‑1 (MQ=255)
                                                                                     cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                 <  1:2866079/101‑1 (MQ=255)
                                                                                      tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                <  2:1833488/101‑1 (MQ=255)
                                                                                      cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                <  1:1572588/97‑1 (MQ=2)
                                                                                       gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg               <  1:513803/98‑1 (MQ=2)
                                                                                         cccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt             <  1:870736/100‑1 (MQ=255)
                                                                                           cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt           >  1:1467419/1‑101 (MQ=255)
                                                                                             gacccgccgCCCCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa         <  1:1641566/98‑1 (MQ=255)
                                                                                             gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa         >  2:4289623/1‑101 (MQ=255)
                                                                                              aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGACAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTcct        >  2:2318498/1‑99 (MQ=255)
                                                                                               tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc       <  1:1105933/101‑1 (MQ=255)
                                                                                                 cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg     >  1:2333012/1‑101 (MQ=255)
                                                                                                  gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCTCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt    >  2:2964995/1‑101 (MQ=255)
                                                                                                   ccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg   <  1:4330497/101‑1 (MQ=255)
                                                                                                    cgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc  <  1:3462265/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGC  >  NC_014034/681344‑681544

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