Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 681,562 G→T 2.4% A682E (GCG→GAG)  RCAP_rcc00630 ← ice nucleation protein repeat family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034681,5620G→T 2.4% 2.8 ‑0.5 ‑0.2167A682E (GCG→GAG) RCAP_rcc00630ice nucleation protein repeat family protein
Reads supporting (aligned to +/- strand):  new base (3/1):  ref base (88/75):  total (91/76)
Fisher's exact test strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.95e-01

Only 100 of 199 total aligned reads displayed.

CTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCT  >  NC_014034/681462‑681662
                                                                                                    |                                                                                                    
 tGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATcgccgc                                                                                                     <  1:1083017/101‑1 (MQ=255)
     tGCTGGTCAGCGCGCAGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTg                                                                                                 >  1:3566822/1‑101 (MQ=255)
      gCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGc                                                                                                >  2:2263147/1‑101 (MQ=255)
       cTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGcc                                                                                               >  2:2417968/1‑101 (MQ=255)
        tGGTCAGCGCGCCGATCTGCCCCGAGCCCAGACGCCCCGCCTGCGCCGTCGTCAGCGTCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGccc                                                                                              <  1:3620554/101‑1 (MQ=255)
           tCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc                                                                                           <  1:3281936/101‑1 (MQ=255)
           tCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgc                                                                                           <  2:3169977/101‑1 (MQ=255)
            cAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAgcc                                                                                          >  1:4102368/1‑100 (MQ=255)
               cgcgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccg                                                                                       >  2:3988488/1‑101 (MQ=255)
               cccgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGcgccg                                                                                       <  2:1861293/99‑1 (MQ=255)
                 cgcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCgccgcc                                                                                     <  1:1056279/101‑1 (MQ=255)
                  gcCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCa                                                                                    <  1:2318619/101‑1 (MQ=255)
                    cGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCAcc                                                                                  <  2:2240735/101‑1 (MQ=255)
                      aTCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGt                                                                                >  1:3597699/1‑101 (MQ=255)
                       tCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTc                                                                               <  2:294444/101‑1 (MQ=255)
                        cTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCg                                                                              >  2:895222/1‑101 (MQ=255)
                         tGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGCCAGCGCCACGATGTCATCGGGGCCAAAGCCGCGGATCGCCGCCGTGCCCAGTGCCGCCACCGTCGc                                                                             >  1:3951717/1‑101 (MQ=255)
                          gcgcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGcc                                                                            >  1:4129272/1‑101 (MQ=255)
                            gcCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGa                                                                          >  1:3336373/1‑101 (MQ=255)
                             ccGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGaa                                                                         <  2:4061610/101‑1 (MQ=255)
                               gAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGc                                                                       >  2:2749677/1‑101 (MQ=255)
                                  cccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAg                                                                    <  1:351609/101‑1 (MQ=255)
                                   ccAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCGAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgc                                                                   >  2:3789613/1‑101 (MQ=255)
                                    cgtgcgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgcg                                                                  <  2:3566569/97‑1 (MQ=255)
                                     aGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAgcgc                                                                 >  1:2600665/1‑101 (MQ=255)
                                      gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc                                                                >  1:2367742/1‑101 (MQ=255)
                                      gACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgcc                                                                <  1:1376051/101‑1 (MQ=255)
                                        cgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccgc                                                              <  2:3350183/101‑1 (MQ=255)
                                        cgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccgc                                                              <  2:130634/101‑1 (MQ=255)
                                        cgcgCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGcgccgc                                                              >  1:2880623/1‑101 (MQ=255)
                                          cgccgcCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGa                                                            >  1:1896006/1‑101 (MQ=255)
                                           gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGcgccgcaag                                                           >  1:3173664/1‑98 (MQ=25)
                                           gccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGac                                                           >  2:1135505/1‑100 (MQ=255)
                                           cccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGAt                                                           <  1:2661722/100‑1 (MQ=255)
                                            ccgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATc                                                          >  2:870953/1‑101 (MQ=255)
                                             cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                                         <  1:533685/101‑1 (MQ=255)
                                             cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                                         <  1:2728077/101‑1 (MQ=255)
                                             cgccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAACCCGCGGATCCCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCt                                                         <  1:2534983/101‑1 (MQ=255)
                                              gccTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTg                                                        <  1:2309538/101‑1 (MQ=255)
                                               ccTGCGCCGTCGTCAGCGCCACGATGTCATCGGAGCCGAGGCCCTTGATCGCATTGGCACTCAGCGCCGTCACGCTGGCCGTATTCAGCGCCGCAAGCTgc                                                       <  1:3673976/101‑1 (MQ=37)
                                               ccTGCGCCGGCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgc                                                       <  1:836827/101‑1 (MQ=255)
                                                cTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTgcg                                                      <  1:3372290/101‑1 (MQ=255)
                                                   cgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCg                                                   >  2:2886928/1‑101 (MQ=255)
                                                   cgcCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCg                                                   >  1:552372/1‑101 (MQ=255)
                                                     ccGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTg                                                 <  2:748842/101‑1 (MQ=255)
                                                      cgtcgtCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGc                                                >  1:2129992/1‑101 (MQ=255)
                                                       gtcgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCt                                               <  1:1220897/101‑1 (MQ=255)
                                                       gtcgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCt                                               >  1:792352/1‑101 (MQ=255)
                                                         cttcAGCGCCCCGTTGTCCCCGGCGCCTAAGCCTCGGGTCGCCACCGTGCCCAGCGCCGCCACCGTCTCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                             <  2:2076227/99‑1 (MQ=25)
                                                         cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                             <  2:3502653/101‑1 (MQ=255)
                                                         cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                             >  2:415676/1‑101 (MQ=255)
                                                         cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCACCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                             <  2:3660709/101‑1 (MQ=255)
                                                         cgtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCACCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCa                                             >  1:4174621/1‑101 (MQ=255)
                                                          gtcAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAg                                            <  2:2460018/101‑1 (MQ=255)
                                                            cAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCACCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAgcg                                          >  2:1680611/1‑101 (MQ=255)
                                                              gcgcCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgcc                                        >  2:2501441/1‑101 (MQ=255)
                                                                  cACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgccg                                    >  1:3081011/1‑101 (MQ=255)
                                                                  cACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGcgccgccg                                    >  2:3893902/1‑101 (MQ=255)
                                                                     gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                 <  2:2348344/101‑1 (MQ=255)
                                                                     gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                 >  2:3263484/1‑101 (MQ=255)
                                                                     gATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCt                                 <  1:2442011/101‑1 (MQ=255)
                                                                      aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                                >  2:2793438/1‑101 (MQ=255)
                                                                      aTGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTg                                >  1:2691492/1‑101 (MQ=255)
                                                                       tGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGc                               >  1:2885222/1‑101 (MQ=255)
                                                                       tGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGc                               >  2:267250/1‑101 (MQ=255)
                                                                        gTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCa                              >  2:500491/1‑101 (MQ=255)
                                                                        gTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCa                              >  1:2954979/1‑101 (MQ=255)
                                                                         tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                             >  1:707031/1‑101 (MQ=255)
                                                                         tCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCAc                             <  1:2234708/101‑1 (MQ=255)
                                                                           aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                           <  1:1787075/101‑1 (MQ=255)
                                                                           aTCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCg                           <  2:3136447/101‑1 (MQ=255)
                                                                             cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                         >  1:3874009/1‑101 (MQ=255)
                                                                             cGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAg                         >  1:156243/1‑101 (MQ=255)
                                                                              cttGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGc                        <  2:2171762/99‑1 (MQ=255)
                                                                               ggcccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGcc                       <  2:1505412/98‑1 (MQ=255)
                                                                               gTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGcc                       >  1:86128/1‑101 (MQ=255)
                                                                                tGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGcag                      >  1:2171907/1‑99 (MQ=255)
                                                                                tGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCa                      >  1:2476392/1‑101 (MQ=255)
                                                                                 gcaaacgCCGCGGACCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCaa                     <  2:969250/95‑1 (MQ=255)
                                                                                 cccaaccccGCGGATCCCCGCCGTTCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCaa                     <  1:622343/94‑1 (MQ=38)
                                                                                  ccAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAg                    >  2:4196249/1‑101 (MQ=255)
                                                                                        ccGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGc              >  2:807441/1‑101 (MQ=255)
                                                                                          gcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGGGCCGTGCTCAGCGCCGCCGCCTGCCCCGAGCCAAGCGCGGCCa            >  2:1065554/1‑101 (MQ=255)
                                                                                          gcgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCa            >  2:2200436/1‑101 (MQ=255)
                                                                                           cgcctCGCCGCCGTGCCCGGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCGGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAc           <  2:1953634/97‑1 (MQ=255)
                                                                                           cgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAc           <  1:509295/101‑1 (MQ=255)
                                                                                           cgGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCAc           >  2:1258087/1‑101 (MQ=255)
                                                                                             gaccgccgcCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCg         <  2:3336128/98‑1 (MQ=255)
                                                                                              aTCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGc        >  1:2750639/1‑101 (MQ=255)
                                                                                              aTCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGc        <  2:706954/101‑1 (MQ=255)
                                                                                               tggccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCa       <  2:637060/99‑1 (MQ=255)
                                                                                                cgccgcCGTGCCCAGCGCCGCCACCGTCGCTGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCaa      >  2:261754/1‑101 (MQ=255)
                                                                                                cgccgcCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCaa      >  1:2700687/1‑101 (MQ=255)
                                                                                                 gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg     <  2:220353/101‑1 (MQ=255)
                                                                                                 gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg     >  1:2228822/1‑101 (MQ=255)
                                                                                                 gccgccGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAg     <  1:256453/101‑1 (MQ=255)
                                                                                                   cgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGc   <  1:2900303/101‑1 (MQ=255)
                                                                                                   cgccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGc   >  2:2697249/1‑101 (MQ=255)
                                                                                                    gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCt  <  1:370174/101‑1 (MQ=255)
                                                                                                    gccgTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCACGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCt  <  2:3329643/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
CTGCGTGCTGGTCAGCGCGCCGATCTGCGCCGAGCCCAGACGCGCCGCCTGCGCCGTCGTCAGCGCCACGATGTCATCGGTGCCAAAGCCGCGGATCGCCGCCGTGCCCAGCGCCGCCACCGTCGCCGAAGCCAGCGCCGCGATCTGCGCCGTGCTCAGCGCCGCCGCCTGCACCGAGCCAAGCCCGGCCACCGCAAGGCT  >  NC_014034/681462‑681662

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