Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,459 G→A 20.7% I716I (ATC→ATT RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,4590G→A20.7% 107.1 ‑2.2 ‑0.1193I716I (ATC→ATTRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (15/25):  ref base (81/72):  total (96/97)
Fisher's exact test strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.41e-01

Only 100 of 215 total aligned reads displayed.

CCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCG  >  NC_014034/681359‑681559
                                                                                                    |                                                                                                      
 cgcgATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGa                                                                                                       >  2:4309452/1‑101 (MQ=255)
     aTCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTg                                                                                                   >  2:1129468/1‑101 (MQ=255)
      tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                                  <  2:2445812/101‑1 (MQ=255)
      tCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGc                                                                                                  <  1:4519515/101‑1 (MQ=255)
       cTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCg                                                                                                 >  1:514676/1‑101 (MQ=255)
        tGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGt                                                                                                <  1:1656132/101‑1 (MQ=255)
         gcgcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTg                                                                                               >  2:2928914/1‑101 (MQ=255)
          cgcgcCTGAAGGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCCGATCCGCCGCCTCGAGCGGCGCGAGCTGCGTGc                                                                                              >  1:4363226/1‑101 (MQ=255)
          cgcgcCTGAAAGCCCAGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGc                                                                                              >  1:2396552/1‑101 (MQ=255)
           gcgcCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCt                                                                                             <  1:2843335/101‑1 (MQ=255)
               cTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTc                                                                                         <  2:2032147/101‑1 (MQ=255)
                tGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCa                                                                                        <  1:4208096/101‑1 (MQ=255)
                 gAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTTCGTGCTGGTCAg                                                                                       >  1:2027285/1‑101 (MQ=255)
                    aGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAgcgc                                                                                    <  2:3048653/101‑1 (MQ=255)
                       ccGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGNGTGCTGGTCAGCGCGcc                                                                                 >  2:1838227/1‑101 (MQ=255)
                        cGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCg                                                                                >  2:1532296/1‑101 (MQ=255)
                          gCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.At                                                                             <  1:690547/101‑1 (MQ=255)
                           ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATc                                                                            <  1:1387547/101‑1 (MQ=255)
                           ccATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATc                                                                            >  1:3072881/1‑101 (MQ=255)
                            cATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCt                                                                           >  1:3875933/1‑101 (MQ=255)
                             aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGCGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                          >  1:4751252/1‑101 (MQ=255)
                             aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                          >  2:5014923/1‑101 (MQ=255)
                             aTCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTg                                                                          >  1:3997497/1‑101 (MQ=255)
                              tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgc                                                                         >  2:2026532/1‑101 (MQ=255)
                              tCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTgc                                                                         <  2:2505273/101‑1 (MQ=255)
                                 ttgCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGTGCGCCGCGATCTGCCTGCTGGTGGGCGCGCCG.ATCTgcgcg                                                                      >  2:932109/3‑100 (MQ=38)
                                  ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGGGCCG.ATCTGCGCCg                                                                     >  2:2361997/1‑101 (MQ=255)
                                  ggCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCg                                                                     >  1:4486761/1‑101 (MQ=255)
                                   gCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGa                                                                    >  1:3513976/1‑101 (MQ=255)
                                    ccGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAg                                                                   <  2:3735592/101‑1 (MQ=255)
                                     cGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGc                                                                  >  2:866258/1‑101 (MQ=255)
                                      gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGcc                                                                 >  2:785457/1‑101 (MQ=255)
                                      gAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGcc                                                                 <  1:547224/101‑1 (MQ=255)
                                        gcccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCa                                                               <  2:3982639/101‑1 (MQ=255)
                                         cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAg                                                              <  2:4143183/101‑1 (MQ=255)
                                         cccagcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAg                                                              >  1:3356525/1‑101 (MQ=255)
                                             gcccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcgc                                                          <  1:3949258/101‑1 (MQ=255)
                                              cccGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGAcgc.g                                                        >  1:4137795/1‑101 (MQ=255)
                                                cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGc.gcc                                                      <  2:4839140/101‑1 (MQ=255)
                                                cGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGc.gcc                                                      <  1:4912142/101‑1 (MQ=255)
                                                 ggAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGc.gccg                                                     >  1:127957/1‑101 (MQ=255)
                                                   agagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.gccgcc                                                   <  1:1019758/101‑1 (MQ=255)
                                                    gagCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCt                                                  <  2:341714/101‑1 (MQ=255)
                                                       cgcCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTgcg                                               >  1:1941853/1‑101 (MQ=255)
                                                        gcCGTCGCGCTGATCTGCGCCAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTgcgc                                              >  1:4490819/1‑101 (MQ=255)
                                                          cgccgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATGTGCGTGCTGGTCAGGGCGCCG.ATCTGCGCCGAGCCAAGACGC.GCCGCCTGCGCcg                                            <  2:619325/98‑1 (MQ=25)
                                                          cGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcg                                            <  2:1850689/101‑1 (MQ=255)
                                                           gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgt                                           >  1:3677838/1‑101 (MQ=255)
                                                           gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgt                                           >  2:48974/1‑101 (MQ=255)
                                                           gTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCATACGCcgt                                           >  2:4007452/1‑101 (MQ=255)
                                                            gcgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtc                                          <  2:1452435/100‑1 (MQ=255)
                                                             cgcgCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcg                                         >  2:4795266/1‑101 (MQ=255)
                                                              gcgcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCcgtcgt                                        <  2:2917503/101‑1 (MQ=255)
                                                                gcTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCa                                      <  1:3279146/101‑1 (MQ=255)
                                                                 cTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGACTGCGCCGTCGCCAg                                     >  1:4445062/1‑101 (MQ=255)
                                                                  tGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GGCGCCTGCGCCGTCGTCAgc                                    >  1:2964804/1‑101 (MQ=255)
                                                                  ggagccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAgc                                    <  1:1093843/92‑1 (MQ=2)
                                                                   gATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAgcg                                   >  2:4562526/1‑101 (MQ=255)
                                                                      ccgcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCa                                <  1:5008976/96‑1 (MQ=2)
                                                                        gcgcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACg                              >  2:3791042/1‑101 (MQ=255)
                                                                        gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCTTGCTCGACAGGGCCAGA.AGCTGCGTCGAGCCCAG.CGCGGCCGCCTGCGCCGTCGTCAGCCCCACg                              >  2:1636690/4‑101 (MQ=2)
                                                                        gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTGGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg                              <  2:3281452/98‑1 (MQ=2)
                                                                        gcaccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACg                              <  2:2598648/98‑1 (MQ=2)
                                                                         cgcCAGATCCGCCGCCTCGAGGGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGa                             >  1:3717894/1‑101 (MQ=255)
                                                                         caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGa                             <  2:1836725/99‑1 (MQ=9)
                                                                         caccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGACGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGGGCCACGa                             >  2:5021035/3‑101 (MQ=9)
                                                                          gcCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGAt                            <  1:18655/101‑1 (MQ=255)
                                                                           ccAGATCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                           <  2:2788658/101‑1 (MQ=9)
                                                                           ccAGATCCGCCGGCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                           >  1:2738498/1‑101 (MQ=255)
                                                                           ccAGATCCGCCGCCTGGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATg                           <  1:889882/101‑1 (MQ=255)
                                                                              gATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                        >  2:4499916/1‑101 (MQ=255)
                                                                               aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                       <  2:3819267/101‑1 (MQ=255)
                                                                               aTCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAt                       >  1:3870981/1‑101 (MQ=255)
                                                                                tCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATc                      <  1:4607355/101‑1 (MQ=255)
                                                                                   gcggtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGt                   >  2:550453/10‑101 (MQ=9)
                                                                                      gtttccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                <  2:2089590/95‑1 (MQ=12)
                                                                                      gccTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGcc                >  2:173597/1‑101 (MQ=255)
                                                                                         tccaGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgag             <  1:3852182/98‑4 (MQ=9)
                                                                                         tCGAGCGCCGCGCTCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa             >  2:3519808/1‑101 (MQ=255)
                                                                                         tCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCaaa             <  1:1532426/101‑1 (MQ=255)
                                                                                          cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg            >  1:852211/1‑101 (MQ=255)
                                                                                          cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg            <  2:4749914/101‑1 (MQ=255)
                                                                                          cGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAg            <  2:3875641/101‑1 (MQ=255)
                                                                                           gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc           >  2:4187818/1‑101 (MQ=255)
                                                                                           gAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGc           <  1:2579101/101‑1 (MQ=255)
                                                                                           caGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGccgaggc           <  2:1892595/100‑6 (MQ=9)
                                                                                            aGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGCCATCGGTGCCAAAGcc          >  1:1315237/1‑101 (MQ=255)
                                                                                            aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc          >  2:3152601/1‑101 (MQ=9)
                                                                                            aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc          <  2:4721859/101‑1 (MQ=9)
                                                                                            aGCGCCGCAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc          <  1:201605/101‑1 (MQ=9)
                                                                                            aGCGCCGCAACCTGCGTGCTCGACAGCGC.CAGAAGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGcc          >  1:2113145/1‑101 (MQ=2)
                                                                                              cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct        >  1:4160371/1‑99 (MQ=2)
                                                                                              cgccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccct        >  2:2498584/1‑99 (MQ=2)
                                                                                               gccgcAACCTGCGTGCTCGACAGCGCCAGA.AGCTGCGTCGAGCCCAGCGCC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGccctt       <  2:1268038/101‑4 (MQ=9)
                                                                                                 cgcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGa     <  1:229757/101‑1 (MQ=255)
                                                                                                  tcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt    <  2:3618756/100‑1 (MQ=255)
                                                                                                  gcgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGAt    <  1:3971375/101‑1 (MQ=255)
                                                                                                   cgATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATc   >  2:1597303/1‑101 (MQ=255)
                                                                                                    gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg  <  2:514597/101‑1 (MQ=255)
                                                                                                    gATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcg  <  2:3072643/101‑1 (MQ=255)
                                                                                                    |                                                                                                      
CCGCGATCTGCGCGCCTGAAAGCCCGGCCATCGTGGCCGAGCCCAGCCCGGAGAGCGCCGTCGCGCTGATCTGCGCCAGATCCGCCGCCTCGAGCGCCGCGATCTGCGTGCTGGTCAGCGCGCCG.ATCTGCGCCGAGCCCAGACGC.GCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCG  >  NC_014034/681359‑681559

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG