Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 678,809 T→C 6.6% S1600G (AGC→GGC)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034678,8090T→C 6.6% 24.3 ‑3.0 ‑0.2214S1600G (AGC→GGC) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/11):  ref base (100/99):  total (104/110)
Fisher's exact test strand distribution p-value = 5.12e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.45e-01

Only 100 of 224 total aligned reads displayed.

GTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTCGACAGC  >  NC_014034/678709‑678909
                                                                                                    |                                                                                                    
gTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCt                                                                                                      <  2:3231938/101‑1 (MQ=255)
  cccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                                                    <  2:2548102/101‑1 (MQ=255)
  cccGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                                                    <  1:2189868/101‑1 (MQ=255)
  cccGATCCCCGCCAGGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTGTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                                                    <  2:3142842/101‑1 (MQ=255)
  cccGAGCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCa                                                                                                    <  2:1429690/101‑1 (MQ=255)
    cGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGc                                                                                                  <  1:3375857/101‑1 (MQ=255)
     gATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGcc                                                                                                 <  2:3533410/101‑1 (MQ=255)
      aTCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCg                                                                                                >  2:3838216/1‑101 (MQ=255)
            gCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccgc                                                                                          >  2:1245864/1‑101 (MQ=255)
             ccATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTgccgcc                                                                                         >  2:567863/1‑101 (MQ=255)
             ccATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCACTCAGCCGTgccgcc                                                                                         <  1:2607950/101‑1 (MQ=255)
              cATGGCCGAGCCCGAAAGCTGGGCGATATCCGCTGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCt                                                                                        <  2:748247/101‑1 (MQ=255)
              cATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCt                                                                                        <  2:2813295/101‑1 (MQ=255)
               aTGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTg                                                                                       <  1:174461/101‑1 (MQ=255)
                tGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTgc                                                                                      >  1:1384655/1‑101 (MQ=255)
                    cGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCg                                                                                  >  2:4596139/1‑101 (MQ=255)
                      agaCCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTg                                                                                <  1:1401538/101‑1 (MQ=255)
                       gaCCGAAAGCTGGGCGGTATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCCCTCAGCCGTCCCGCCTGCGCCGTgg                                                                               <  1:776863/101‑1 (MQ=255)
                         ccGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTc                                                                             <  2:4163993/101‑1 (MQ=255)
                          cGAAATCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                                                            <  2:424041/101‑1 (MQ=255)
                          cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                                                            <  1:4543010/101‑1 (MQ=255)
                          cGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCa                                                                            <  1:94752/101‑1 (MQ=255)
                             aaGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgcg                                                                         >  1:3684962/1‑101 (MQ=255)
                              aGCTTGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAgcgc                                                                        >  1:235204/1‑101 (MQ=255)
                                cTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCa                                                                      <  2:775255/101‑1 (MQ=255)
                                 tGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCAc                                                                     >  1:4065907/1‑101 (MQ=255)
                                 tGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCGCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCAc                                                                     >  1:4047868/1‑101 (MQ=255)
                                  gggCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACg                                                                    >  1:2638954/1‑101 (MQ=255)
                                  gggCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACg                                                                    <  1:1234920/101‑1 (MQ=255)
                                  gggCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACg                                                                    <  1:822837/101‑1 (MQ=255)
                                    gCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGAt                                                                  <  1:4735033/101‑1 (MQ=255)
                                      gATTTCCGCGGTTTCCAGCGCCGCGATCGCCGCCGTCCTCAGCCCCGCGATCTGCGTCACGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                                                                <  2:1102812/101‑1 (MQ=25)
                                      gATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGt                                                                >  1:1021239/1‑101 (MQ=255)
                                      cagatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGt                                                                >  2:3403579/4‑101 (MQ=2)
                                       atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                                               <  2:3493566/101‑1 (MQ=255)
                                       atatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTc                                                               <  1:2102398/101‑1 (MQ=255)
                                       agatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTc                                                               >  2:3073693/3‑101 (MQ=2)
                                        tatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCa                                                              <  2:2651056/101‑1 (MQ=255)
                                        tatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCa                                                              >  2:4720270/1‑101 (MQ=255)
                                        tatCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCa                                                              >  1:2531505/1‑101 (MQ=255)
                                        gatCGGCGGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCa                                                              <  1:3468438/100‑1 (MQ=2)
                                         atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                                                             >  2:4546108/1‑101 (MQ=255)
                                         atCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCAt                                                             <  1:3838451/101‑1 (MQ=255)
                                          tcggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                                                            <  1:2026796/98‑2 (MQ=9)
                                          tcggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                                                            <  1:4662516/98‑2 (MQ=9)
                                          tcggcgGTTTCCAGCGCCGCAACCTGCGTGCTCGACAGCGCCAGAAGCTGCGTCGAGCCCAGCGCCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCAtc                                                            <  1:3275503/98‑2 (MQ=9)
                                           ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGTTGAGCGCGACCATGACCTTc                                                           >  2:1778357/1‑101 (MQ=25)
                                           ccGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTc                                                           >  1:38866/1‑101 (MQ=255)
                                            cgcgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCg                                                          >  1:1435986/1‑101 (MQ=255)
                                             gcgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGt                                                         <  1:1946983/101‑1 (MQ=255)
                                              cgGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACCATGTCATTCGTg                                                        <  1:2565181/101‑1 (MQ=255)
                                               ggTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGc                                                       >  2:2917815/1‑101 (MQ=255)
                                                   tCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGGGTCGCCCTCAGCCGTGCCGCCTGCGCCGTGGTAAGCGCCACGATGTCATTCGTGCCCac                                                   >  2:793276/1‑100 (MQ=255)
                                                   tCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAg                                                   >  1:4653867/1‑101 (MQ=255)
                                                    ccAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGc                                                  >  1:1974822/1‑101 (MQ=255)
                                                     cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcc                                                 <  1:4333634/101‑1 (MQ=255)
                                                     cAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGcc                                                 >  1:1157146/1‑101 (MQ=255)
                                                      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc                                                <  1:3615588/101‑1 (MQ=255)
                                                      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc                                                <  2:1012580/101‑1 (MQ=255)
                                                      aGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGccc                                                >  2:1567097/1‑101 (MQ=255)
                                                        cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                                              >  2:86203/1‑101 (MQ=255)
                                                        cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                                              <  2:288489/101‑1 (MQ=255)
                                                        cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                                              >  2:647786/1‑101 (MQ=255)
                                                        cgccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCt                                              <  1:2904491/101‑1 (MQ=255)
                                                         gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt                                             <  1:3088379/101‑1 (MQ=255)
                                                         gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt                                             <  2:2265594/101‑1 (MQ=255)
                                                         gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt                                             <  2:2014442/101‑1 (MQ=255)
                                                         gccgcGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCtt                                             >  2:1750493/1‑101 (MQ=255)
                                                              gATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAccg                                        >  1:1509507/1‑101 (MQ=255)
                                                               aTCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAgcgc                                       >  2:2752963/1‑101 (MQ=255)
                                                                tCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGcgcc                                      <  2:38859/101‑1 (MQ=255)
                                                                   gcgcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCcgcg                                   >  1:1598106/1‑101 (MQ=255)
                                                                     gcCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGt                                 <  2:4348022/101‑1 (MQ=255)
                                                                      ccGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTg                                <  1:633917/101‑1 (MQ=255)
                                                                        gTGCTCAGCGCCGCGATCTGCGTCGCGCTCGGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCt                              >  2:1945794/1‑101 (MQ=255)
                                                                         tGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTc                             <  1:4677350/101‑1 (MQ=255)
                                                                          gCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCa                            <  1:2051173/101‑1 (MQ=255)
                                                                           cTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAg                           <  1:2445970/101‑1 (MQ=255)
                                                                           cTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAg                           >  2:4324259/1‑101 (MQ=255)
                                                                             cAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgcg                         >  2:4568070/1‑101 (MQ=255)
                                                                             cAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgcg                         >  1:799631/1‑101 (MQ=255)
                                                                              aGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgcgc                        >  1:228242/1‑101 (MQ=255)
                                                                              aGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgcgc                        <  1:1373598/101‑1 (MQ=255)
                                                                              aGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgcgc                        >  2:3350234/1‑101 (MQ=255)
                                                                                cgccgcGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAgcgcgc                      >  1:4150207/1‑101 (MQ=255)
                                                                                  ccgcGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCg                    <  2:1371892/101‑1 (MQ=255)
                                                                                     cgATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATc                 <  1:2740243/101‑1 (MQ=255)
                                                                                      gATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCg                <  2:2462427/101‑1 (MQ=255)
                                                                                       aTCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGc               >  1:4252860/1‑101 (MQ=255)
                                                                                       aTCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGc               <  1:4480643/101‑1 (MQ=255)
                                                                                         cTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCCCGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGccc             <  1:867357/101‑1 (MQ=255)
                                                                                         cTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGccc             >  2:3326460/1‑101 (MQ=255)
                                                                                          tGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGcccc            >  2:4026181/1‑101 (MQ=255)
                                                                                          tGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGcccc            <  2:4047626/101‑1 (MQ=255)
                                                                                           gCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCg           >  2:3595419/1‑101 (MQ=255)
                                                                                              tCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTc        <  2:4680539/101‑1 (MQ=255)
                                                                                              tCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTc        <  2:4645057/101‑1 (MQ=255)
                                                                                              tCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTc        <  1:4324506/101‑1 (MQ=255)
                                                                                                atcttCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTCGa      <  2:2709155/97‑1 (MQ=255)
                                                                                                 cgcTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTCGAc     <  1:2667783/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GTCCCGATCCCCGCCATGGCCGAGACCGAAAGCTGGGCGATATCCGCGGTTTCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCGATCTGCGTCGCGCTCAGCCGTGCCGCCTGCGCCGTGGTCAGCGCCACGATGTCATTCGTGCCCAGCCCCTTCAGCGCCGCGGTGCTCAGCGCGCCGATCGCCCCGGTCGACAGC  >  NC_014034/678709‑678909

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