Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,945 Δ1 bp 11.1% coding (827/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9450G→.11.1% 138.8 ‑0.6 ‑1.6412coding (827/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (26/19):  ref base (171/156):  total (210/202)
Fisher's exact test strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.38e-02

Only 100 of 447 total aligned reads displayed.

CCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCC  >  NC_014034/1396845‑1397045
                                                                                                         |                                                                                                            
 ccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGGAACA..GC.CT.GCGCG.Gc                                                                                                             >  2:300366/1‑101 (MQ=255)
      ggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGGCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..t                                                                                                     <  1:4976689/101‑1 (MQ=255)
        ggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCg                                                                                                   >  2:4972547/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)
             ggATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.Tc                                                                                             <  2:1460631/101‑1 (MQ=255)
              gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.Tcc                                                                                            <  1:646784/101‑1 (MQ=255)
              gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.Tcc                                                                                            <  1:1788990/101‑1 (MQ=255)
                    aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..c                                                                                   >  1:4145564/1‑101 (MQ=255)
                    aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..c                                                                                   <  1:2135722/101‑1 (MQ=255)
                     aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..Cg                                                                                  >  2:1096149/1‑101 (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)
                                gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GGCGCTccc                                                                      >  2:1724329/1‑101 (MQ=255)
                                 aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTC..CGTG.GGCGCTcccc                                                                     >  2:1226999/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.GTCGCTCCCCCcg                                                                  <  1:2897212/101‑1 (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: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)
                                               ctcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                                        <  1:4689136/91‑1 (MQ=0)
                                                tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                       <  2:409452/92‑1 (MQ=0)
                                                ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                       >  1:3474085/11‑101 (MQ=9)
                                                ggtcaacactgTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                       <  2:3587270/91‑1 (MQ=9)
                                                ctcctcAGACGTCTCGCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                      <  2:2635324/101‑1 (MQ=255)
                                                 ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                     >  1:4088813/10‑101 (MQ=2)
                                                 tcctcAGGCGTCCCGCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                     <  1:3091438/101‑1 (MQ=255)
                                                  tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                    <  1:3217404/93‑1 (MQ=2)
                                                  tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                    >  1:1501469/9‑101 (MQ=2)
                                                  cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                    <  1:936471/101‑1 (MQ=255)
                                                    tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                  <  2:1290338/101‑1 (MQ=255)
                                                     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                                                 >  1:3260658/1‑101 (MQ=255)
                                                     cAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                 <  2:3157727/101‑1 (MQ=255)
                                                     acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                  >  1:4947849/6‑101 (MQ=9)
                                                     acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                  >  1:1521890/6‑101 (MQ=9)
                                                      gcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                   >  2:4175570/4‑101 (MQ=2)
                                                      aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGGGGGGGTCGaaa                                                >  1:20549/1‑101 (MQ=255)
                                                      aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                >  2:185027/1‑101 (MQ=255)
                                                      aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGaaa                                                >  1:429884/1‑101 (MQ=255)
                                                       tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAT.GG.CTATGGCG.GC.AGTG..TCAGGT.T.TTC.CAG..CGGC.GGCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                >  2:290430/3‑101 (MQ=2)
                                                       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:3058242/3‑100 (MQ=2)
                                                        cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                              <  2:4723137/99‑2 (MQ=2)
                                                         cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                <  1:2838087/101‑1 (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..CGAGCGGCGCTCCCCCCGCGGGGCCGTGGATCGaaa                                                >  2:1677022/1‑101 (MQ=2)
                                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATc                                            >  2:1004889/1‑101 (MQ=255)
                                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActac                                            >  1:3824646/1‑98 (MQ=2)
                                                          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)
                                                           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: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:3393917/1‑95 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACA...A.CT.GACC...C.GGCA..ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                         <  1:44290/101‑7 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAAC...GC.TG.GCGCGCGC.CGCG..CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                          >  1:1927665/1‑96 (MQ=9)
                                                             ccTCCCCCCGCGGGGGCGTGGATCGAAACT..GC.GT.GCGCG.CCATGAAT.TCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                            <  2:3235949/100‑3 (MQ=2)
                                                             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)
                                                               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)
                                                                ccccccGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCCCGGGGGCGTGGGTCGAAACATCAGCtcg                                      >  2:1958367/1‑101 (MQ=255)
                                                                 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:1394581/101‑1 (MQ=255)
                                                                         ggggGGTGGATCCAAACA..GC.CT.GCGCG.GC.CGCA..TCGGTT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                             <  1:4875229/101‑1 (MQ=255)
                                                                         ggggCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                             <  2:3998185/101‑1 (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:4006472/1‑101 (MQ=255)
                                                                                ggATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                      >  1:3777720/1‑101 (MQ=255)
                                                                                  aTCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                    <  1:2350078/101‑1 (MQ=255)
                                                                                    cGAAACA..GC.CG.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                  <  2:2069792/101‑1 (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:2536770/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)
                                                                                             gC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc           <  1:2753155/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)
                                                                                                   gggcg.GC.CGCA..TCGGCT.TCCGC.CTG..CGGG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc       <  2:3748452/99‑1 (MQ=255)
                                                                                                   gcgcG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc       >  2:1358561/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:1131640/1‑101 (MQ=255)
                                                                                                         gC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc  <  1:2427407/101‑1 (MQ=255)
                                                                                                         |                                                                                                            
CCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..GC.CT.GCGCG.GC.CGCA..TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCC  >  NC_014034/1396845‑1397045

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