Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,951 T→C 20.1% I278T (ATC→GCC)  RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9510T→C20.1% 137.6 ‑3.7 ‑0.0361I278T (ATC→ACC) RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (34/21):  ref base (117/146):  total (178/183)
Fisher's exact test strand distribution p-value = 2.54e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.19e-01

Only 100 of 377 total aligned reads displayed.

GGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGG  >  NC_014034/1396851‑1397051
                                                                                                         |                                                                                                          
     gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt                                                                                                       <  2:3407159/101‑1 (MQ=255)
      tGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCt.t                                                                                                     >  2:3586799/1‑101 (MQ=255)
            gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.Ct                                                                                             >  1:2122624/1‑101 (MQ=255)
            gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.Ct                                                                                             >  1:245020/1‑101 (MQ=255)
             aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTg                                                                                            >  2:1212239/1‑101 (MQ=255)
             aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTg                                                                                            >  2:2104088/1‑101 (MQ=255)
             aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CCC.CTg                                                                                            >  1:916672/1‑101 (MQ=255)
              aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..c                                                                                         >  1:2806305/1‑101 (MQ=255)
               aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..Cg                                                                                        <  2:3943762/101‑1 (MQ=255)
                 cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTg                                                                                      <  1:2962214/101‑1 (MQ=255)
                   gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGCGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGCg.gg                                                                                   >  1:3764269/1‑100 (MQ=255)
                    tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTc                                                                                  <  2:953218/101‑1 (MQ=255)
                      cAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGc                                                                                >  1:4271552/1‑101 (MQ=255)
                       agagATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCt                                                                               <  1:1781672/101‑1 (MQ=255)
                         agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcc                                                                             <  1:818834/101‑1 (MQ=255)
                         agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcc                                                                             >  1:4059258/1‑101 (MQ=255)
                         agaTAGACCGCGCGGACCTCCTCAGACGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcc                                                                             <  1:2754403/101‑1 (MQ=255)
                           aTAGACCGCGCGGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcccc                                                                           <  1:988952/101‑1 (MQ=255)
                            tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTccccc                                                                          >  1:3333695/1‑101 (MQ=255)
                            tAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTccccc                                                                          <  2:2610443/101‑1 (MQ=255)
                             aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTcccccc                                                                         <  1:1364585/101‑1 (MQ=255)
                              gACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCcg                                                                        <  2:244993/101‑1 (MQ=255)
                                ccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCcgcg                                                                      >  2:1613305/1‑101 (MQ=255)
                                   cgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCgggg                                                                   <  1:3551321/101‑1 (MQ=255)
                                    gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCggggg                                                                  >  1:1047280/1‑101 (MQ=255)
                                        aCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTg                                                              >  2:74058/1‑101 (MQ=255)
                                          ctcctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                            <  1:1580185/101‑1 (MQ=255)
                                          attggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  2:2386524/11‑101 (MQ=0)
                                           ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  2:3119270/92‑1 (MQ=2)
                                           gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  2:2125338/10‑101 (MQ=9)
                                           gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  1:3687396/10‑101 (MQ=9)
                                           gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                            >  1:100646/10‑101 (MQ=9)
                                           cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGa                                                            >  1:657/9‑101 (MQ=0)
                                            tggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  1:2274427/9‑101 (MQ=2)
                                            tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:3612837/93‑1 (MQ=9)
                                            acccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                             <  1:457664/94‑1 (MQ=0)
                                             ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  1:2994627/94‑1 (MQ=2)
                                             ctcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         <  2:1561235/101‑1 (MQ=255)
                                             caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                          >  1:3953331/8‑101 (MQ=9)
                                              tcAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        >  2:1725907/1‑101 (MQ=255)
                                              aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         >  1:1629128/7‑101 (MQ=9)
                                                cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  2:3522513/5‑101 (MQ=2)
                                                 tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGGGGGTCGaaa                                                      >  2:3763312/3‑101 (MQ=2)
                                                 gACGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                     <  1:1882133/101‑1 (MQ=255)
                                                 actgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  1:3241675/4‑101 (MQ=9)
                                                  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)
                                                  aCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                    >  1:4414090/1‑101 (MQ=255)
                                                  aCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                                    >  2:2773836/1‑101 (MQ=255)
                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                    >  1:4316482/2‑100 (MQ=9)
                                                   tgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAct                                                    <  2:1330382/100‑2 (MQ=9)
                                                   cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                   <  1:1512467/101‑1 (MQ=255)
                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.ACCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATc                                                  <  1:3795696/101‑1 (MQ=255)
                                                     tctgtCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccg                                                  <  1:512214/97‑4 (MQ=2)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                                <  2:3677026/101‑6 (MQ=2)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                                <  1:2622314/101‑6 (MQ=2)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TT.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGGGGGTCGAAACCTCAGc                                               >  1:2409259/1‑101 (MQ=37)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                               >  2:1767107/1‑101 (MQ=255)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCG.GCGCC.GC.AGCG.AGTGCG.TCTGGA.ACA..GGAT.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActatgc                                                <  1:308799/101‑6 (MQ=2)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                               >  1:168507/1‑95 (MQ=2)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                               <  1:1729196/101‑7 (MQ=2)
                                                       gCTCCCCCCGCGGGGGCGTGGAGCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GGCGCTCCCCCCGCGGGGGCGTGGATCGGAActacaac                                               >  2:4167770/1‑95 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                                 >  1:1117362/1‑101 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGTTCGAAActaccccg                                              >  1:4213785/1‑94 (MQ=2)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                               >  2:3795821/1‑95 (MQ=9)
                                                        cTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcc                                               <  2:2905661/101‑7 (MQ=9)
                                                         tCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgccac                                              >  1:957034/1‑94 (MQ=2)
                                                         tCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                             <  2:14609/101‑1 (MQ=255)
                                                         tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTCGATCGAAActcggcct                                              >  2:671292/1‑94 (MQ=9)
                                                         cccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                             <  1:3585781/100‑1 (MQ=255)
                                                           cccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgt                                           <  2:4308950/101‑1 (MQ=255)
                                                           cccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgggg                                           >  1:2316220/1‑91 (MQ=2)
                                                           cccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgggg                                           >  1:195512/1‑91 (MQ=2)
                                                           cccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcg                                            <  2:3233692/101‑10 (MQ=9)
                                                            ccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                          <  2:4443363/101‑1 (MQ=255)
                                                               cgcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                       <  1:1837398/101‑1 (MQ=255)
                                                               cgcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                       <  2:548894/101‑1 (MQ=255)
                                                                gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                      <  1:588639/101‑1 (MQ=255)
                                                                gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                      >  1:4106072/1‑101 (MQ=255)
                                                                gcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTg                                      >  2:2021740/1‑101 (MQ=255)
                                                                  gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCCCGGGGGCGTGGATCGAAACATCAGCTCGTCTTTGGt                                    >  2:4393506/1‑101 (MQ=255)
                                                                    gggGGTGGATGGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCGCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                  <  1:2083712/101‑1 (MQ=255)
                                                                    gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                  <  2:852037/101‑1 (MQ=255)
                                                                    gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                                  >  2:3873928/1‑101 (MQ=255)
                                                                      gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                                >  1:3916595/1‑101 (MQ=255)
                                                                      gCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGc                                <  2:4080902/101‑1 (MQ=255)
                                                                        gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                              <  2:3333458/101‑1 (MQ=255)
                                                                         tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                             >  1:1326538/1‑101 (MQ=255)
                                                                          ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                            <  2:613299/101‑1 (MQ=255)
                                                                             tCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGa                         >  2:1418349/1‑101 (MQ=255)
                                                                              cGAAACA.GCCT.GCCCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGAt                        <  1:3425929/101‑1 (MQ=255)
                                                                                 aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                     <  2:1101399/101‑1 (MQ=255)
                                                                                  accatctt.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                     <  2:3505768/94‑1 (MQ=255)
                                                                                    a.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                  <  2:3988085/101‑1 (MQ=255)
                                                                                         t.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt              >  2:2231599/1‑101 (MQ=255)
                                                                                          ggcgcG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt              <  1:3134398/100‑1 (MQ=255)
                                                                                           gcgcG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc             <  2:2402747/101‑1 (MQ=255)
                                                                                             gcg.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc           <  1:1186559/101‑1 (MQ=255)
                                                                                                    cGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgc      >  2:2324486/1‑101 (MQ=255)
                                                                                                         tCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg  <  1:139769/101‑1 (MQ=255)
                                                                                                         |                                                                                                          
GGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGG  >  NC_014034/1396851‑1397051

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