Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,444 C→G 8.0% 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 8.0% 22.6 ‑2.5 ‑0.0149A721A (GCG→GCC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/9):  ref base (72/65):  total (75/74)
Fisher's exact test strand distribution p-value = 7.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.28e-01

Only 100 of 166 total aligned reads displayed.

GCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTG  >  NC_014034/681344‑681543
                                                                                                    |                                                                                                   
gCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgcc                                                                                                     >  2:2015036/1‑101 (MQ=255)
gCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATccgcc                                                                                                     >  2:2516799/1‑101 (MQ=255)
  cgcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCcgccgc                                                                                                   <  1:1690664/101‑1 (MQ=255)
   gcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                  <  1:791619/101‑1 (MQ=255)
   gcgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCgccgcc                                                                                                  <  2:2449312/101‑1 (MQ=255)
    cgTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCt                                                                                                 <  1:4707354/101‑1 (MQ=255)
         aCAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAgc                                                                                            <  1:882672/101‑1 (MQ=255)
              gccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGcgccgc                                                                                       >  2:4128185/1‑101 (MQ=255)
               ccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCcgcg                                                                                      >  2:2003041/1‑101 (MQ=255)
               ccgcGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCcgcg                                                                                      <  1:1439781/101‑1 (MQ=255)
                 gcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGAt                                                                                    <  1:2166742/101‑1 (MQ=255)
                  cgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATc                                                                                   <  1:2028818/101‑1 (MQ=255)
                  cgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCAGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATc                                                                                   >  1:3303914/1‑101 (MQ=255)
                   gATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCt                                                                                  >  1:4579442/1‑101 (MQ=255)
                    aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                 >  1:4449247/1‑101 (MQ=255)
                     tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                >  1:2531617/1‑101 (MQ=255)
                      cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                               >  1:3274409/1‑101 (MQ=255)
                      cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                               >  1:4117742/1‑101 (MQ=255)
                      cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                               >  1:4341432/1‑101 (MQ=255)
                      cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                               >  2:509883/1‑101 (MQ=255)
                        gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                             <  1:1322777/101‑1 (MQ=255)
                               tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                      <  2:1013222/101‑1 (MQ=255)
                               tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                      <  2:1514218/101‑1 (MQ=255)
                                gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                     <  2:69896/101‑1 (MQ=255)
                                 aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                    >  1:3685386/1‑101 (MQ=255)
                                 aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                    <  1:1960564/101‑1 (MQ=255)
                                 aaaGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgc                                                                    <  1:4642062/101‑1 (MQ=255)
                                   aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                  >  1:4718585/1‑101 (MQ=255)
                                     cccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgcgc                                                                >  1:2381332/1‑101 (MQ=255)
                                      ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGcc                                                               >  2:3458996/1‑101 (MQ=255)
                                       cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                              <  2:3111290/101‑1 (MQ=255)
                                         gCCATGGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                            <  1:4668574/101‑1 (MQ=255)
                                         gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGAt                                                            <  1:2003158/101‑1 (MQ=255)
                                          cgaTCTTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                           <  1:1636810/99‑1 (MQ=255)
                                          ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATc                                                           >  1:2815165/1‑101 (MQ=255)
                                            aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                         >  1:1130143/1‑101 (MQ=255)
                                            aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTg                                                         <  1:391185/101‑1 (MQ=255)
                                            aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTAGTCAGCGCGCCGATCTg                                                         >  2:4265993/1‑101 (MQ=255)
                                             tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGagatcg                                                        <  2:4024827/101‑6 (MQ=255)
                                              cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcg                                                       >  2:4632271/1‑101 (MQ=255)
                                              cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcg                                                       >  2:1700033/1‑101 (MQ=255)
                                               gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTgcgc                                                      >  2:4254391/1‑101 (MQ=255)
                                               gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGCCAGCGCGCCGATCTgcgc                                                      >  1:114340/1‑101 (MQ=255)
                                                tGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATcagatcg                                                     >  2:1636710/1‑95 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGTCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                    <  1:811363/101‑1 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                    >  2:304737/1‑101 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCg                                                    >  1:1016974/1‑101 (MQ=255)
                                                 ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGAGCCg                                                    <  2:4017818/101‑1 (MQ=255)
                                                  gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                   >  1:2906909/1‑101 (MQ=255)
                                                  gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                   >  2:2654058/1‑101 (MQ=255)
                                                  gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGa                                                   <  1:3878016/101‑1 (MQ=255)
                                                    cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                 >  1:1709621/1‑101 (MQ=255)
                                                    cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGc                                                 >  2:2485234/1‑101 (MQ=255)
                                                     gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTTCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                <  2:2668232/101‑1 (MQ=255)
                                                     gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGcc                                                >  1:126515/1‑101 (MQ=255)
                                                       gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCa                                              <  1:1746313/101‑1 (MQ=255)
                                                       gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCa                                              >  2:1895284/1‑101 (MQ=255)
                                                            gcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgc                                         <  2:4718329/101‑1 (MQ=255)
                                                             cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgcg                                        >  2:3316963/1‑101 (MQ=255)
                                                             cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGAcgcg                                        >  2:3833035/1‑101 (MQ=255)
                                                              ccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACgcgc                                       <  1:3848647/101‑1 (MQ=255)
                                                                   gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCt                                  <  1:1614020/101‑1 (MQ=255)
                                                                   gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCt                                  <  2:4529032/101‑1 (MQ=255)
                                                                     gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTgc                                >  1:2628432/1‑101 (MQ=255)
                                                                        ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                             <  1:4343181/101‑1 (MQ=255)
                                                                         cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGACCCCAGACGCGCCGCCTGCGCcg                            >  2:3100594/1‑101 (MQ=255)
                                                                           tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgcc                          <  2:4579189/101‑3 (MQ=255)
                                                                            cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                         <  1:3983329/101‑1 (MQ=255)
                                                                            cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGCGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcg                         >  1:3789320/1‑101 (MQ=255)
                                                                              cggtGATTTGCGCCAGTTCCGCCGCCTCGAGCGCCGCTATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgttgtc                       <  2:2435431/98‑1 (MQ=21)
                                                                              cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                       >  1:2637116/1‑101 (MQ=255)
                                                                               gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                      >  1:4658288/1‑101 (MQ=255)
                                                                               gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                      >  2:1369879/1‑101 (MQ=255)
                                                                               gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                      <  2:3293948/101‑1 (MQ=255)
                                                                                gggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAg                     <  2:3790946/91‑1 (MQ=2)
                                                                                 tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                    <  2:4753330/101‑1 (MQ=255)
                                                                                 tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                    >  2:1672570/1‑101 (MQ=255)
                                                                                 tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                    >  2:359491/1‑101 (MQ=255)
                                                                                 ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                    <  2:3440122/92‑1 (MQ=2)
                                                                                  gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                   <  2:1264514/93‑1 (MQ=2)
                                                                                  gagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgcg                   <  1:3384721/93‑1 (MQ=2)
                                                                                  gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                   >  1:666851/1‑101 (MQ=255)
                                                                                  gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                   >  2:3848803/1‑101 (MQ=255)
                                                                                    gccgccccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                 <  2:2086493/95‑1 (MQ=2)
                                                                                    gccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGcc                 <  2:1131612/95‑1 (MQ=2)
                                                                                     cTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                <  2:2906733/101‑1 (MQ=255)
                                                                                      tGCGCCAGATCCGCCGCCTCGAGCGCCGCGGTCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc               <  2:4594053/101‑1 (MQ=255)
                                                                                      cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc               >  2:1236656/5‑101 (MQ=2)
                                                                                       gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg              <  2:631363/98‑1 (MQ=2)
                                                                                        caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa             <  2:3302928/99‑1 (MQ=9)
                                                                                         accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt            <  1:984012/100‑1 (MQ=2)
                                                                                         accAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt            >  1:1243664/2‑101 (MQ=2)
                                                                                          ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg           <  2:666787/101‑1 (MQ=255)
                                                                                          ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATg           <  1:4414092/101‑1 (MQ=255)
                                                                                           cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt          >  1:2452113/1‑101 (MQ=255)
                                                                                            aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc         >  1:2691703/1‑101 (MQ=255)
                                                                                            aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc         >  1:2873626/1‑101 (MQ=255)
                                                                                            aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc         >  2:2433586/1‑101 (MQ=255)
                                                                                                 cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg    >  1:1219475/1‑101 (MQ=255)
                                                                                                   ccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg  >  2:2750991/1‑101 (MQ=255)
                                                                                                    |                                                                                                   
GCCGCGTCGACAGCGCCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTG  >  NC_014034/681344‑681543

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