Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,961 C→T 12.3% R281R (CGC→TTT RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9610C→T12.3% 109.8 ‑0.7 ‑0.7333R281R (CGC→CGTRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (23/18):  ref base (143/148):  total (167/166)
Fisher's exact test strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.81e-01

Only 100 of 360 total aligned reads displayed.

TCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1396861‑1397061
                                                                                                           |                                                                                                        
  gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGTCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.Ct                                                                                                       <  2:2383937/101‑1 (MQ=255)
  gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.Ct                                                                                                       <  1:1094541/101‑1 (MQ=255)
  gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.Ct                                                                                                       <  1:1872164/101‑1 (MQ=255)
  gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.Ct                                                                                                       <  2:2247220/101‑1 (MQ=255)
    aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..c                                                                                                   <  2:1532053/101‑1 (MQ=255)
    aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..c                                                                                                   <  2:1890588/101‑1 (MQ=255)
     aCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..Cg                                                                                                  >  2:162702/1‑101 (MQ=255)
        ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTg.g                                                                                              <  1:2691227/101‑1 (MQ=255)
        ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTg.g                                                                                              >  2:1955516/1‑101 (MQ=255)
         gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCGCGCGGGGGCGTGGATCGAAGCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CCG..CGGg.gc                                                                                             >  1:3985993/1‑100 (MQ=38)
         gTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.Gt                                                                                             <  1:2045554/101‑1 (MQ=255)
          tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTc                                                                                            <  1:2797972/101‑1 (MQ=255)
          tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTc                                                                                            <  1:1926985/101‑1 (MQ=255)
          tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTc                                                                                            <  2:463615/101‑1 (MQ=255)
           ccAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCg                                                                                           <  1:3079118/101‑1 (MQ=255)
            cAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGc                                                                                          >  2:3595621/1‑101 (MQ=255)
              gagaTAGACCGCGCGGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTc                                                                                        <  1:1585332/101‑1 (MQ=255)
              gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTc                                                                                        >  2:2317944/1‑101 (MQ=255)
                gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTccc                                                                                      >  2:4221109/1‑101 (MQ=255)
                   aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCCcccccc                                                                                   >  2:210502/1‑101 (MQ=255)
                    gtccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCACCCcg                                                                                  <  2:944517/99‑1 (MQ=255)
                      ccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCCCCCCCcgcg                                                                                >  1:351994/1‑101 (MQ=255)
                        gcgcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCggg                                                                              >  1:1477450/1‑101 (MQ=255)
                          gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCACTCCCCCCGCggggg                                                                            >  1:2714616/1‑101 (MQ=255)
                           cgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGc                                                                           <  1:1034042/101‑1 (MQ=255)
                                 ttggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     >  1:1786666/10‑101 (MQ=2)
                                 gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      <  2:3143444/92‑1 (MQ=9)
                                 gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      <  1:477244/92‑1 (MQ=9)
                                 cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      >  2:4593709/9‑101 (MQ=2)
                                 cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                      <  2:2701087/93‑1 (MQ=2)
                                  tcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     <  1:4507913/93‑1 (MQ=9)
                                  gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTgggt                                                                     >  2:4362149/8‑99 (MQ=0)
                                  gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                                     <  2:826403/94‑1 (MQ=2)
                                   ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   >  2:3939412/8‑101 (MQ=2)
                                    aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                   <  1:4392637/95‑1 (MQ=9)
                                      ggtggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCC.GC.AGCG.AGGGCG.TCTGGA.ACA..GGAT.GTCGCTCCCCCCGCCGGGGCGTGGAtggga                                                                 >  1:3690788/5‑97 (MQ=0)
                                      cccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                <  1:3159427/97‑1 (MQ=2)
                                      aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                >  2:1756737/1‑101 (MQ=255)
                                       tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGGGTGGATCGaaa                                                                >  2:1783993/3‑101 (MQ=2)
                                         cGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAt                                                             <  2:563786/101‑1 (MQ=255)
                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                                             <  1:2917010/101‑3 (MQ=9)
                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGCATCGAAACATCa                                                           <  1:2521070/101‑1 (MQ=255)
                                           ccccTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCa                                                           <  1:1659474/98‑1 (MQ=255)
                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACac                                                             >  1:2232993/1‑100 (MQ=2)
                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcgg                                                           >  1:572720/1‑97 (MQ=9)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACacc                                                            <  1:1346730/101‑3 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                                          <  2:346236/101‑6 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGc                                                         >  1:97918/1‑101 (MQ=255)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                         <  2:3673995/101‑7 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                                         <  2:1113765/101‑7 (MQ=2)
                                             gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                                          >  1:3204648/1‑96 (MQ=9)
                                              cTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaaaccgcca                                                         >  2:1186824/1‑94 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                                          <  2:3521576/101‑2 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                                       >  2:2967744/1‑93 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TA.TGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                                       >  2:2429733/1‑93 (MQ=2)
                                               tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                                        >  1:2141991/1‑94 (MQ=9)
                                               gccccccGCGGGGGCGTGGATCGAACCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                       <  1:210526/100‑1 (MQ=255)
                                               cccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                                       <  1:3017310/100‑1 (MQ=255)
                                                ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                                         <  1:247915/101‑3 (MQ=2)
                                                ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCTCCCCCCGGGGGGGCGTGGATCGAAAcccgccacg                                                       >  1:3922382/1‑93 (MQ=0)
                                                ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GGCGCTCCCCCCGCGGGGGCGTGGATCGTAActcggcgtc                                                       >  2:3349999/1‑93 (MQ=9)
                                                 cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                                        <  1:920291/101‑4 (MQ=2)
                                                  cccgGCGGGGGCGTGGATCGAAACATGGCTAGGGCG.GC.AGTG.TCAGGT.T..TTC.CAG..CGGC.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgtc                                                     <  2:3250934/101‑11 (MQ=2)
                                                  ccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgc                                                       <  2:2927158/101‑5 (MQ=2)
                                                  ccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCTGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                                    >  1:855830/1‑101 (MQ=255)
                                                  ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTT.TGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                                     <  2:4425292/101‑11 (MQ=2)
                                                   cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AG.CGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                                      <  2:1212956/101‑6 (MQ=0)
                                                     cgggggggCGTGCATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                                 <  2:1426817/98‑1 (MQ=255)
                                                     cgcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                                 >  2:2638554/1‑101 (MQ=255)
                                                     cgcgGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGtt                                                 >  1:378779/1‑101 (MQ=255)
                                                       ggggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                               <  2:2189925/100‑1 (MQ=255)
                                                       cgGGGGGGTGGATCGAAACA.GCCT.GCCCG.GN.CGCA.TCGGCT.TC.CTC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTgg                                               <  1:1231237/101‑1 (MQ=255)
                                                              gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                        >  1:4241087/1‑101 (MQ=255)
                                                              gTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCg                                        >  1:2625043/1‑101 (MQ=255)
                                                                ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                      <  2:3776703/101‑1 (MQ=255)
                                                                ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                                      >  2:1046852/1‑101 (MQ=255)
                                                                 gATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGAc                                     <  1:3735938/101‑1 (MQ=255)
                                                                  aTCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                                    >  2:4014005/1‑101 (MQ=255)
                                                                      aaaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGt                                >  1:3480855/1‑101 (MQ=255)
                                                                       aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc                               >  2:132336/1‑101 (MQ=255)
                                                                         cA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCg                             <  2:1251902/101‑1 (MQ=255)
                                                                          a.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt                            <  2:2037908/101‑1 (MQ=255)
                                                                              cT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGGTCGTCGTTGGTCAGCCGGACGATGTCCGTCGc                         >  2:2538739/1‑101 (MQ=255)
                                                                              cT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGGCAGCCGGACGAGGTCCGTCGc                         >  2:414347/1‑101 (MQ=255)
                                                                               t.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCt                        >  1:809899/1‑101 (MQ=255)
                                                                                  cgcg.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcc                      >  2:701828/1‑101 (MQ=255)
                                                                                   gcg.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc                     <  2:2483171/101‑1 (MQ=255)
                                                                                    cg.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccc                    <  1:2018958/101‑1 (MQ=255)
                                                                                       gC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc                  <  1:2289942/101‑1 (MQ=255)
                                                                                       gC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcccccc                  <  2:4511398/101‑1 (MQ=255)
                                                                                           gCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCcgcg               >  1:2411195/1‑101 (MQ=255)
                                                                                            cA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgg              >  1:652429/1‑101 (MQ=255)
                                                                                             a.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggg             >  2:3046191/1‑101 (MQ=255)
                                                                                             a.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggg             >  1:2390105/1‑101 (MQ=255)
                                                                                               tCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCgggg            <  1:4161544/101‑1 (MQ=255)
                                                                                                cgccT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCCGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCggggg           <  2:1688783/98‑1 (MQ=255)
                                                                                                  gCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCg         <  2:855754/101‑1 (MQ=255)
                                                                                                    t.tC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCggg       >  2:1386681/1‑99 (MQ=255)
                                                                                                    gtttgcGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCg         >  1:4317809/6‑101 (MQ=255)
                                                                                                      tC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGGgg      >  1:4574773/1‑101 (MQ=255)
                                                                                                           |                                                                                                        
TCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TC.CGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCCCGCGGGGGCGTGGATCG  >  NC_014034/1396861‑1397061

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG