Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,396,944 +C 12.3% coding (826/1578 nt) RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,396,9441.→C12.3% 130.6 ‑0.7 0.0332coding (826/1578 nt)RCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (23/18):  ref base (146/145):  total (169/163)
Fisher's exact test strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 364 total aligned reads displayed.

CCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCC  >  NC_014034/1396844‑1397044
                                                                                                       |                                                                                                           
ccccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.Gcgcg                                                                                                              <  1:1006551/101‑1 (MQ=255)
  ccGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACA.GCCT.GGGCA.Gc                                                                                                           >  2:790807/1‑101 (MQ=255)
      gggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCa                                                                                                      >  1:3603195/1‑101 (MQ=255)
      gggggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCa                                                                                                      >  1:1083970/1‑101 (MQ=255)
         ggCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCg                                                                                                  >  2:4602678/1‑101 (MQ=255)
            gTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCC.TCGGCt                                                                                               >  1:2112765/1‑101 (MQ=255)
               gATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.Tcc                                                                                           <  2:2558992/101‑1 (MQ=255)
                aTCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCg                                                                                          <  1:4016466/101‑1 (MQ=255)
                   gAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.Ct                                                                                      <  2:2247220/101‑1 (MQ=255)
                    aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTg                                                                                     >  1:1712842/1‑101 (MQ=255)
                    aaaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCCC.CTg                                                                                     >  1:1872347/1‑101 (MQ=255)
                     aaCCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..c                                                                                  <  2:1890588/101‑1 (MQ=255)
                        cGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..Cggg                                                                               >  1:627415/1‑99 (MQ=255)
                         ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg.g                                                                             >  2:1955516/1‑101 (MQ=255)
                         ggTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTg.g                                                                             <  1:2691227/101‑1 (MQ=255)
                           tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTc                                                                           <  1:1926985/101‑1 (MQ=255)
                           tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTc                                                                           <  2:463615/101‑1 (MQ=255)
                           tCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTc                                                                           <  1:2797972/101‑1 (MQ=255)
                               gagaTAGACCGCGCGGCCCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                       <  1:1585332/101‑1 (MQ=255)
                               gagaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTc                                                                       <  1:527195/101‑1 (MQ=255)
                                agaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CCG..CGTg.gggcctcc                                                                      >  1:990866/1‑94 (MQ=255)
                                 gaTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTccc                                                                     >  1:1359517/1‑101 (MQ=255)
                                  aTAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccc                                                                    <  2:4145934/101‑1 (MQ=255)
                                    aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTcccccc                                                                  >  1:2023556/1‑101 (MQ=255)
                                    aGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCCcccccc                                                                  >  2:210502/1‑101 (MQ=255)
                                     gtccGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCACCCcg                                                                 <  2:944517/99‑1 (MQ=255)
                                           gcgGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCACTCCCCCCGCggggg                                                           >  1:2714616/1‑101 (MQ=255)
                                              gACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGt                                                        >  2:2728398/1‑101 (MQ=255)
                                                cctcctCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTgg                                                      <  2:841411/101‑1 (MQ=255)
                                                 ggtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTgg                                                      >  2:2164731/11‑101 (MQ=9)
                                                  gtcaacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGa                                                     <  2:3607657/92‑1 (MQ=9)
                                                  cgcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                     >  1:3768479/9‑101 (MQ=2)
                                                   tgccccccGTCCCTCCCCGCGCGGGGGCGTGGATCGAAGCA..CCT.GTCC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                   <  2:3133455/94‑1 (MQ=2)
                                                   gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                    <  2:685449/94‑1 (MQ=2)
                                                   gcggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                                    <  2:24096/94‑1 (MQ=2)
                                                   acccgcccGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTgg                                                      <  2:2478642/94‑1 (MQ=0)
                                                    ggccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  >  2:3939412/8‑101 (MQ=2)
                                                    caacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGGTc                                                   >  1:4261340/8‑101 (MQ=9)
                                                     ggttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGGTCg                                                  >  2:439236/6‑101 (MQ=2)
                                                     gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                 <  1:1320863/95‑1 (MQ=2)
                                                     gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                 <  2:2904922/95‑1 (MQ=2)
                                                     gccccggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                 <  1:3724628/95‑1 (MQ=2)
                                                     aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  >  2:1455556/7‑101 (MQ=9)
                                                     aacactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  >  1:3026687/7‑101 (MQ=9)
                                                      gttgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                 <  1:3430769/97‑1 (MQ=2)
                                                      acactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                 >  1:782169/6‑101 (MQ=9)
                                                       cactgTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                >  1:1808681/5‑101 (MQ=9)
                                                       aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               >  2:1756737/1‑101 (MQ=255)
                                                       aGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               <  2:4192446/101‑1 (MQ=255)
                                                        tgggTCGCCCCCCCGGCGGGGGCGTGGATCGAAACATGTCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               <  2:2528082/98‑1 (MQ=2)
                                                        tgcGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGGGTGGATCGaaa                                               >  2:1783993/3‑101 (MQ=2)
                                                          ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActa                                            >  2:112522/2‑99 (MQ=2)
                                                           gTCGCTCCCCTCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAccc                                            >  1:149586/1‑99 (MQ=2)
                                                           gTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActc                                            >  2:2680248/1‑99 (MQ=9)
                                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                             >  2:1403036/1‑101 (MQ=2)
                                                            tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACa                                             <  1:4358326/101‑1 (MQ=2)
                                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACac                                            >  2:4065884/1‑100 (MQ=2)
                                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACac                                            >  1:2232993/1‑100 (MQ=2)
                                                             cGCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacac                                         >  1:3830550/1‑96 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                         <  2:346236/101‑6 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcccgcc                                         <  1:2911496/101‑6 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacacc                                        <  2:523761/101‑7 (MQ=2)
                                                              gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggc                                         >  1:3204648/1‑96 (MQ=9)
                                                               ccccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccg                                       <  2:4447132/99‑8 (MQ=2)
                                                               cTCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCa                                          >  1:2046855/1‑101 (MQ=2)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaa                                         <  2:3521576/101‑2 (MQ=2)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                      >  2:2967744/1‑93 (MQ=2)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccgg                                      >  1:2398381/1‑93 (MQ=2)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                       >  1:2141991/1‑94 (MQ=9)
                                                                tCCCCCCGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcct                                       >  1:2581699/1‑94 (MQ=9)
                                                                cccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                      <  1:4602979/100‑1 (MQ=255)
                                                                cccccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtc                                      <  2:1188477/100‑1 (MQ=255)
                                                                 ccccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACACCaat                                        <  1:247915/101‑3 (MQ=2)
                                                                 ccccccGCGGGGGCGTGGATCGAAACATGGCTATGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCTCCCCCCCCCGGGGCGTGGATCGGAAcccgccacg                                      >  2:1574265/1‑93 (MQ=0)
                                                                 ccccccGCGGGGGCGTGGATCGAAACA..ACT.GACC...C.GGCA.ACGGCC.TATGCGCTGGACGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActacaccggg                                     <  1:4485809/101‑10 (MQ=2)
                                                                 ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctc                                      >  1:2991329/1‑93 (MQ=9)
                                                                 ccccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GGCGCTCCCCCCGCGGGGGCGTGGATCGTAActcggcgtc                                      >  2:3349999/1‑93 (MQ=9)
                                                                  cccccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatg                                       <  1:920291/101‑4 (MQ=2)
                                                                   cccgGCGGGGGCGTGGATCGAAACATGGCTAGGGCG.GC.AGTG.TCAGGT.T.TTC.CAG..CGGC.GTCGCCCCCCCCGCGGGGGCGTGGATCGAAAcccgccacgtc                                    <  2:3250934/101‑11 (MQ=2)
                                                                   ccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCTGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                   >  1:855830/1‑101 (MQ=255)
                                                                   ccccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCtcgtc                                   >  1:2293306/1‑101 (MQ=255)
                                                                   ccccGCGGGGGCGTGGATCGAAAC..GCTG.GCGCGCGC.CGCG.CCGGATCTTTGT.CTG..CGCG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAActcggcctcga                                    >  1:2507218/1‑91 (MQ=2)
                                                                    cccGCGGGGGCGTGGATCGAAACT.GCGT.GCGCG.CCATGAATTCATCG.AGCGC.CAT..CGAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCaatgcc                                     <  2:1212956/101‑6 (MQ=0)
                                                                     ccGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTcgtcgt                                 >  1:4200435/1‑101 (MQ=255)
                                                                         gggggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGCGGGCGTGGATCGAAACATCAGCTCGTCGTTGGt                             >  1:3965679/1‑101 (MQ=255)
                                                                           gggCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCa                           <  2:2369390/101‑1 (MQ=255)
                                                                              cGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGcc                        >  2:4261865/1‑101 (MQ=255)
                                                                                tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                      >  2:3703090/1‑101 (MQ=255)
                                                                                tGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCgg                      <  2:3195365/101‑1 (MQ=255)
                                                                                 ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                     <  2:3776703/101‑1 (MQ=255)
                                                                                 ggATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGa                     >  2:1046852/1‑101 (MQ=255)
                                                                                   aTCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACg                   >  2:4014005/1‑101 (MQ=255)
                                                                                        aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc              >  2:132336/1‑101 (MQ=255)
                                                                                        aaCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTc              >  2:3023395/1‑101 (MQ=255)
                                                                                         aCA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTcc             <  1:3150124/101‑1 (MQ=255)
                                                                                           a.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGt           <  2:2037908/101‑1 (MQ=255)
                                                                                             gCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTc          >  1:3512471/1‑101 (MQ=255)
                                                                                                  gcgcG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTc      <  2:1653807/101‑1 (MQ=255)
                                                                                                   cgcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTcc     >  1:1182362/1‑101 (MQ=255)
                                                                                                    gcg.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTccc    >  2:1256528/1‑101 (MQ=255)
                                                                                                       |                                                                                                           
CCCCGCGGGGGCGTGGATCGAAACCGGTCCAGAGATAGACCGCGCGGACCTCCTCAGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GCCT.GCGCG.GC.CGCA.TCGGCT.TCCGC.CTG..CGTG.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCAGCTCGTCGTTGGTCAGCCGGACGATGTCCGTCGCTCCCCC  >  NC_014034/1396844‑1397044

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