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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9360G→.10.9% 122.9 ‑2.0 ‑0.1368coding (818/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (25/15):  ref base (161/163):  total (188/180)
Fisher's exact test strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.12e-01

Only 100 of 386 total aligned reads displayed.

GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTC  >  NC_014034/1396836‑1397036
                                                                                                     |                                                                                                             
gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.g                                                                                                               <  2:3019613/101‑1 (MQ=255)
gTCGCCCCCCCCGCGGGGGCGTGGATCGAANCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.g                                                                                                               <  1:398097/101‑1 (MQ=255)
   ccTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCt                                                                                                            <  2:1796631/100‑1 (MQ=255)
    cTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.g                                                                                                          <  1:2376743/101‑1 (MQ=255)
        ccccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.Gcgcg                                                                                                      >  1:800807/1‑101 (MQ=255)
          ccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.Gc                                                                                                   >  2:1720113/1‑101 (MQ=255)
              gggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCa                                                                                              >  2:465649/1‑101 (MQ=255)
               ggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.t                                                                                            <  2:77436/101‑1 (MQ=255)
                 ggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCg                                                                                          >  2:1512362/1‑101 (MQ=255)
                   cGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGc                                                                                        >  1:14634/1‑101 (MQ=255)
                     tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt.t                                                                                     >  2:3586799/1‑101 (MQ=255)
                     tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt.t                                                                                     <  1:427397/101‑1 (MQ=255)
                      ggAGCGAAACCGGTCCAGAGATCGACCGCGCGGACCTCCTCAGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.Tc                                                                                    <  2:1749438/101‑1 (MQ=255)
                          cGAAACCGGTCCAGAGATAGCCCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGc.c                                                                               <  1:3064231/101‑1 (MQ=255)
                           gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                              >  1:2122624/1‑101 (MQ=255)
                            aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                             >  2:318769/1‑101 (MQ=255)
                             aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..c                                                                          >  1:2806305/1‑101 (MQ=255)
                              aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..Cg                                                                         <  1:1370727/101‑1 (MQ=255)
                                cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg                                                                       >  1:1571533/1‑101 (MQ=255)
                                cGGTCCAGAGAGAGACCGCGCGGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg                                                                       <  2:4421564/101‑1 (MQ=255)
                                 ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg.g                                                                     <  2:55069/101‑1 (MQ=255)
                                  gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGCGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGCg.gg                                                                    >  1:3764269/1‑100 (MQ=255)
                                   tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTc                                                                   <  2:953218/101‑1 (MQ=255)
                                     cAGAGATAGGCCGCGCGGCCCTCCTCAGCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGc                                                                 <  1:2306858/101‑1 (MQ=255)
                                     cAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGc                                                                 >  1:897592/1‑101 (MQ=255)
                                       gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                               <  2:3834154/101‑1 (MQ=255)
                                        agaTAGACCGCGCGGACCTCCTCAGACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcc                                                              <  1:2754403/101‑1 (MQ=255)
                                           tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                           <  2:2610443/101‑1 (MQ=255)
                                           tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccccc                                                           <  1:1212332/101‑1 (MQ=255)
                                            aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                          <  1:3524373/101‑1 (MQ=255)
                                             gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcg                                                         <  2:244993/101‑1 (MQ=255)
                                               ccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCcgcg                                                       >  2:1613305/1‑101 (MQ=255)
                                                   gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCggggg                                                   >  1:1047280/1‑101 (MQ=255)
                                                        ctcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                               <  2:2174672/91‑1 (MQ=0)
                                                         ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                              <  2:619402/91‑1 (MQ=9)
                                                         attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                             >  1:2879910/11‑101 (MQ=0)
                                                         attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                             >  1:4402835/11‑101 (MQ=0)
                                                          gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                             >  1:100646/10‑101 (MQ=9)
                                                          gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                             >  1:3687396/10‑101 (MQ=9)
                                                          gtcaacactgTCGCCCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                             <  2:1450330/92‑1 (MQ=9)
                                                          cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGa                                             >  1:657/9‑101 (MQ=0)
                                                           tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                           >  2:1302922/9‑101 (MQ=2)
                                                           cctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  2:2408563/101‑1 (MQ=255)
                                                             ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                          >  2:4317839/6‑101 (MQ=2)
                                                             gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                         >  1:4119223/7‑101 (MQ=2)
                                                             aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCATGGATCg                                          <  1:3871980/95‑1 (MQ=9)
                                                              acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                         >  1:331521/6‑101 (MQ=9)
                                                               cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       >  2:3522513/5‑101 (MQ=2)
                                                               aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       <  1:318804/101‑1 (MQ=255)
                                                                tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGGTCGaaa                                       >  2:3763312/3‑101 (MQ=2)
                                                                cccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                         <  2:228415/99‑1 (MQ=2)
                                                                 gcggcGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      <  1:2006636/97‑1 (MQ=2)
                                                                 gcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      <  2:1907003/100‑1 (MQ=2)
                                                                 cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                     <  1:4081188/99‑2 (MQ=2)
                                                                 cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                     >  1:298260/3‑100 (MQ=2)
                                                                 aCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                     >  2:2773836/1‑101 (MQ=255)
                                                                 aCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                     <  1:3091156/101‑1 (MQ=255)
                                                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.ACCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                   <  1:3795696/101‑1 (MQ=255)
                                                                    tctgtCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccg                                   <  1:512214/97‑4 (MQ=2)
                                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcg                                   <  1:3408607/101‑4 (MQ=9)
                                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgc                                  >  1:456505/1‑97 (MQ=2)
                                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                  >  2:3791207/1‑97 (MQ=9)
                                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGAGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAActcgg                                  >  2:694300/1‑97 (MQ=9)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGTGGGCGTGGATCGAAACacc                                   >  1:2431591/1‑99 (MQ=2)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                >  2:1983309/1‑95 (MQ=2)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                >  1:168507/1‑95 (MQ=2)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                <  1:1729196/101‑7 (MQ=2)
                                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                 >  2:2381441/1‑96 (MQ=9)
                                                                      gCTCCCCCCGCGGGGGCGTGGAGCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GGCGCTCCCCCCGCGGGGGCGTGGATCGGAActacaac                                >  2:4167770/1‑95 (MQ=2)
                                                                      ccTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                <  2:2041455/100‑1 (MQ=255)
                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                  >  1:1117362/1‑101 (MQ=2)
                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTCGATCGAAACACCa                                  >  2:3707640/1‑101 (MQ=2)
                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGTTCGAAActaccccg                               >  1:4213785/1‑94 (MQ=2)
                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                <  2:2905661/101‑7 (MQ=9)
                                                                       cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                <  2:1071546/101‑7 (MQ=9)
                                                                        tCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcccc                               >  2:3383085/1‑94 (MQ=2)
                                                                        tCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccac                               >  1:957034/1‑94 (MQ=2)
                                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                              >  2:1396058/1‑101 (MQ=255)
                                                                         ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                <  1:2065473/101‑3 (MQ=2)
                                                                         ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                              <  2:3466092/101‑9 (MQ=9)
                                                                           ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                            >  1:3585290/1‑91 (MQ=2)
                                                                            cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                             <  1:1610809/101‑6 (MQ=0)
                                                                            cccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtcg                          <  2:3300753/101‑1 (MQ=255)
                                                                              cgcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                        <  1:1837398/101‑1 (MQ=255)
                                                                               gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                       <  1:588639/101‑1 (MQ=255)
                                                                                cgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.TGCA.TCGGCT.TCCCC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGATCGAAACATCAGCTCGTCGTTgg                      >  1:3950927/1‑101 (MQ=255)
                                                                                   gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                   <  1:2687169/101‑1 (MQ=255)
                                                                                   gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                   <  2:852037/101‑1 (MQ=255)
                                                                                     gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                 <  2:4080902/101‑1 (MQ=255)
                                                                                      cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGGCAGcc                >  1:1964785/1‑101 (MQ=255)
                                                                                        tGGATCGAAACA.GCCT.GCGCG.GG.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg              <  1:1917606/101‑1 (MQ=255)
                                                                                        tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg              >  1:1326538/1‑101 (MQ=255)
                                                                                           aTCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg           >  1:2527941/1‑101 (MQ=255)
                                                                                             cGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt         <  1:3648629/101‑1 (MQ=255)
                                                                                              gAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg        >  1:1780540/1‑101 (MQ=255)
                                                                                              gAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg        >  2:3908200/1‑101 (MQ=255)
                                                                                              gAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATg        >  2:1112564/1‑101 (MQ=255)
                                                                                                aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc      >  2:21219/1‑101 (MQ=255)
                                                                                                aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc      <  1:2265870/101‑1 (MQ=255)
                                                                                                   a.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt   <  2:3988085/101‑1 (MQ=255)
                                                                                                     |                                                                                                             
GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTC  >  NC_014034/1396836‑1397036

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