Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,947 C→G 13.4% R277G (CGC→GGC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9470C→G13.4% 108.3 ‑0.6 ‑1.7417R277G (CGC→GGC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (26/19):  ref base (156/146):  total (210/207)
Fisher's exact test strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.97e-02

Only 100 of 448 total aligned reads displayed.

CGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGC  >  NC_014034/1396849‑1397047
                                                                                                        |                                                                                                           
cgGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGc                                                                                                           >  1:4157232/1‑101 (MQ=255)
  ggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..t                                                                                                       <  1:4668657/101‑1 (MQ=255)
    ggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCg                                                                                                     >  2:4972547/1‑101 (MQ=255)
      cGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGc                                                                                                   >  2:3322895/1‑101 (MQ=255)
      cGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGc                                                                                                   >  1:2307705/1‑101 (MQ=255)
        tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCt.t                                                                                                >  2:958840/1‑101 (MQ=255)
          gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.Tcc                                                                                              <  1:1788990/101‑1 (MQ=255)
             cGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGc.c                                                                                          >  1:4975406/1‑101 (MQ=255)
                aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..c                                                                                     <  1:2135722/101‑1 (MQ=255)
                    ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CCG..CGTg.g                                                                                >  2:1399264/1‑101 (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                                                                             >  2:2677452/1‑101 (MQ=255)
                         agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCt                                                                           <  2:3512569/101‑1 (MQ=255)
                           agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcc                                                                         <  1:2909133/101‑1 (MQ=255)
                           agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcc                                                                         >  2:2929718/1‑101 (MQ=255)
                            gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GGCGCTccc                                                                        >  2:1724329/1‑101 (MQ=255)
                              tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCCC.CTG..CGTG.GTCGCTccccc                                                                      >  1:3955416/1‑101 (MQ=255)
                               aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                     >  2:3407247/1‑101 (MQ=255)
                                gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GGCGCTCCCCCcg                                                                    >  2:4440185/1‑101 (MQ=255)
                                     cgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCgggg                                                               <  1:1988471/101‑1 (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                                                            >  2:2226839/1‑101 (MQ=255)
                                         gccctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGTTCCCCCCGCGGGGGCGt                                                           <  1:1854325/99‑1 (MQ=255)
                                         gccctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                           <  2:3633955/99‑1 (MQ=255)
                                            tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                         <  2:409452/92‑1 (MQ=0)
                                            tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                         <  2:1455824/92‑1 (MQ=0)
                                            ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                         >  1:3474085/11‑101 (MQ=9)
                                            ctcctcAGACGTCTCGCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                        <  2:2635324/101‑1 (MQ=255)
                                            ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                        <  1:2808515/101‑1 (MQ=255)
                                            attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                        >  2:3242318/11‑101 (MQ=0)
                                             ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                       >  1:4088813/10‑101 (MQ=2)
                                              tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                      >  1:1501469/9‑101 (MQ=2)
                                              tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                      <  1:3217404/93‑1 (MQ=2)
                                                tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                    <  2:1290338/101‑1 (MQ=255)
                                                ggttgcGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                     <  2:1666920/96‑1 (MQ=2)
                                                 ccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   <  2:4591282/96‑1 (MQ=2)
                                                 cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   <  2:3157727/101‑1 (MQ=255)
                                                 cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                   >  1:3260658/1‑101 (MQ=255)
                                                 acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                    >  1:1521890/6‑101 (MQ=9)
                                                 acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                    >  1:4947849/6‑101 (MQ=9)
                                                  gcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                     >  2:4175570/4‑101 (MQ=2)
                                                  cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  <  1:903933/97‑1 (MQ=2)
                                                  aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGGGGGGGTCGaaa                                                  >  1:20549/1‑101 (MQ=255)
                                                   gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                 <  1:622073/101‑1 (MQ=255)
                                                   cccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                    <  1:4788053/99‑1 (MQ=2)
                                                    gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                 <  1:716399/100‑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)
                                                     cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  >  2:477790/1‑101 (MQ=2)
                                                     cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                  <  1:2838087/101‑1 (MQ=2)
                                                      gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                               <  2:1090364/101‑3 (MQ=2)
                                                      gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATc                                              >  2:1004889/1‑101 (MQ=255)
                                                      gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAActc                                               >  1:1057255/1‑99 (MQ=9)
                                                      gTCCCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                              <  1:3648862/101‑4 (MQ=2)
                                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                <  2:3735872/101‑1 (MQ=2)
                                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                >  1:4031554/1‑101 (MQ=2)
                                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATCa                                             >  1:3577883/1‑101 (MQ=255)
                                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                             >  1:4860951/1‑101 (MQ=255)
                                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                             <  1:2667094/101‑1 (MQ=255)
                                                       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.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                               >  1:1221739/1‑100 (MQ=2)
                                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                           <  1:44290/101‑7 (MQ=2)
                                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                           >  1:3393917/1‑95 (MQ=2)
                                                         gCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                            >  1:1927665/1‑96 (MQ=9)
                                                         ccTCCCCCCGCGGGGGCGGGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCacc                                           <  2:4909765/100‑3 (MQ=255)
                                                          ccccccccGCGGGGGCGTGGGTCCAAACATTGC.TT.TGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                           <  2:2766090/99‑7 (MQ=2)
                                                          cTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                             <  2:2053069/101‑1 (MQ=2)
                                                           tCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                            <  1:2944103/101‑2 (MQ=2)
                                                           cccccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                          <  1:4517929/100‑8 (MQ=9)
                                                            ccccccGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCTCCCGCGGGGGCGTGGATCGAAACACCaat                                           >  2:3077635/1‑99 (MQ=2)
                                                            ccccccGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCCCGGGGGCGTGGGTCGAAACATCAGCtcg                                        >  2:1958367/1‑101 (MQ=255)
                                                            ccccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                         >  2:3037009/1‑93 (MQ=9)
                                                             cccccGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                        <  1:1301550/101‑10 (MQ=9)
                                                              ccccGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CC.ATGAGTTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgc                                         <  2:4299223/101‑5 (MQ=0)
                                                              ccccGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                      <  2:1537314/101‑1 (MQ=255)
                                                                 cgcgGGGGCGTGGAGCGAACCA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                   <  2:799502/101‑1 (MQ=255)
                                                                  gggggggCGTGGATCGAAACA..GC.CT.GCGCG.GC.CCCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                  <  1:2195171/99‑1 (MQ=255)
                                                                  gcgGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                  >  1:3748730/1‑101 (MQ=255)
                                                                   cgGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                 >  1:639669/1‑101 (MQ=255)
                                                                   cgGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                 <  1:1394581/101‑1 (MQ=255)
                                                                    gggggCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                                >  1:4378169/1‑101 (MQ=255)
                                                                     ggggCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                               <  2:3998185/101‑1 (MQ=255)
                                                                     ggggCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATCAGCTCGTCGTTGGTc                               >  1:3835259/1‑101 (MQ=255)
                                                                          gTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCCTCGGTGGGCAGCCg                          >  2:4335797/1‑101 (MQ=255)
                                                                            ggATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                        >  1:3777720/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                   <  1:771636/101‑1 (MQ=255)
                                                                                 gAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg                   >  1:309201/1‑101 (MQ=255)
                                                                                   aaCA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                 <  1:1694266/101‑1 (MQ=255)
                                                                                    aCA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc                <  1:1439366/101‑1 (MQ=255)
                                                                                     cA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg               <  1:2595785/101‑1 (MQ=255)
                                                                                          c.cT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTcc            >  2:2398575/1‑100 (MQ=255)
                                                                                          c.cT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCg            >  1:2236265/1‑101 (MQ=255)
                                                                                            cT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc           >  2:3896688/1‑101 (MQ=255)
                                                                                              cgcgcG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGGCGTTGGTCAGCCGGACGATGTCCGTCGCt          <  1:2355598/100‑1 (MQ=255)
                                                                                                   g.gC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccccc     >  2:1131640/1‑101 (MQ=255)
                                                                                                   g.gC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccccc     >  2:121791/1‑101 (MQ=255)
                                                                                                     gC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc    <  1:2427407/101‑1 (MQ=255)
                                                                                                        |                                                                                                           
CGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGC  >  NC_014034/1396849‑1397047

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG