Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 682,593 C→G 7.8% E338D (GAG→GAC RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034682,5930C→G 7.8% 27.9 ‑0.2 ‑0.1193E338D (GAG→GACRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (6/9):  ref base (80/97):  total (87/106)
Fisher's exact test strand distribution p-value = 7.91e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.70e-01

Only 100 of 204 total aligned reads displayed.

TCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGCC  >  NC_014034/682493‑682692
                                                                                                     |                                                                                                   
  tGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.GGCTCAGCCCGTCCATGTCCTc                                                                                                   >  1:1589040/1‑101 (MQ=255)
   gCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCg                                                                                                  <  2:936180/101‑1 (MQ=255)
     cgcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                                                <  1:1128282/101‑1 (MQ=255)
     cgcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGt                                                                                                <  2:391157/101‑1 (MQ=255)
      gcgTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTc                                                                                               <  2:1432816/101‑1 (MQ=255)
        gTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGa                                                                                             <  1:3336628/101‑1 (MQ=255)
         tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                                                            <  2:184440/101‑1 (MQ=255)
         tGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaa                                                                                            >  1:1021938/1‑101 (MQ=255)
          gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                                           <  1:2765948/101‑1 (MQ=255)
          gCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGaaa                                                                                           >  2:473556/1‑101 (MQ=255)
           cccAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAg                                                                                          <  1:2615838/101‑1 (MQ=255)
              aGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGccc                                                                                       <  1:32098/101‑1 (MQ=255)
               gAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGcccc                                                                                      <  2:2071888/101‑1 (MQ=255)
                aaCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCt                                                                                     <  2:3633457/101‑1 (MQ=255)
                 aCCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCtt                                                                                    <  1:2639510/101‑1 (MQ=255)
                   cgccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCttga                                                                                  >  1:3597997/3‑99 (MQ=37)
                   ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                                                  <  2:4902726/101‑1 (MQ=255)
                   ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                                                  <  2:4642083/101‑1 (MQ=255)
                   ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                                                  <  2:1320517/101‑1 (MQ=255)
                   ccccAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCa                                                                                  <  2:1694016/101‑1 (MQ=255)
                     ccAGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgc                                                                                <  1:1611283/101‑1 (MQ=25)
                     ccACCTGGGCCGTGCTGCGCCCCGCCATCTGCGCCCCGGGCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgc                                                                                <  1:1421511/101‑1 (MQ=255)
                       aGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAgcgc                                                                              >  2:1011348/1‑101 (MQ=255)
                       aGCTGCGCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAgcgc                                                                              <  2:2466930/101‑1 (MQ=25)
                        gCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGcgcg                                                                             >  2:2495135/1‑101 (MQ=255)
                          tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                                           >  2:4221088/1‑101 (MQ=255)
                          tGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggc                                                                           >  2:4998914/1‑101 (MQ=255)
                           gggCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCgcggcg                                                                          <  2:738270/101‑1 (MQ=255)
                             gCCGTGCTCAGCGCCGCCACCTGCTGCGAGGTCATCGCCGCCA.TCTGATCGCCGGTCAGCGCGACGATGTCGCCGGTCGACAGCCCCTTGAGcgcgccact                                                                        <  2:4358505/101‑6 (MQ=21)
                              ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                                                       <  2:2592166/101‑1 (MQ=255)
                              ccGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTg                                                                       >  2:2823665/1‑101 (MQ=255)
                                gTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGcc                                                                     >  2:4419686/1‑101 (MQ=255)
                                  gCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGa                                                                   >  2:2326604/1‑101 (MQ=255)
                                   cTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGAt                                                                  >  2:4218573/1‑101 (MQ=255)
                                    tCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                                                                 <  1:1774271/101‑1 (MQ=255)
                                    tCAGCGCCGCCATCTGCCCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATc                                                                 <  2:1588908/101‑1 (MQ=255)
                                     cAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcg                                                                >  2:3056828/1‑101 (MQ=255)
                                       gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgcc                                                              <  2:4579640/101‑1 (MQ=255)
                                       gcgcCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCACCCCGGCCATGTCCTCGGTCGAAAGCCCGTTCAGCGCGGCGGTGCCGATcccc                                                              >  2:297821/1‑101 (MQ=255)
                                        cgccgcCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATcgccg                                                             <  1:2660942/101‑1 (MQ=255)
                                          ccggcATCCGCGCCCCGGCCCGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                                           <  1:3564818/101‑1 (MQ=255)
                                          ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                                           <  1:716736/101‑1 (MQ=255)
                                          ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                                           <  1:3424451/101‑1 (MQ=255)
                                          ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                                           <  1:2163734/101‑1 (MQ=255)
                                          ccgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCa                                                           >  1:4183555/1‑101 (MQ=255)
                                           cgccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaa                                                          >  2:1482470/1‑101 (MQ=255)
                                            gccATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCaaa                                                         >  1:3684156/1‑101 (MQ=255)
                                                tCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGcc                                                     <  2:4470436/101‑1 (MQ=255)
                                                 cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                                                    >  1:4852117/1‑101 (MQ=255)
                                                 cTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCcg                                                    >  1:3456272/1‑101 (MQ=255)
                                                      ccccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCg                                               <  1:475025/101‑1 (MQ=255)
                                                        ccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGa                                             >  1:4354544/1‑101 (MQ=255)
                                                        ccGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGa                                             <  2:543005/101‑1 (MQ=255)
                                                         cGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaa                                            <  1:629904/101‑1 (MQ=255)
                                                          ggTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGAGCCGGaaa                                           >  2:832547/1‑101 (MQ=255)
                                                           gTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGaaaa                                          <  1:897490/101‑1 (MQ=255)
                                                             cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                                        >  2:3138130/1‑101 (MQ=255)
                                                             cAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAgg                                        <  1:453346/101‑1 (MQ=255)
                                                              aGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGgc                                       <  2:3403013/101‑1 (MQ=255)
                                                              aGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGgc                                       <  1:2053398/101‑1 (MQ=255)
                                                                 gccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGcc                                    <  1:5024800/101‑1 (MQ=255)
                                                                  ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                   <  1:305085/101‑1 (MQ=255)
                                                                  ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                   <  2:1133609/101‑1 (MQ=255)
                                                                  ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                   <  1:807958/101‑1 (MQ=255)
                                                                  ccgccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCg                                   <  1:4420005/101‑1 (MQ=255)
                                                                    gccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGAt                                 <  1:3738794/101‑1 (MQ=255)
                                                                     ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATc                                <  2:3049152/101‑1 (MQ=255)
                                                                     ccAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCCCCCGGGAAAAGGCCCCgttt                                >  2:3761368/1‑98 (MQ=255)
                                                                      cAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt                               >  2:4155948/1‑101 (MQ=255)
                                                                      cAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCt                               <  1:2419120/101‑1 (MQ=255)
                                                                       aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCGGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg                              <  1:3991324/101‑1 (MQ=255)
                                                                       aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg                              <  1:4648690/101‑1 (MQ=255)
                                                                       aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg                              <  1:476802/101‑1 (MQ=255)
                                                                       aGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTg                              >  2:3861406/1‑101 (MQ=255)
                                                                         ttCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCt                            <  1:2924968/101‑1 (MQ=255)
                                                                          tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc                           <  2:3290996/101‑1 (MQ=255)
                                                                          tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc                           <  1:536372/101‑1 (MQ=255)
                                                                          tCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTc                           <  1:3861691/101‑1 (MQ=255)
                                                                            gATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCgg                         >  1:1885267/1‑101 (MQ=255)
                                                                             aTCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGc                        >  2:4859325/1‑101 (MQ=255)
                                                                               cG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCg                      <  2:713029/101‑1 (MQ=255)
                                                                                g.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCga                     <  2:4530313/101‑1 (MQ=255)
                                                                                g.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCga                     <  1:1573805/101‑1 (MQ=255)
                                                                                  tGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCgag                    <  2:1343922/101‑1 (MQ=255)
                                                                                   gCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCgaga                   <  1:2495309/101‑1 (MQ=255)
                                                                                   gCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCgaga                   >  1:4946052/1‑101 (MQ=255)
                                                                                   gCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCCATCTGCTCGGCCgagg                   >  1:3092329/1‑100 (MQ=255)
                                                                                     tCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAgc                 <  1:4999244/101‑1 (MQ=255)
                                                                                      cAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGCGAgcg                >  1:298008/1‑101 (MQ=255)
                                                                                        gCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGcc              >  1:1652843/1‑101 (MQ=255)
                                                                                           cGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGaa           >  1:2561767/1‑101 (MQ=255)
                                                                                             gCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAg         >  1:4280127/1‑101 (MQ=255)
                                                                                               cATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGaaaccc       >  2:3158432/1‑98 (MQ=255)
                                                                                                aTGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTg      >  1:19405/1‑101 (MQ=255)
                                                                                                aTGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTg      >  2:1073053/1‑101 (MQ=255)
                                                                                                  gTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTgcg    >  1:781919/1‑101 (MQ=255)
                                                                                                   tCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTgcgc   >  1:2370445/1‑101 (MQ=255)
                                                                                                    ccTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGcc  <  1:2232838/101‑1 (MQ=255)
                                                                                                    ccTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGcc  <  1:2823878/101‑1 (MQ=255)
                                                                                                    ccTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGcc  >  2:3809763/1‑101 (MQ=255)
                                                                                                     |                                                                                                   
TCTGCCGCGTGCCCAGAACCCCCAGCTGGGCCGTGCTCAGCGCCGCCATCTGCGCCCCGGTCAGCGCCGCCAGTTCGATCG.TGCTCAGCCCGGCCATGTCCTCGGTCGAAAGCCCCTTCAGCGCGGCGGTGCCGATCGCCGCAAAGGCCGCGCCGGAAAAGGCGCCGATCTGCTCGGCCGAGAGCGCCGAAAGCTGCGCC  >  NC_014034/682493‑682692

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