Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,951 T→A 14.8% I278N (ATC→AAC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9510T→A14.8% 285.2 ‑0.8 ‑0.3702I278N (ATC→AAC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (42/55):  ref base (256/282):  total (342/360)
Fisher's exact test strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.49e-01

Only 100 of 729 total aligned reads displayed.

GGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGG  >  NC_014034/1396851‑1397051
                                                                                                         |                                                                                                          
     gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.Ct                                                                                                      <  1:3783600/101‑1 (MQ=255)
      tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.Ct.t                                                                                                    >  1:4377929/1‑101 (MQ=255)
      tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.Ct.t                                                                                                    >  1:5604044/1‑101 (MQ=255)
      tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.Ct.t                                                                                                    >  1:1616203/1‑101 (MQ=255)
      tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.Ct.t                                                                                                    <  1:4210677/101‑1 (MQ=255)
       ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GT.CGCA.TCGG.CT.Tc                                                                                                   >  1:9078144/1‑101 (MQ=255)
       ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.Tc                                                                                                   >  1:6385334/1‑101 (MQ=255)
           cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGc.c                                                                                              >  2:67661/1‑101 (MQ=255)
           cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGc.c                                                                                              <  2:7516233/101‑1 (MQ=255)
             aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTg                                                                                            >  2:6863829/1‑101 (MQ=255)
                  ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CCTg.g                                                                                    >  2:1911263/1‑101 (MQ=255)
                   gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.Gt                                                                                   <  1:6221110/101‑1 (MQ=255)
                    tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTc                                                                                  <  1:8470796/101‑1 (MQ=255)
                     ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCg                                                                                 >  1:5110565/1‑101 (MQ=255)
                     ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGC.GGGGCCTGGATCGAAACA.GCCC.GCGCG.GG.CGCC.TCGG.CC.TCCCG.CTG..CGGg.ggccc                                                                                >  2:1361235/1‑97 (MQ=25)
                            tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                          <  2:757795/101‑1 (MQ=255)
                            tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTccccc                                                                          >  1:6748457/1‑101 (MQ=255)
                             aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                         >  2:7756054/1‑101 (MQ=255)
                              gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCcg                                                                        <  1:659122/101‑1 (MQ=255)
                                 cgcgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCgg                                                                     <  1:8148808/101‑1 (MQ=255)
                                    gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCCCggggg                                                                  >  2:5260768/1‑101 (MQ=255)
                                       gACCTCCTCAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                               <  2:4587213/101‑1 (MQ=255)
                                         cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                             <  1:9201562/101‑1 (MQ=255)
                                          tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                             >  2:8200495/10‑101 (MQ=0)
                                          attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  1:7074626/11‑101 (MQ=0)
                                           ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:3930523/92‑1 (MQ=2)
                                            tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  2:3962932/93‑1 (MQ=2)
                                            tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  1:7547085/9‑101 (MQ=2)
                                            tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:3232805/93‑1 (MQ=9)
                                            cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  2:8595222/1‑101 (MQ=255)
                                            cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  2:6965279/101‑1 (MQ=255)
                                            cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  1:7910785/1‑101 (MQ=255)
                                             cggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  2:1037128/95‑1 (MQ=2)
                                              ccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGac                                                           >  2:2442617/6‑100 (MQ=2)
                                              ccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  2:6397868/96‑1 (MQ=2)
                                              ccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:3376433/96‑1 (MQ=2)
                                               gttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        >  1:77640/5‑101 (MQ=2)
                                               gttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GGCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        >  2:4124103/5‑101 (MQ=2)
                                               cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       >  1:5777606/1‑101 (MQ=255)
                                               cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       >  1:1736816/1‑101 (MQ=255)
                                                gcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  1:6088371/98‑1 (MQ=2)
                                                cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      <  2:5983736/97‑1 (MQ=2)
                                                  cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                    <  1:8271636/99‑2 (MQ=2)
                                                  cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                    <  2:5659945/99‑2 (MQ=2)
                                                  ccgccGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.T.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                   <  1:6856401/97‑1 (MQ=25)
                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                    >  1:8050018/2‑100 (MQ=9)
                                                   cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcc                                                    >  1:2312214/1‑100 (MQ=2)
                                                   cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                   <  1:5923102/101‑1 (MQ=255)
                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAAccc                                                   >  2:8765360/1‑99 (MQ=2)
                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                                  <  1:4579230/101‑1 (MQ=255)
                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                                  <  1:2093368/101‑4 (MQ=2)
                                                    gTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                   <  1:7344510/101‑3 (MQ=9)
                                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                 <  2:2304327/101‑1 (MQ=255)
                                                     tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg                                                  >  1:8763255/1‑98 (MQ=9)
                                                     tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                 <  2:5389740/101‑1 (MQ=255)
                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                   <  2:891981/101‑2 (MQ=2)
                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                 <  2:4967170/101‑5 (MQ=9)
                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                 >  2:7185895/1‑97 (MQ=9)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                               >  2:1462616/1‑95 (MQ=2)
                                                       ccTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                  <  2:9108524/100‑3 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                 <  1:5311414/101‑1 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                 <  2:7265493/101‑1 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAG.GT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                               <  1:1834552/101‑7 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCt                                              <  2:5789628/101‑1 (MQ=255)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                               >  1:2771859/1‑95 (MQ=9)
                                                         tCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                             >  2:6569527/1‑101 (MQ=255)
                                                         tCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACATG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                             <  2:4368828/101‑9 (MQ=0)
                                                         cccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                             <  2:359820/100‑1 (MQ=255)
                                                          ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                               >  1:5499544/1‑99 (MQ=2)
                                                          ccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGG.CC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                            >  2:3111111/1‑92 (MQ=2)
                                                          ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                             >  1:3481028/1‑93 (MQ=9)
                                                          ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGG.ATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                             <  2:8300755/101‑9 (MQ=9)
                                                            ccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCAT.CG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgc                                             >  1:857663/1‑97 (MQ=2)
                                                             cccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtcg                                         >  1:4839196/1‑101 (MQ=255)
                                                             cccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCCTATGCGC.CGGA.CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACAACaccggggcc                                         >  2:2979755/1‑93 (MQ=2)
                                                                 cgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGCTgg                                     >  2:2633654/1‑101 (MQ=255)
                                                                  gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                    <  1:1942063/101‑1 (MQ=255)
                                                                    gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                  <  1:1013270/101‑1 (MQ=255)
                                                                      gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                                <  2:1652766/101‑1 (MQ=255)
                                                                       cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                               >  2:3400048/1‑101 (MQ=255)
                                                                        gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                              <  2:5855154/101‑1 (MQ=255)
                                                                        gTGCATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATAAGCTCGTCGTTGGTCAGCCg                              >  1:7202431/1‑101 (MQ=255)
                                                                         tGGATCGAAACA.GCCT.GCGCG.GC.C.CA.TCGG.CT.TCCGC.CCG..CCGG.GGCGCCCCCCCCGCGGGGGCGGGGGTCGGAACATCAGCTCGGCGTTGGGCAGCCCGa                            >  1:1538122/1‑101 (MQ=255)
                                                                          ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCCGCCGGa                            >  2:6146303/1‑101 (MQ=255)
                                                                          ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGGCAGCCGGa                            >  2:6696335/1‑101 (MQ=255)
                                                                            aTCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGACCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                          <  1:6473366/101‑1 (MQ=255)
                                                                              cGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                        >  1:2434568/1‑101 (MQ=255)
                                                                               gAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                       >  1:4092171/1‑101 (MQ=255)
                                                                               gAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCCTTGGTCAGGCGGACGCTg                       >  2:3744144/1‑101 (MQ=255)
                                                                                aaaCA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                      >  2:1190333/1‑101 (MQ=255)
                                                                                  aCA.GCCT.GCACG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                    >  2:2102330/1‑101 (MQ=255)
                                                                                    a.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                  <  2:7876987/101‑1 (MQ=255)
                                                                                    a.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGGTGGTCACACGGACGATGTCCGt                  >  2:7356877/1‑101 (MQ=255)
                                                                                            cgcg.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcc            >  2:332854/1‑101 (MQ=255)
                                                                                              cg.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccc          >  1:1227867/1‑101 (MQ=255)
                                                                                                  c.cGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcg       >  1:4986870/1‑101 (MQ=255)
                                                                                                  c.cGCA.GCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcg       <  1:2633736/101‑1 (MQ=255)
                                                                                                       a.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCCTCGCGCCCCCCGCggg   >  1:8949293/1‑101 (MQ=255)
                                                                                                         tCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg  >  2:4843964/1‑101 (MQ=255)
                                                                                                         tCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg  <  1:3057276/101‑1 (MQ=255)
                                                                                                         |                                                                                                          
GGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGG.CT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGG  >  NC_014034/1396851‑1397051

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