Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,485 T→G 18.0% M708L (ATC→CTC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4850T→G18.0% 197.4 ‑3.0 ‑0.1414M708L (ATC→CTC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (23/51):  ref base (147/191):  total (171/243)
Fisher's exact test strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.67e-01

Only 100 of 441 total aligned reads displayed.

GCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTC  >  NC_014034/681385‑681585
                                                                                                     |                                                                                                     
 ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.c                                                                                                     <  1:791302/101‑1 (MQ=255)
  cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.Ct                                                                                                    <  2:8467935/101‑1 (MQ=255)
    tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgc                                                                                                  <  1:4743009/101‑1 (MQ=255)
    tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgc                                                                                                  <  2:1826129/101‑1 (MQ=255)
    tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTgc                                                                                                  <  1:8739075/101‑1 (MQ=255)
       tGGCCGAGCCCAGCCCGGAGCGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGcc                                                                                               <  1:8820805/101‑1 (MQ=255)
        ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCg                                                                                              <  1:5125290/101‑1 (MQ=255)
         gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGa                                                                                             >  1:5540303/1‑101 (MQ=255)
         gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGa                                                                                             >  1:8371467/1‑101 (MQ=255)
            gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGcc                                                                                          >  2:3558256/1‑101 (MQ=255)
            gAGCCCAGCCCGGAGAGCGCCGTAGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCACCGGGCCG.AT.CTTGGCCGAGca                                                                                          >  2:2963427/1‑100 (MQ=38)
               cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAg                                                                                       <  2:7279503/101‑1 (MQ=255)
                  agccaGGGGAGCGCCGTGGCGCTGATCTGCGCCCGATCCGCCGCTTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGAcg                                                                                    <  2:6989282/101‑1 (MQ=37)
                      cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgcc                                                                                >  2:2264079/1‑101 (MQ=255)
                      cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgcc                                                                                >  2:7335782/1‑101 (MQ=255)
                       ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGcgccg                                                                               <  2:7650308/101‑1 (MQ=255)
                         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCgccgcc                                                                             <  2:6083718/101‑1 (MQ=255)
                         agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCgccgcc                                                                             >  1:1468991/1‑101 (MQ=255)
                          gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCt                                                                            <  2:1941648/101‑1 (MQ=255)
                          gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCt                                                                            <  2:7119944/101‑1 (MQ=255)
                          gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCt                                                                            >  2:908490/1‑101 (MQ=255)
                           agCGGCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                           <  1:2916031/101‑1 (MQ=255)
                           agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTg                                                                           >  2:138578/1‑101 (MQ=255)
                           agCGCCGTCGCGCTGATCTGCGCCAGACCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.GTGCGCCGTGCCCAGACGCGCCGCCTg                                                                           <  2:2518776/101‑1 (MQ=255)
                                  tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                                    <  1:3834428/101‑1 (MQ=255)
                                  gtggtgtGATCCGCACCATATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCcgtc                                                                    <  2:8399893/95‑1 (MQ=2)
                                    gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtcgt                                                                  >  2:7767872/1‑101 (MQ=255)
                                      gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                                >  1:633313/1‑101 (MQ=255)
                                        tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCcggggcattt                                                              <  2:3788813/101‑9 (MQ=255)
                                        tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                              >  1:2064066/1‑101 (MQ=255)
                                         gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                             >  2:1026708/1‑101 (MQ=255)
                                         gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                             >  2:4424966/1‑101 (MQ=255)
                                         gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                                             <  2:3935649/101‑1 (MQ=255)
                                          aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                            >  2:6534962/1‑101 (MQ=255)
                                          aTCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcgc                                                            >  2:5707473/1‑101 (MQ=255)
                                           tCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGcc                                                           <  2:6374336/101‑1 (MQ=255)
                                            ccgcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCa                                                          <  2:4549563/99‑1 (MQ=255)
                                             tGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                         >  1:4371361/1‑101 (MQ=255)
                                             ccggcCAGATCCGCCGCCTCGAGCGCCGCGCTCTGCGTGCTGGGCAGCGCGCCG.AT.CTGCGCCGAGCCCAGCCGCGCCGCCTGCGCCGTCGTCAGCGCCAc                                                         <  2:8853564/98‑1 (MQ=21)
                                              gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                        <  1:4492502/101‑1 (MQ=255)
                                              gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTAGTCAGCGCCACg                                                        <  1:7476520/101‑1 (MQ=255)
                                              gcaccAGATCGGCGGGTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACg                                                        <  1:2961108/98‑1 (MQ=9)
                                                gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                                      <  1:192119/101‑1 (MQ=255)
                                                  cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                    <  1:7957285/101‑1 (MQ=255)
                                                  cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                    <  2:4662058/101‑1 (MQ=255)
                                                  cAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                    >  1:4947848/1‑101 (MQ=255)
                                                  cAGATCCGCCGCCTCACGCGCCGCGATCTGCTGACTGGATAGCGCGCCG.AT.CTGCGCCGAGCCCTGCTGTGCCGTCTGCGGCGTCGTCTGTGCCACGaaca                                                    >  2:972495/1‑98 (MQ=16)
                                                    gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                  >  2:7908291/1‑101 (MQ=255)
                                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                 <  1:2453450/101‑1 (MQ=255)
                                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                 <  2:5540096/101‑1 (MQ=255)
                                                     aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                                 <  1:6381518/101‑1 (MQ=255)
                                                      tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                <  1:3091636/101‑1 (MQ=255)
                                                      tCCCCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                                <  1:1011637/101‑1 (MQ=255)
                                                       ccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCg                                               >  1:9171611/1‑101 (MQ=255)
                                                        ggcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                              <  2:2177711/91‑1 (MQ=9)
                                                         gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgga                                             >  1:767248/1‑100 (MQ=255)
                                                           ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                           <  2:6155462/94‑1 (MQ=9)
                                                           ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                           >  1:584162/8‑101 (MQ=9)
                                                           ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                                           >  1:7121834/8‑101 (MQ=9)
                                                            gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                          <  1:6932790/101‑1 (MQ=255)
                                                            gccGCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                                          <  1:1074174/101‑1 (MQ=255)
                                                             ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                                         <  1:574162/101‑1 (MQ=255)
                                                              cTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaa                                        >  1:3240182/1‑101 (MQ=255)
                                                               tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                       <  2:8558090/101‑1 (MQ=255)
                                                               tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCCCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                                       >  2:3681176/1‑101 (MQ=255)
                                                                ccaGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgagg                                      >  1:1770656/3‑97 (MQ=2)
                                                                  aGCGCCGCGATCTGCGTGCTGGTCGGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGcc                                    <  2:5841875/101‑1 (MQ=255)
                                                                  aGCGCCGCACCCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                    <  1:705644/101‑1 (MQ=9)
                                                                  aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAG.AAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                                    >  1:1259597/1‑101 (MQ=9)
                                                                   gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccc                                   >  2:3094648/1‑100 (MQ=2)
                                                                   gccgcgtgATCTGCGTGCCGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcg                                   <  2:2923858/94‑1 (MQ=255)
                                                                    cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                                  >  1:8172284/1‑99 (MQ=2)
                                                                     gccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgcg                                 >  2:7461390/1‑101 (MQ=255)
                                                                     gccgcAACCTGCGTGCTCGACAGCGCCAGA.AG.CTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                                 <  1:8539796/101‑4 (MQ=9)
                                                                      ccgcGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCgg                                <  2:5084072/101‑1 (MQ=255)
                                                                        gcgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt                              <  1:2544415/101‑1 (MQ=255)
                                                                         cgATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc                             <  2:2037065/101‑1 (MQ=255)
                                                                          gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                            <  1:5461293/101‑1 (MQ=255)
                                                                          gATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg                            >  2:1229520/1‑101 (MQ=255)
                                                                           aTCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                           >  2:8336724/1‑101 (MQ=255)
                                                                           aTCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgc                           >  2:1696301/1‑101 (MQ=255)
                                                                           aTCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGGGCCACAGCCGCGGCTccc                           >  1:4285097/1‑99 (MQ=255)
                                                                                cGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg                      <  2:5211326/101‑1 (MQ=255)
                                                                                cGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGccgccg                      >  1:9020362/1‑101 (MQ=255)
                                                                                 gTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt                     >  1:2328972/1‑101 (MQ=255)
                                                                                    cTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCCCCGCCTGCGCCGTCGTCATCGCCACGATATCATCGGTGCCAAACCCACCGAGcaccgcgccgcc                  >  2:4257220/1‑95 (MQ=18)
                                                                                     tGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc                 <  1:3609888/101‑1 (MQ=255)
                                                                                      ggTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                >  1:3935235/1‑101 (MQ=255)
                                                                                      ggTCAGCGCGCCA.AT.CTGCGCAGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa                <  2:4884992/101‑1 (MQ=255)
                                                                                       gTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAg               >  2:7877737/1‑101 (MQ=255)
                                                                                        tCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc              <  1:6067648/101‑1 (MQ=255)
                                                                                         cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg             <  1:3777312/101‑1 (MQ=255)
                                                                                         cAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcg             >  1:8788796/1‑101 (MQ=255)
                                                                                          ggcgcgcCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcgc            <  1:7335117/100‑1 (MQ=255)
                                                                                                 cgaaT.CTGCGCCGGGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAc      <  2:2684192/98‑1 (MQ=255)
                                                                                                 cG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc     <  2:810616/101‑1 (MQ=255)
                                                                                                 cG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc     <  2:8187661/101‑1 (MQ=255)
                                                                                                 cG.AT.CGGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc     <  1:5317330/101‑1 (MQ=255)
                                                                                                  g.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCg    <  1:7785506/101‑1 (MQ=255)
                                                                                                     t.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTc  <  2:2090684/101‑1 (MQ=255)
                                                                                                     |                                                                                                     
GCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.AT.CTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTC  >  NC_014034/681385‑681585

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