Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,953 G→A 16.5% G279S (GGC→ATA)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9530G→A16.5% 190.2 ‑2.8 ‑1.7415G279S (GGC→AGC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (27/41):  ref base (181/163):  total (209/206)
Fisher's exact test strand distribution p-value = 6.29e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.92e-02

Only 100 of 447 total aligned reads displayed.

GGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1396853‑1397053
                                                                                                           |                                                                                                         
    tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCt.t                                                                                                      >  2:958840/1‑101 (MQ=255)
    tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCt.t                                                                                                      <  1:1484668/101‑1 (MQ=255)
     tgATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.Tc                                                                                                     <  2:3670818/100‑1 (MQ=255)
     ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.Tc                                                                                                     >  2:1468206/1‑101 (MQ=255)
          gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGGCGCTCCTCCCGCGGGGACGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                               >  2:1654952/1‑101 (MQ=255)
          gAAACCGGCGCAGAGATAGACCGCGCGTGCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                               <  2:3774023/101‑1 (MQ=38)
           aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                              >  2:1988327/1‑101 (MQ=255)
            aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..c                                                                                           <  1:2043143/101‑1 (MQ=255)
                   ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TGCGC.CCG..CGTG.GTCg                                                                                   >  2:2844613/1‑101 (MQ=255)
                   ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                   >  1:3157978/1‑101 (MQ=255)
                   ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                   >  2:2677452/1‑101 (MQ=255)
                   ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCg                                                                                   >  2:4513459/1‑101 (MQ=255)
                     agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                                 >  2:4417943/1‑101 (MQ=255)
                      gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                                >  2:1116184/1‑101 (MQ=255)
                       agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcc                                                                               >  2:2929718/1‑101 (MQ=255)
                             aCCGCGCGGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgc                                                                         <  2:3948293/101‑1 (MQ=255)
                                   cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTC..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                   >  1:2130045/1‑101 (MQ=255)
                                    ggACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCg                                                                  >  2:1847137/1‑101 (MQ=255)
                                    ggACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCg                                                                  >  1:2885563/1‑101 (MQ=255)
                                      aCCTCCCCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                                <  2:3189345/101‑1 (MQ=255)
                                       cctcctCAGACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CCCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                               <  2:4993851/101‑1 (MQ=255)
                                        tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                               <  2:4033074/92‑1 (MQ=0)
                                        ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                               >  1:4002526/11‑101 (MQ=9)
                                        ctcctcAGACGTCTCGCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                              <  2:2635324/101‑1 (MQ=255)
                                        attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                              >  2:3242318/11‑101 (MQ=0)
                                         tcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                             <  1:2447335/101‑1 (MQ=255)
                                         cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                              >  1:308969/9‑101 (MQ=2)
                                          tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGGGTgggc                                                             >  2:4596172/9‑99 (MQ=9)
                                           ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGTT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAGCg                                                           >  1:3701711/1‑101 (MQ=255)
                                           ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           <  1:2816395/101‑1 (MQ=255)
                                           cggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                            <  2:4132904/95‑1 (MQ=2)
                                            ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           >  1:2327536/6‑101 (MQ=2)
                                            gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                          >  2:3462617/7‑101 (MQ=2)
                                            aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                           >  1:4290267/7‑101 (MQ=9)
                                            aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCCCCCCCCGCGGGGGGGGGGGTCg                                                           >  1:3902051/7‑101 (MQ=2)
                                             cgcccggcGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                            <  2:4786259/94‑1 (MQ=2)
                                             cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                         <  1:1845182/101‑1 (MQ=255)
                                             cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                                         >  1:635910/1‑101 (MQ=255)
                                              cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                        <  1:903933/97‑1 (MQ=2)
                                               gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTTGCTCCCCCCGCGGGGGCGGGGATCGAAAc                                                       >  1:747289/1‑101 (MQ=255)
                                               gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                       <  1:3771865/101‑1 (MQ=255)
                                               cccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                          <  1:4788053/99‑1 (MQ=2)
                                                ctgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                       <  2:230965/99‑1 (MQ=9)
                                                cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                      <  2:4723137/99‑2 (MQ=2)
                                                aCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                      <  2:3589981/101‑1 (MQ=255)
                                                 tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                      >  1:1780662/2‑100 (MQ=9)
                                                 ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                                     >  1:4124903/2‑99 (MQ=2)
                                                 ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                                     <  1:2707310/100‑3 (MQ=2)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                       >  1:440269/1‑101 (MQ=2)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                                     <  2:1090364/101‑3 (MQ=2)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAActc                                                     >  1:1057255/1‑99 (MQ=9)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                     <  1:989258/101‑3 (MQ=9)
                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGCGGATCGAAActc                                                     >  1:1664245/1‑99 (MQ=9)
                                                   tcccTCCCCCCGCGGGGGCGTGGATCGAGACC..GC.CC.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                   <  1:2922142/98‑1 (MQ=25)
                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                      <  2:1261210/101‑1 (MQ=2)
                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCACCGCGGGGGCGTGGATCGAAActaca                                                   >  1:4238981/1‑97 (MQ=2)
                                                   tCGCTCCACCCGCGGGGTCGTGGGCCGAAACA..ACGCG.GCGCG.GG.CCTA.T.GGCTGTCGCC.CCC..CCCC.GTCGCTCCCCCCCCGGGGGCGTGGATCGAAActgc                                                    <  2:710065/101‑4 (MQ=2)
                                                   gcGCGCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                   <  2:624089/100‑1 (MQ=255)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                     >  1:1221739/1‑100 (MQ=2)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                                     <  2:1704308/101‑2 (MQ=2)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                   >  1:2176174/1‑97 (MQ=9)
                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCC.CCGGATCTTTGT.CTG..CGCG.GGCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                   >  1:2235468/1‑97 (MQ=9)
                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                    <  1:1913275/101‑3 (MQ=2)
                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                 >  1:3242347/1‑101 (MQ=255)
                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                 <  2:1447142/101‑7 (MQ=2)
                                                      ccccccccGCGGGGGCGTGGGTCCAAACATTGC.TT.TGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                                 <  2:2766090/99‑7 (MQ=2)
                                                      cacccccgGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCt                                                <  1:4995760/99‑1 (MQ=255)
                                                       tCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                               >  1:3457895/1‑93 (MQ=2)
                                                       tCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGGGTGGATCGAAActcggcct                                                >  1:4080572/1‑94 (MQ=9)
                                                       cccccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                <  1:4517929/100‑8 (MQ=9)
                                                        ccccccGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCTCCCGCGGGGGCGTGGATCGAAACACCaat                                                 >  2:3077635/1‑99 (MQ=2)
                                                        ccccccGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                 <  1:1519451/101‑3 (MQ=2)
                                                        ccccccGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcg                                              <  1:2642444/101‑1 (MQ=255)
                                                        ccccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                               >  2:3037009/1‑93 (MQ=9)
                                                         cccccGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgc                                              >  1:2290327/1‑92 (MQ=2)
                                                         cccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                              <  2:4921676/101‑10 (MQ=9)
                                                         cccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                              <  1:1301550/101‑10 (MQ=9)
                                                          ccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                             <  1:3534010/101‑11 (MQ=2)
                                                            ccGGGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTcgtcgt                                          <  1:682592/101‑1 (MQ=255)
                                                             cgcgGGGGCGTGGAGCGAACCA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                         <  2:799502/101‑1 (MQ=255)
                                                              gcgGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGATg                                        >  2:4255975/1‑101 (MQ=255)
                                                               cgGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                       <  1:1394581/101‑1 (MQ=255)
                                                               cgGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGGG.GTCGCTCCCCCCGCGGGGGCGTGGATCGGAACATCAGCTCGTCGTTgg                                       >  1:3792349/1‑101 (MQ=255)
                                                                 ggggGGTGGATCCAAACA..GC.CT.GCGCG.GC.CGCA.TCGGTT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                                     <  1:4875229/101‑1 (MQ=255)
                                                                      gTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                <  1:3623409/101‑1 (MQ=255)
                                                                      gTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                >  1:1970120/1‑101 (MQ=255)
                                                                           tCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGa                           >  2:1353300/1‑101 (MQ=255)
                                                                             gAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                         >  2:4320729/1‑101 (MQ=255)
                                                                              aaaCA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                        >  1:4963719/1‑101 (MQ=255)
                                                                                aCA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                      >  2:1913806/1‑101 (MQ=255)
                                                                                aCA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                      <  1:1439366/101‑1 (MQ=255)
                                                                                  a..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                    >  1:2400336/1‑101 (MQ=255)
                                                                                      c.cT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTcc                  >  2:2398575/1‑100 (MQ=255)
                                                                                             gcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc             >  1:338754/1‑101 (MQ=255)
                                                                                             gcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc             >  2:3328568/1‑101 (MQ=255)
                                                                                               g.gC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccccc           >  2:121791/1‑101 (MQ=255)
                                                                                                     gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg       >  2:1396510/1‑101 (MQ=255)
                                                                                                     gCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg       <  1:4434090/101‑1 (MQ=255)
                                                                                                        ctCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggg     <  2:2185965/100‑1 (MQ=255)
                                                                                                          cGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg   >  1:4604600/1‑101 (MQ=255)
                                                                                                           |                                                                                                         
GGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGC  >  NC_014034/1396853‑1397053

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