Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,937 C→A 9.4% S273R (AGC→AGA RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9370C→A 9.4% 90.9 ‑0.3 ‑0.7465S273R (AGC→AGARCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/20):  ref base (209/178):  total (248/217)
Fisher's exact test strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.08e-01

Only 100 of 492 total aligned reads displayed.

TCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTC  >  NC_014034/1396837‑1397036
                                                                                                     |                                                                                                             
tCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.Gc                                                                                                               <  2:5352753/101‑1 (MQ=255)
tCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.Gc                                                                                                               <  2:49277/101‑1 (MQ=255)
 cGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.Gcc                                                                                                              <  1:5601343/101‑1 (MQ=255)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCt                                                                                                             >  2:1585550/1‑101 (MQ=255)
   cTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.g                                                                                                           <  2:3331613/101‑1 (MQ=255)
    tCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.gc                                                                                                          <  1:1680242/101‑1 (MQ=255)
     ccccccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.gcg                                                                                                         >  1:5357460/1‑101 (MQ=255)
     ccccccGCGGGGGCGTGGATCCAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.gcg                                                                                                         <  2:5738093/101‑1 (MQ=255)
      cccccGCGGGGGCGTGGATCGAAACCGGGCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GGCT.ccgc                                                                                                        >  1:154513/1‑101 (MQ=255)
              ggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.t                                                                                             >  1:1956333/1‑101 (MQ=255)
              ggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.t                                                                                             <  1:5162695/101‑1 (MQ=255)
                ggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCg                                                                                           >  1:5483944/1‑101 (MQ=255)
                 gCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCgg                                                                                          >  1:2401842/1‑101 (MQ=255)
                   gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt                                                                                        >  2:2107233/1‑101 (MQ=255)
                    tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt.t                                                                                      <  1:2292620/101‑1 (MQ=255)
                    tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt.t                                                                                      >  1:620574/1‑101 (MQ=255)
                      gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATNGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.Tc.c                                                                                   >  1:793522/1‑101 (MQ=255)
                           aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTg                                                                             >  2:5154856/1‑101 (MQ=255)
                           aaaCCGATCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTg                                                                             >  2:4398971/1‑101 (MQ=255)
                              ccGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGt                                                                        >  1:4852741/1‑101 (MQ=255)
                                  tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTc                                                                   <  1:2156084/101‑1 (MQ=255)
                                  tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTc                                                                   <  1:1122374/101‑1 (MQ=255)
                                   ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCg                                                                  >  1:3979871/1‑101 (MQ=255)
                                     agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCt                                                                <  1:868810/101‑1 (MQ=255)
                                      gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTc                                                               <  2:1490887/101‑1 (MQ=255)
                                       agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcc                                                              <  2:735996/101‑1 (MQ=255)
                                         aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcccc                                                            <  2:2318731/101‑1 (MQ=255)
                                           aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcccccc                                                          >  1:4276099/1‑101 (MQ=255)
                                              ccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GGCGCTCCCCCcgcg                                                       >  1:1269027/1‑101 (MQ=255)
                                                 cgcgGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCgggg                                                    <  1:1016808/101‑1 (MQ=255)
                                                 cgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCgggg                                                    >  1:268738/1‑101 (MQ=255)
                                                   cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGc                                                  <  2:4755575/101‑1 (MQ=255)
                                                     gACCTCCTCAGACGTCGCTACCCCCTCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCgg                                                >  2:178476/1‑100 (MQ=255)
                                                      aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                               <  1:1369689/101‑1 (MQ=255)
                                                       ctcgcggttgcGTCACTCCCCCCGCGGGGGCGTAGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTg                                               <  1:4994562/91‑1 (MQ=0)
                                                       cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                              <  2:620516/101‑1 (MQ=255)
                                                        tcgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgg                                              <  1:5265435/92‑1 (MQ=0)
                                                         gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCACTCCCCCCGCGGGGGCGTGGa                                             >  1:421215/10‑101 (MQ=9)
                                                          tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  2:1152388/93‑1 (MQ=9)
                                                          tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  2:4342947/93‑1 (MQ=9)
                                                            tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGGTCGa                                         >  1:1996142/1‑101 (MQ=255)
                                                            ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                          >  1:4526889/6‑101 (MQ=2)
                                                             cccccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  2:500830/97‑1 (MQ=2)
                                                             acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                         >  2:946426/6‑101 (MQ=9)
                                                              gcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                          >  1:2948433/4‑101 (MQ=2)
                                                              aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       >  2:683994/1‑101 (MQ=255)
                                                               gACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      <  1:411050/101‑1 (MQ=255)
                                                                cggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                     <  2:3559995/99‑2 (MQ=2)
                                                                ccGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  2:3803703/100‑1 (MQ=2)
                                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                      <  2:3060844/101‑1 (MQ=2)
                                                                  gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                   >  2:4454542/1‑101 (MQ=255)
                                                                  gTCCCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                   <  2:3895656/101‑1 (MQ=255)
                                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                     >  2:5553029/1‑101 (MQ=2)
                                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCC.GC.AGCG.AGGGCG.TCTGGC.ACA..GGAC.GTCTCTCCCCCCGCGGGGGCGTGGGTCGaaagaat                                   >  2:5623443/1‑97 (MQ=2)
                                                                   tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                  <  2:1704773/101‑1 (MQ=255)
                                                                   tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                  <  2:3614776/101‑1 (MQ=255)
                                                                   tCGCCCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActaca                                  <  2:1286268/101‑5 (MQ=2)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCCCTCCCCCCGCGGGGGCGGGGATCCAAAcccgc                                  >  1:4373381/1‑97 (MQ=2)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAg                                 >  2:2508000/1‑101 (MQ=255)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacac                                 >  2:4646957/1‑96 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                 <  1:3888855/101‑6 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                <  2:239843/101‑1 (MQ=255)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                >  2:2760915/1‑95 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                 >  1:1445904/1‑96 (MQ=9)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcca                                <  1:2217695/101‑7 (MQ=2)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                               >  1:3513749/1‑94 (MQ=2)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                <  1:4640327/101‑7 (MQ=9)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                 >  1:1842704/1‑100 (MQ=2)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                 >  2:3445870/1‑100 (MQ=2)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccac                               >  1:3759943/1‑94 (MQ=2)
                                                                       cccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacctg                              <  1:4776734/100‑9 (MQ=2)
                                                                       cccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                              <  2:4852853/100‑9 (MQ=2)
                                                                        ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                >  2:2795979/1‑99 (MQ=2)
                                                                        ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                              <  2:3240912/101‑9 (MQ=9)
                                                                        ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                              <  2:4283306/101‑9 (MQ=9)
                                                                         cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                               <  2:5061901/101‑4 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgt                             <  1:1639134/101‑10 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                             <  2:311322/101‑10 (MQ=9)
                                                                          ccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgtc                            <  2:785092/101‑11 (MQ=2)
                                                                          ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                            >  2:4225878/1‑91 (MQ=2)
                                                                            ccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCACCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccg                            <  1:36876/101‑7 (MQ=0)
                                                                            ccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGAATCGAAACATCAGCTcgtcgt                         <  2:2907424/101‑1 (MQ=255)
                                                                             cgcgGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccga                           >  2:1969327/1‑94 (MQ=0)
                                                                             cgcgGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgccga                           >  2:1374990/1‑94 (MQ=0)
                                                                             cgcgGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACACCaatgccga                           >  2:2243042/1‑94 (MQ=0)
                                                                                gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                     >  2:4581550/1‑101 (MQ=255)
                                                                                 ggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTc                    <  1:1753461/101‑1 (MQ=255)
                                                                                  gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                   >  2:3984824/1‑101 (MQ=255)
                                                                                    gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                 >  2:10536/1‑101 (MQ=255)
                                                                                    gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGAT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                 >  2:5228674/1‑101 (MQ=255)
                                                                                     cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                <  1:3035518/101‑1 (MQ=255)
                                                                                      gTGGATCGAAACA.GCCG.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg               <  1:1189524/101‑1 (MQ=255)
                                                                                        ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa             <  2:1862340/101‑1 (MQ=255)
                                                                                        ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCCCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTTGTTGGTGAACCggc             >  1:1799396/1‑100 (MQ=38)
                                                                                         gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATGGAAACATCAACTCGTCGTTGGTCAGCCGGAc            >  1:3614762/1‑101 (MQ=255)
                                                                                         gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc            <  2:871324/101‑1 (MQ=255)
                                                                                         gATCAAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc            <  2:4188171/101‑1 (MQ=255)
                                                                                            cGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt         >  1:2833236/1‑101 (MQ=255)
                                                                                                 cA.GCCT.GCGCG.GC.CGCA.TCAGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg    <  1:3265856/101‑1 (MQ=255)
                                                                                                    gCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc  <  2:3691992/101‑1 (MQ=255)
                                                                                                     |                                                                                                             
TCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTC  >  NC_014034/1396837‑1397036

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