Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,471 G→C 20.8% T712T (ACC→ACG RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4710G→C20.8% 164.5 ‑1.8 ‑1.8321T712T (ACC→ACGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (30/37):  ref base (138/116):  total (168/153)
Fisher's exact test strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.75e-02

Only 100 of 333 total aligned reads displayed.

CGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCA  >  NC_014034/681371‑681571
                                                                                                    |                                                                                                     
cgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTg                                                                                                       <  2:3481604/101‑1 (MQ=255)
 gcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTgg                                                                                                      >  1:6517200/1‑101 (MQ=255)
   cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                                                    <  1:4262922/101‑1 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                   >  1:4787788/1‑101 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                   <  2:1293308/101‑1 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                   >  1:2033095/1‑101 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                   >  2:853969/1‑101 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                   >  1:1507404/1‑101 (MQ=255)
    tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                                   >  2:6444112/1‑101 (MQ=255)
     gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAg                                                                                                  <  2:3901702/101‑1 (MQ=255)
        aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGGGCTGGTCAgcgc                                                                                               <  1:6320364/101‑1 (MQ=255)
           ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGcc                                                                                            <  1:2913150/101‑1 (MQ=255)
            cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                                                           >  1:4055435/1‑101 (MQ=255)
            cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                                                           >  1:562620/1‑101 (MQ=255)
             ggCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.a                                                                                         >  2:2623441/1‑101 (MQ=255)
                 aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                                     >  2:1507793/1‑101 (MQ=255)
                 aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                                     >  2:5894934/1‑101 (MQ=255)
                   cGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgcg                                                                                   >  2:1328201/1‑101 (MQ=255)
                    gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgcgc                                                                                  <  1:1654731/101‑1 (MQ=255)
                    gTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgcgc                                                                                  >  1:60321/1‑101 (MQ=255)
                      ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCg                                                                                <  1:3905556/101‑1 (MQ=255)
                         cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                             >  1:3396797/1‑101 (MQ=255)
                         cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                             >  2:1722509/1‑101 (MQ=255)
                         cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                             >  2:4129041/1‑101 (MQ=255)
                           agcccagcccGGAGAGCGCCGTCGCGTTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGccc                                                                           <  2:4627686/101‑1 (MQ=255)
                            gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCa                                                                          >  2:2498268/1‑101 (MQ=255)
                            gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCa                                                                          <  1:4406728/101‑1 (MQ=255)
                            gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCa                                                                          <  1:235004/101‑1 (MQ=255)
                             cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAg                                                                         <  1:1044410/101‑1 (MQ=255)
                              ccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGa                                                                        >  1:3007125/1‑101 (MQ=255)
                                     ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgccg                                                                 >  2:5052794/1‑101 (MQ=255)
                                      gagaGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGcgccgc                                                                >  1:732015/1‑101 (MQ=255)
                                       agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCgccgcc                                                               <  2:2083177/101‑1 (MQ=255)
                                        gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCt                                                              >  1:6726803/1‑101 (MQ=255)
                                         agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTg                                                             <  1:4129236/101‑1 (MQ=255)
                                         agCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTg                                                             <  1:3795200/101‑1 (MQ=255)
                                          gcgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCCCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTgc                                                            <  1:3699806/101‑1 (MQ=255)
                                             ccGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGcc                                                         <  1:3138634/101‑1 (MQ=255)
                                                tCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCcgtc                                                      <  2:311680/101‑1 (MQ=255)
                                                   cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                   >  2:211815/1‑101 (MQ=255)
                                                   cgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCgtcgtc                                                   <  1:5917029/101‑1 (MQ=255)
                                                    gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                  >  2:5641927/1‑101 (MQ=255)
                                                    gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCa                                                  <  2:6743848/101‑1 (MQ=255)
                                                      tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgc                                                <  1:6094901/101‑1 (MQ=255)
                                                      ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAgc                                                <  2:4450011/92‑1 (MQ=2)
                                                       gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                               >  1:5960451/1‑101 (MQ=255)
                                                       cagccgcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAgcg                                               >  2:5820085/6‑101 (MQ=255)
                                                           cgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCAc                                           <  2:5526735/97‑1 (MQ=2)
                                                              gcCAGCTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGGCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                        <  1:2733073/101‑1 (MQ=255)
                                                              gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                        >  2:632557/1‑101 (MQ=255)
                                                              gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                        <  1:529540/101‑1 (MQ=255)
                                                               ccAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGGCAGCGCCACGATg                                       >  2:7057178/1‑101 (MQ=255)
                                                                 aGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGGCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTTCGCCGTCGCACGCGTCTCGATGTc                                     >  2:1110867/1‑101 (MQ=25)
                                                                  gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                    >  1:3313505/1‑101 (MQ=255)
                                                                   aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCCCGATGTCAt                                   >  2:4343735/1‑101 (MQ=255)
                                                                   aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                                   <  1:2498059/101‑1 (MQ=255)
                                                                    tCCGCCCCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCAAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                  <  1:7029558/101‑1 (MQ=255)
                                                                    cccgccgCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                                  <  1:991507/100‑1 (MQ=255)
                                                                      cgcggcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                <  2:5009225/101‑1 (MQ=255)
                                                                      cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                >  1:1815489/1‑101 (MQ=255)
                                                                      cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCgg                                >  2:616016/1‑101 (MQ=255)
                                                                      cgccgcCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGGCATCgg                                >  2:4027826/1‑101 (MQ=255)
                                                                       gcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                               <  1:1772036/92‑1 (MQ=9)
                                                                       gccgccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCGTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                               >  1:4740623/1‑101 (MQ=255)
                                                                        cggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTg                              >  1:2628802/9‑101 (MQ=9)
                                                                         ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                             <  2:4929670/94‑1 (MQ=9)
                                                                         ggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGc                             >  2:2907028/8‑101 (MQ=9)
                                                                          gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                            >  2:3307890/7‑101 (MQ=9)
                                                                          gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                            <  1:3981151/95‑1 (MQ=9)
                                                                          gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                            <  2:1879065/95‑1 (MQ=12)
                                                                          gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACTCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                            <  2:3394130/101‑1 (MQ=255)
                                                                           ccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCa                           >  2:5983090/1‑101 (MQ=255)
                                                                            ttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccga                          >  1:5463988/5‑99 (MQ=2)
                                                                             tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                         <  1:6960136/98‑4 (MQ=9)
                                                                             tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag                         <  2:4850458/98‑4 (MQ=9)
                                                                             tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                         <  1:6253503/101‑1 (MQ=255)
                                                                             tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa                         >  2:2378513/1‑101 (MQ=255)
                                                                               caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc                       >  1:1293356/2‑96 (MQ=9)
                                                                                aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGC.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                      <  2:5259344/101‑1 (MQ=9)
                                                                                aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                      >  2:3525382/1‑101 (MQ=2)
                                                                                aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc                      >  1:3492027/1‑101 (MQ=2)
                                                                                 gcgcCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGCGCCGAGGccc                     <  1:6624869/101‑2 (MQ=9)
                                                                                  cgccgcGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCcgc                    >  1:167022/1‑101 (MQ=255)
                                                                                  cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct                    <  1:5089339/101‑3 (MQ=9)
                                                                                   gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt                   >  1:1261454/1‑98 (MQ=2)
                                                                                     cgcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa                 >  1:4432873/1‑101 (MQ=255)
                                                                                        tgtCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg              <  2:6922615/99‑1 (MQ=255)
                                                                                        gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg              >  1:678763/1‑101 (MQ=255)
                                                                                        gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAGcg              >  1:2271622/1‑101 (MQ=255)
                                                                                          tCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc            <  2:540994/101‑1 (MQ=255)
                                                                                          tCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc            <  2:3007016/101‑1 (MQ=255)
                                                                                          tCTGCGTGCGGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc            <  1:5698647/101‑1 (MQ=255)
                                                                                          tCTCCGCGCTGGTCAGCGCACCG.ATCTGCGCCGAGCCCAGACGCCCCGCCTCCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgcc            <  1:2767074/101‑1 (MQ=37)
                                                                                          ccTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGTCACGATGTCATCGGTGCCGAGGCACCTGATcgca            >  2:2222327/2‑100 (MQ=2)
                                                                                           cTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg           >  2:63945/1‑101 (MQ=255)
                                                                                           cTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccg           <  1:4957250/101‑1 (MQ=255)
                                                                                           cTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCACCGGTGCCGAGGCCCTTGATcgcat           >  1:4553372/1‑99 (MQ=2)
                                                                                             gCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCgccgcc         <  2:370777/101‑1 (MQ=255)
                                                                                               gggCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGt       <  1:131929/99‑1 (MQ=255)
                                                                                                    ggTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCa  <  1:63946/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
CGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCA  >  NC_014034/681371‑681571

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