Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 678,783 G→C 8.4% S1608R (AGC→ACC RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034678,7830G→C 8.4% 31.5 ‑1.9 ‑0.1231S1608R (AGC→AGGRCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (8/12):  ref base (124/87):  total (132/99)
Fisher's exact test strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.65e-01

Only 100 of 246 total aligned reads displayed.

CCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCA  >  NC_014034/678683‑678883
                                                                                                    |                                                                                                    
ccGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTg                                                                                                      >  1:4434075/1‑101 (MQ=255)
ccGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTg                                                                                                      <  1:2765585/101‑1 (MQ=255)
 cGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGc                                                                                                     <  1:3327522/101‑1 (MQ=255)
  gTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCt                                                                                                    >  2:4380430/1‑101 (MQ=255)
   tGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTc                                                                                                   <  1:1944334/101‑1 (MQ=255)
    gCTCAGGCTTGGCATCAGCGCTGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCa                                                                                                  >  1:2366120/1‑101 (MQ=255)
    gCTCAGGCTTGGCATCAGCGCCATCCAGATCTCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCTGCACCGCAACCTACGCCGTCCTCa                                                                                                  >  2:4383677/1‑101 (MQ=255)
     cTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAg                                                                                                 >  1:1004151/1‑101 (MQ=255)
      tCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgc                                                                                                >  2:2069986/1‑101 (MQ=255)
       cAGGCTTGGCATCGGCGCCGTCCCGACCCCCGCCATGGCCGAGGCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                                               <  1:1819664/101‑1 (MQ=255)
       cAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                                                                               >  2:3917082/1‑101 (MQ=255)
          gCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccg                                                                                            <  1:3995561/101‑1 (MQ=255)
           cTTGGCATCAGCGCCGTCCCGGTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATCTCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                                           <  2:2308292/101‑1 (MQ=255)
           cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                                           <  1:3828410/101‑1 (MQ=255)
           cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                                           >  2:1101326/1‑101 (MQ=255)
           cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                                           >  2:3014909/1‑101 (MQ=255)
           cTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGAGCCGCGATCTGCGCCGTGCTCAGcgccgc                                                                                           >  2:4740848/1‑101 (MQ=255)
            ttGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCcgcg                                                                                          >  2:2416628/1‑101 (MQ=255)
             tGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGa                                                                                         >  1:3289169/1‑101 (MQ=255)
              ggCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGAt                                                                                        >  1:4973590/1‑101 (MQ=255)
              ggCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGAt                                                                                        <  2:636654/101‑1 (MQ=255)
                cATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCt                                                                                      <  1:4752942/101‑1 (MQ=255)
                 aTCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTg                                                                                     >  2:1073485/1‑101 (MQ=255)
                 aTCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTg                                                                                     >  2:1466747/1‑101 (MQ=255)
                 aTCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTg                                                                                     >  2:3995410/1‑101 (MQ=255)
                  tCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGc                                                                                    <  1:4389165/101‑1 (MQ=255)
                   cAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCg                                                                                   >  1:2559110/1‑101 (MQ=255)
                   cAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCg                                                                                   <  1:4294142/101‑1 (MQ=255)
                    aGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGt                                                                                  >  1:1123827/1‑101 (MQ=255)
                    aGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGt                                                                                  >  1:1632562/1‑101 (MQ=255)
                     gcgcCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTc                                                                                 <  2:4173205/101‑1 (MQ=255)
                       gcCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTcgc                                                                               >  1:2893193/1‑101 (MQ=255)
                           tCCCGCTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTc                                                                           <  1:1018428/101‑1 (MQ=255)
                               gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                                       >  2:1405623/1‑101 (MQ=255)
                               gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                                       <  1:2828708/101‑1 (MQ=255)
                               gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                                       >  1:2427986/1‑101 (MQ=255)
                                aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                      >  2:3305932/1‑101 (MQ=255)
                                aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                      >  2:2254459/1‑101 (MQ=255)
                                aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                      >  1:303726/1‑101 (MQ=255)
                                aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                      >  1:1073056/1‑101 (MQ=255)
                                   cccGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgc                                                                   <  2:4669758/101‑1 (MQ=255)
                                       ccATTGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccgcc                                                               >  1:607837/1‑101 (MQ=255)
                                        cATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCt                                                              >  1:4186279/1‑101 (MQ=255)
                                           ggCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgcg                                                           >  1:3505347/1‑101 (MQ=255)
                                              cGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCg                                                        >  2:677546/1‑101 (MQ=255)
                                               gagaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGt                                                       >  2:1743963/1‑101 (MQ=255)
                                                  aCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTCAGCCGTggc                                                    >  2:1237982/1‑100 (MQ=255)
                                                   ccGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTc                                                   <  1:571765/101‑1 (MQ=255)
                                                    cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                                  >  1:3709515/1‑101 (MQ=255)
                                                    cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                                  >  2:1242778/1‑101 (MQ=255)
                                                          cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCa                                            >  2:1831965/1‑101 (MQ=255)
                                                          cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCa                                            >  2:278468/1‑101 (MQ=255)
                                                          cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCa                                            >  2:4930940/1‑101 (MQ=255)
                                                           tGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCAc                                           >  2:4019432/1‑101 (MQ=255)
                                                            gggCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACg                                          >  1:2779911/1‑101 (MQ=255)
                                                             taccagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGa                                         >  2:3763456/7‑101 (MQ=2)
                                                              accagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGAt                                        >  1:2546915/6‑101 (MQ=2)
                                                               cGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATg                                       <  2:4083351/101‑1 (MQ=255)
                                                               cGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATg                                       <  2:1855715/101‑1 (MQ=255)
                                                                gATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                                      >  2:3003578/1‑101 (MQ=255)
                                                                 atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGGCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                     >  1:343626/1‑101 (MQ=255)
                                                                 atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                     <  1:4931270/101‑1 (MQ=255)
                                                                 atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                     <  1:1405750/101‑1 (MQ=255)
                                                                 atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                     >  1:1339708/1‑101 (MQ=255)
                                                                 agatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                     <  2:36773/99‑1 (MQ=2)
                                                                  gatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                    <  2:2969633/100‑1 (MQ=2)
                                                                  gatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                    >  1:4624860/2‑101 (MQ=2)
                                                                   atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                                   <  1:1101440/101‑1 (MQ=255)
                                                                   atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                                   <  2:1632430/101‑1 (MQ=255)
                                                                   atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                                   <  1:2075483/101‑1 (MQ=255)
                                                                   atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                                   >  1:2096005/1‑101 (MQ=255)
                                                                    tCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAtt                                  <  2:4603452/101‑1 (MQ=255)
                                                                     ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTc                                 <  2:2365956/101‑1 (MQ=255)
                                                                     ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTc                                 <  1:535249/101‑1 (MQ=255)
                                                                         ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                             >  1:5016304/1‑101 (MQ=255)
                                                                          gTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGcc                            >  1:679118/1‑101 (MQ=255)
                                                                           tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                           <  1:1466877/101‑1 (MQ=255)
                                                                           tttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGccc                           <  2:1004041/101‑1 (MQ=255)
                                                                            ttCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCa                          <  1:220645/101‑1 (MQ=255)
                                                                             tCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAg                         >  1:316009/1‑101 (MQ=255)
                                                                              ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGc                        >  2:3519762/1‑101 (MQ=255)
                                                                              ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGc                        >  1:3211485/1‑101 (MQ=255)
                                                                               cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcc                       <  2:257146/101‑1 (MQ=255)
                                                                                aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc                      <  2:2095857/101‑1 (MQ=255)
                                                                                aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCCGCGCCGTGGCCAGCGCCACGATGTCATTCGTGCCCAGccc                      >  2:841827/1‑101 (MQ=255)
                                                                                  cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                    <  2:524066/101‑1 (MQ=255)
                                                                                    ccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTc                  <  1:1512969/101‑1 (MQ=255)
                                                                                     cgcgATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCa                 <  2:4431951/101‑1 (MQ=255)
                                                                                         aTCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGCTGTCATTCGTGCCCAGCCCCTTCAgcgc             >  1:2287467/1‑101 (MQ=255)
                                                                                         aTCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAgcgc             >  2:3348809/1‑101 (MQ=255)
                                                                                           cTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGcgccg           <  2:4965254/101‑1 (MQ=255)
                                                                                            tGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGcgccgc          >  2:782685/1‑101 (MQ=255)
                                                                                             gcgcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCACTTCAGCGCcgcg         >  2:4463504/1‑101 (MQ=255)
                                                                                              cgcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCgg        >  2:3893792/1‑101 (MQ=255)
                                                                                              cgcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCgg        <  2:553929/101‑1 (MQ=255)
                                                                                              cgcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCgg        >  2:3189083/1‑101 (MQ=255)
                                                                                               gcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGt       <  1:3726844/101‑1 (MQ=255)
                                                                                                   tGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTc   >  1:3956932/1‑101 (MQ=255)
                                                                                                   tGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTc   <  1:586262/101‑1 (MQ=255)
                                                                                                    gCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCa  <  2:218249/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CCGTGCTCAGGCTTGGCATCAGCGCCGTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCA  >  NC_014034/678683‑678883

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