Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 678,773 T→C 5.6% M1612V (ATC→GTC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034678,7730T→C 5.6% 21.1 ‑3.2 ‑0.2230M1612V (ATC→GTC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/10):  ref base (120/97):  total (123/107)
Fisher's exact test strand distribution p-value = 4.13e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.85e-01

Only 100 of 237 total aligned reads displayed.

CGGATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCC  >  NC_014034/678673‑678873
                                                                                                    |                                                                                                    
 gggtCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                     <  2:2351655/98‑1 (MQ=255)
 ggATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                     >  1:3231164/1‑101 (MQ=255)
 ggATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                     <  1:1395432/101‑1 (MQ=255)
 ggATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                     >  1:2620618/1‑101 (MQ=255)
 ggATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                     >  1:2813460/1‑101 (MQ=255)
 ggATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATc                                                                                                     >  1:1755867/1‑101 (MQ=255)
   aTCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTg                                                                                                   >  2:2986249/1‑101 (MQ=255)
   aTCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTg                                                                                                   >  1:2137239/1‑101 (MQ=255)
    tCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgc                                                                                                  <  2:133541/101‑1 (MQ=255)
     cTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCCAGACCGAAAGATGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgcg                                                                                                 >  1:2409025/1‑101 (MQ=255)
      tGCGCCGTGCTCAGGCTTGCCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTgcgc                                                                                                <  2:577484/101‑1 (MQ=255)
        cgcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCg                                                                                              <  2:2216646/101‑1 (MQ=255)
        cgcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCg                                                                                              >  2:4289206/1‑101 (MQ=255)
        cgcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCg                                                                                              >  2:767710/1‑101 (MQ=255)
         gcCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCTGt                                                                                             <  1:2706037/101‑1 (MQ=255)
          ccGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTg                                                                                            >  1:2283981/1‑101 (MQ=255)
          ccGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTg                                                                                            >  1:2516956/1‑101 (MQ=255)
           cGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGc                                                                                           >  2:3877543/1‑101 (MQ=255)
            gTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGTCGTGCt                                                                                          >  1:769756/1‑101 (MQ=255)
             tGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTc                                                                                         >  2:2538418/1‑101 (MQ=255)
              gCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                        >  2:4411419/1‑101 (MQ=255)
              gCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                        >  2:1946866/1‑101 (MQ=255)
               cTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAg                                                                                       <  2:1704511/101‑1 (MQ=255)
                 cAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                                     >  2:3619140/1‑101 (MQ=255)
                  aGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                                                                    >  2:633853/1‑101 (MQ=255)
                  aGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcgc                                                                                    >  2:1234830/1‑101 (MQ=255)
                     cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                                 >  2:1385612/1‑101 (MQ=255)
                      ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCcgcg                                                                                >  1:1576930/1‑101 (MQ=255)
                      ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCcgcg                                                                                >  2:4514470/1‑101 (MQ=255)
                      ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCcgcg                                                                                >  2:4663604/1‑101 (MQ=255)
                      ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCcgcg                                                                                >  1:3630944/1‑101 (MQ=255)
                        ggctcccgccccGTCCCGGTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCCATCTGCGCCGTGCTCAGCGCCGCGAt                                                                              <  2:744785/91‑1 (MQ=38)
                         gCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATc                                                                             <  1:3489431/101‑1 (MQ=255)
                           aTCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTg                                                                           >  2:822764/1‑101 (MQ=255)
                              aGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGt                                                                        <  2:771817/101‑1 (MQ=255)
                               gcgcCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTc                                                                       <  2:3967457/101‑1 (MQ=255)
                                    gTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCt                                                                  <  2:3231938/101‑1 (MQ=255)
                                    gTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCt                                                                  >  2:4325308/1‑101 (MQ=255)
                                      cccGATCCCCGCCAGGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTGTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                <  2:3142842/101‑1 (MQ=255)
                                      cccGAGCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                <  2:1429690/101‑1 (MQ=255)
                                        cGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGc                                                              >  2:342944/1‑101 (MQ=255)
                                           tCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGt                                                           >  1:1281813/1‑101 (MQ=255)
                                              ccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgcc                                                        <  2:405332/101‑1 (MQ=255)
                                              ccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgcc                                                        <  2:3310275/101‑1 (MQ=255)
                                               cGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccg                                                       >  1:2742695/1‑101 (MQ=255)
                                                 ccATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCACTCAGCCGTgccgcc                                                     <  1:2607950/101‑1 (MQ=255)
                                                    tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                  >  1:1384655/1‑101 (MQ=255)
                                                    tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                  <  2:2408890/101‑1 (MQ=255)
                                                     ggCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgcg                                                 >  1:655924/1‑101 (MQ=255)
                                                        cGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCg                                              >  2:4596139/1‑101 (MQ=255)
                                                        cGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCg                                              >  1:2569612/1‑101 (MQ=255)
                                                          agaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTg                                            <  1:1401538/101‑1 (MQ=255)
                                                             ccGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTc                                         <  2:4163993/101‑1 (MQ=255)
                                                              cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                        <  1:4543010/101‑1 (MQ=255)
                                                               gAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAg                                       <  1:4289453/101‑1 (MQ=255)
                                                                aaaGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgc                                      >  2:2667629/1‑101 (MQ=255)
                                                                  aGCTTGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgcgc                                    >  1:235204/1‑101 (MQ=255)
                                                                     tGGGCGATATCCGCGGTTTCCAGCGCCGCTATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCAc                                 >  1:3698057/1‑101 (MQ=255)
                                                                     tGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCAc                                 >  1:4065907/1‑101 (MQ=255)
                                                                     tGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCGCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCAc                                 >  1:4047868/1‑101 (MQ=255)
                                                                      gggCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACg                                >  2:3848503/1‑101 (MQ=255)
                                                                        gCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGAt                              <  1:4735033/101‑1 (MQ=255)
                                                                        gCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGAt                              <  1:2683092/101‑1 (MQ=255)
                                                                        gCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTCCGCCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGAt                              <  1:3900795/101‑1 (MQ=255)
                                                                         cGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATg                             <  2:1281722/101‑1 (MQ=255)
                                                                          gATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                            >  1:1021239/1‑101 (MQ=255)
                                                                           atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                           <  2:3493566/101‑1 (MQ=255)
                                                                           agatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                           >  2:3073693/3‑101 (MQ=2)
                                                                           agatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                           >  1:1514285/3‑101 (MQ=2)
                                                                            tatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGTTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCa                          <  2:509437/101‑1 (MQ=255)
                                                                            tatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCa                          <  2:2651056/101‑1 (MQ=255)
                                                                             atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                         <  2:1006866/101‑1 (MQ=255)
                                                                              tcggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                        <  1:2026796/98‑2 (MQ=9)
                                                                              tcggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                        <  1:4662516/98‑2 (MQ=9)
                                                                              tCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAtt                        >  2:3743081/1‑101 (MQ=255)
                                                                               ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGTTGAGCGCGACCATGACCTTc                       >  2:1778357/1‑101 (MQ=25)
                                                                               ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTc                       >  2:4705365/1‑101 (MQ=255)
                                                                                 gcgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGt                     <  1:1946983/101‑1 (MQ=255)
                                                                                 gcgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGt                     <  2:1755759/101‑1 (MQ=255)
                                                                                   ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                   >  2:3540454/1‑101 (MQ=255)
                                                                                   ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                   <  1:2538563/101‑1 (MQ=255)
                                                                                    gTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGcc                  >  1:3410460/1‑101 (MQ=255)
                                                                                     tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                 >  2:2761975/1‑101 (MQ=255)
                                                                                      ttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCa                <  2:3401371/101‑1 (MQ=255)
                                                                                       tCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGGGTCGCCCTCAGCCGTGCCGCCTGCGCCGTGGTAAGCGCCACGATGTCATTCGTGCCCac               >  2:793276/1‑100 (MQ=255)
                                                                                       tCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGGCACGATGTCATTCGTGCCCAg               >  2:4554679/1‑101 (MQ=255)
                                                                                       tCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAg               >  1:4653867/1‑101 (MQ=255)
                                                                                         cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcc             >  1:1157146/1‑101 (MQ=255)
                                                                                         cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcc             <  1:4333634/101‑1 (MQ=255)
                                                                                          aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc            <  1:3615588/101‑1 (MQ=255)
                                                                                          aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc            >  2:1567097/1‑101 (MQ=255)
                                                                                          aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc            >  2:2848389/1‑101 (MQ=255)
                                                                                            cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt          <  1:2075605/101‑1 (MQ=255)
                                                                                            cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt          <  1:2904491/101‑1 (MQ=255)
                                                                                            cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt          <  2:288489/101‑1 (MQ=255)
                                                                                             gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt         <  1:86219/101‑1 (MQ=255)
                                                                                             gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt         >  1:811738/1‑101 (MQ=255)
                                                                                             gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt         <  2:2014442/101‑1 (MQ=255)
                                                                                               cgcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCa       <  2:1021160/101‑1 (MQ=255)
                                                                                               cgcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCa       <  1:1481648/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CGGATCTGCGCCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCC  >  NC_014034/678673‑678873

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