Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,968 G→C 8.4% G284R (GGT→CGT)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9680G→C 8.4% 108.2 ‑1.1 0.0694G284R (GGT→CGT) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (30/29):  ref base (280/354):  total (311/383)
Fisher's exact test strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 787 total aligned reads displayed.

CGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGC  >  NC_014034/1396868‑1397068
                                                                                                               |                                                                                                       
 ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTg.g                                                                                                       >  2:6841430/1‑101 (MQ=255)
 ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGGg.g                                                                                                       >  2:8228939/1‑101 (MQ=255)
        agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcc                                                                                                <  2:5389410/101‑1 (MQ=255)
        agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcc                                                                                                <  1:8229289/101‑1 (MQ=255)
         gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CCC.CTC..CGTG.GTCGCTccc                                                                                               >  1:2231997/1‑101 (MQ=255)
           tAGGCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTccccc                                                                                             <  1:4850070/101‑1 (MQ=255)
            aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcccccc                                                                                            >  1:2080588/1‑101 (MQ=255)
             gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCcg                                                                                           <  1:659122/101‑1 (MQ=255)
             gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCcg                                                                                           <  2:1459281/101‑1 (MQ=255)
                       aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                                 <  1:3977567/101‑1 (MQ=255)
                        cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                <  1:4468753/101‑1 (MQ=255)
                         ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                >  2:6640914/11‑101 (MQ=9)
                         ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.ATG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               >  1:4990518/1‑101 (MQ=255)
                         attggccccggTCGATCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               >  1:7177468/11‑101 (MQ=0)
                          ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                              <  2:7752115/92‑1 (MQ=2)
                          tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                              >  2:7050003/1‑101 (MQ=255)
                          tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                              <  1:7522927/101‑1 (MQ=255)
                          gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               >  1:7048394/10‑101 (MQ=9)
                          cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGG.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                               >  2:4605635/9‑101 (MQ=0)
                           tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  2:3962932/93‑1 (MQ=2)
                           tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  1:9159490/93‑1 (MQ=2)
                           gcccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                                                <  2:1013263/94‑1 (MQ=0)
                           cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGGG.GTCGCTCCCCCCGCGGGGGGGGGGATc                                                                             >  2:2934296/1‑101 (MQ=255)
                            ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                            <  2:3599676/94‑1 (MQ=2)
                            ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                            <  2:3416006/94‑1 (MQ=2)
                            cggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  2:495877/95‑1 (MQ=2)
                            cggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                             <  1:1694852/95‑1 (MQ=2)
                             tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                           >  2:7575275/1‑101 (MQ=255)
                             ccAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                           <  1:2598986/100‑1 (MQ=255)
                              gttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                           <  2:7267444/97‑1 (MQ=2)
                              gttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                           <  1:6315184/97‑1 (MQ=2)
                              ccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                          >  1:8125942/6‑101 (MQ=2)
                              ccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GCCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                          <  2:9026809/96‑1 (MQ=2)
                                gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  2:1626215/101‑1 (MQ=255)
                                 gcGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACATTGCTA.TGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  2:1359520/100‑1 (MQ=2)
                                 cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                       <  2:5659945/99‑2 (MQ=2)
                                 cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                                       >  2:4495934/3‑100 (MQ=2)
                                 ccgccGCCCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.T..CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                                      <  1:6856401/97‑1 (MQ=25)
                                  ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                                                      >  2:8619251/2‑99 (MQ=2)
                                  ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                                                      <  2:1139272/100‑3 (MQ=2)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcc                                                                       >  1:3376552/1‑100 (MQ=2)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcct                                                                      <  1:8780067/101‑3 (MQ=255)
                                  cGTCGCCGCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GCCC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                                                      <  1:1811975/101‑3 (MQ=2)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                        <  1:8415003/101‑1 (MQ=2)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                                                      >  1:7982959/1‑99 (MQ=2)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..CCCG.GCGCC.GC.AGCG.AGGGCG.TCTGGA.ACA..GGAT.GTCGCTCCCCCCGCGGNGGCGGGGATCGAAActa                                                                      >  2:2521233/1‑99 (MQ=2)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                                     >  1:8975475/1‑98 (MQ=2)
                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                                    >  2:7185895/1‑97 (MQ=9)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                                  >  2:2163257/1‑101 (MQ=255)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                                  >  1:3249739/1‑101 (MQ=255)
                                      ccTCCCCCCGCGGGGGCGTGGATCGAAACT..GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                                     <  2:9108524/100‑3 (MQ=2)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACT..GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                                    >  1:3131852/1‑101 (MQ=2)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACAT.GGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                                                  <  2:6348510/101‑7 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAACT..GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                                  <  1:2712219/101‑3 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                               >  2:8384948/1‑92 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAACA...ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                                               >  2:3111111/1‑92 (MQ=2)
                                         ccccccGCGGGGGCGTGGATCGAAAC...GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                                                <  2:2371519/101‑9 (MQ=9)
                                                 gggggCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGt                                                      >  1:8516486/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGt                                                      <  1:158091/101‑1 (MQ=255)
                                                   gggCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCa                                                    <  1:1013270/101‑1 (MQ=255)
                                                     gCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGc                                                  <  1:67661/101‑1 (MQ=255)
                                                      cGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGGT.TC.CGC.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGcc                                                 >  1:6464517/1‑101 (MQ=255)
                                                        gggATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCgg                                               <  2:7419924/100‑1 (MQ=255)
                                                          gATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAc                                             >  1:532960/1‑101 (MQ=255)
                                                           aTCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.g                                           >  2:581259/1‑101 (MQ=255)
                                                            tCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.Ga                                          >  2:2015337/1‑101 (MQ=255)
                                                             cGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GAt                                         <  1:7333994/101‑1 (MQ=255)
                                                              gAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATg                                        >  1:4717149/1‑101 (MQ=255)
                                                              gAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATg                                        <  1:7766983/101‑1 (MQ=255)
                                                               aaaCA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGt                                       >  2:8020223/1‑101 (MQ=255)
                                                                aaCA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTc                                      >  1:4187229/1‑101 (MQ=255)
                                                                aaCA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTc                                      >  1:8688806/1‑101 (MQ=255)
                                                                  cA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCg                                    >  1:2417826/1‑101 (MQ=255)
                                                                  cA..GCCG.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCg                                    <  2:607930/101‑1 (MQ=255)
                                                                       ccT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCg                                 >  2:878221/1‑101 (MQ=255)
                                                                         t.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCt                               <  1:5514456/101‑1 (MQ=255)
                                                                            cgcg.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTcc                             >  2:332854/1‑101 (MQ=255)
                                                                                  c.cGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCcg                        >  1:6176697/1‑101 (MQ=255)
                                                                                    cGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCcgc                       <  1:6016707/101‑1 (MQ=255)
                                                                                     gCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCcgcg                      >  1:5787840/1‑101 (MQ=255)
                                                                                       a.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCggg                    >  1:385417/1‑101 (MQ=255)
                                                                                             cT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGt               <  2:4271472/101‑1 (MQ=255)
                                                                                             cT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGt               <  2:4839036/101‑1 (MQ=255)
                                                                                             cT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGt               <  1:5668058/101‑1 (MQ=255)
                                                                                              t.tC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTg              <  2:5110396/101‑1 (MQ=255)
                                                                                                tC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTgg             <  1:5202468/101‑1 (MQ=255)
                                                                                                tC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTgg             <  1:1911323/101‑1 (MQ=255)
                                                                                                tC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTgg             >  2:410356/1‑101 (MQ=255)
                                                                                                 c.cGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGa            >  2:2709715/1‑101 (MQ=255)
                                                                                                   cGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGAt           <  1:1189256/101‑1 (MQ=255)
                                                                                                    gC.CTG..GGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATc          <  1:2934394/101‑1 (MQ=255)
                                                                                                     ggatcccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTg              >  2:3834152/11‑101 (MQ=2)
                                                                                                        tcccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt           <  2:7839174/94‑1 (MQ=2)
                                                                                                        tcccgccg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGAt           <  2:180247/94‑1 (MQ=2)
                                                                                                            cGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     <  2:6736504/101‑1 (MQ=255)
                                                                                                             gTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGGAACt    >  1:4290046/1‑101 (MQ=255)
                                                                                                             gTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACt    <  2:6519875/101‑1 (MQ=255)
                                                                                                              tG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg   <  1:356718/101‑1 (MQ=255)
                                                                                                              cg.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     >  2:5931798/2‑101 (MQ=2)
                                                                                                               g.gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTTGACCGCATCGAAGCTGGTCACATCGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcg    >  1:1346309/1‑100 (MQ=2)
                                                                                                               |                                                                                                       
CGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCG.TCGTTGGTCAGCCGGAC.GATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGC  >  NC_014034/1396868‑1397068

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