Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,048 A→G 2.1% D203D (GAT→GAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0480A→G 2.1% 48.8 ‑3.8 ‑0.51616D203D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/30):  ref base (684/891):  total (694/922)
Fisher's exact test strand distribution p-value = 2.29e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.96e-01

Only 100 of 2126 total aligned reads displayed.

GCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAC  >  NC_014034/221948‑222148
                                                                                                     |                                                                                                    
  tcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                                    <  2:2454578/101‑1 (MQ=255)
  tccgtgggatcggtcgggtccgtgggatcggtcgggtccgt.gGGATCGGTCGGGTCCGTGGGATCCGTAGGGTCGGTGGGATCGGGCGGATCGGTGGGCtc                                                                                                    >  2:2884276/7‑101 (MQ=11)
   gcggggggtcggtcgggtccgtgggatcggtcggggccgt.cggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                                   <  1:670406/93‑1 (MQ=11)
    cgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                                  >  2:504117‑M1/5‑101 (MQ=255)
    cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGACGGTGGGAtcgt                                                                                                  >  2:477159/2‑100 (MQ=11)
      tggaagaggtcgggtccggcggatcggtcgggtccgt.cagatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtc                                                                                                <  2:2535777/94‑1 (MQ=11)
      tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                                                                >  2:3090242‑M1/3‑101 (MQ=255)
      tcggATGGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                <  2:608128/101‑1 (MQ=11)
       gGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGT.GGGATCTGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                               >  1:1399627‑M1/2‑101 (MQ=255)
       gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                               <  1:3498502‑M1/100‑1 (MQ=255)
         gatcggtcgggtccgtgggatcggtcgggtccgt.ggGATCGGTCGGTTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgga                                                                                             <  2:4418604/101‑2 (MQ=11)
          agcggtcgggcccgtgggttcggtggggggcgt.ggggtCGGTCGGGTCCCTGGGCTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcggat                                                                                            <  1:4714244/99‑3 (MQ=9)
           tcggtcgggtccgtgggatcggtcggggccgt.cgggTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                                                           <  2:341901/101‑1 (MQ=11)
           tcggtagggtccgggggatcggtcgggcccgt.cggaTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                                                                                           <  1:1467259/101‑1 (MQ=11)
            gggcgcgggccgtgggctcgggccggtccgt.gggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                                          <  2:4213338/92‑1 (MQ=11)
            cggtggggtccgtgggatcggtcgggtccgt.gggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                                                                          <  1:759189/101‑1 (MQ=14)
            cggtcgggtccgtgggatcggtcgggtacgt.gcggtCCGTCGGATCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                                          <  2:458473/101‑1 (MQ=14)
            cgctcgggtccgtgagatcggtcgggtccgt.cggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                          <  2:4141322/98‑1 (MQ=11)
            ccgtcgggtccgtgggatcggtcgggtccgt.cggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTcc                                                                                          <  2:1269246/99‑1 (MQ=14)
             ggtcggatccgtgggatcggtcggatccgt.aggatcGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                                         >  1:3000167/1‑101 (MQ=11)
              gtgggatccgggggatcggtcgggtccgt.gggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                                        <  2:4326875/95‑1 (MQ=11)
              gtcgggtccgtgggatcggtcgggtccgt.gggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCgg                                                                                        >  1:1060985/1‑100 (MQ=14)
               tcgggtccgtgggatcggtcgggtccgt.ggggtccgTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGATCCGTGGGACCGGTCGGGTCCgcg                                                                                       >  2:4683126/1‑99 (MQ=14)
               gcgtggccgtgggatcggtcggatccgg.gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTa                                                                                       <  2:665355/95‑1 (MQ=11)
                cgcgccagtgggatcggtcgggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtgg                                                                                      <  2:3498232/94‑3 (MQ=11)
                 gggtccgtgggatcggtcgggtccgt.ggtatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                     <  2:2198415/101‑1 (MQ=11)
                 gggtccgtgggatcggtcgggtccgt.gggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                                     <  1:3125997/101‑1 (MQ=11)
                 gggtccgtgggatcggtcgggtccgt.cggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGACGCGTCTGTGCGCTCGACCGGGCCCGtgga                                                                                     >  1:1964793/1‑98 (MQ=2)
                  ggttcgtgggatcggtcgtgtccgt.ggtattggttgGGTTCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                    <  1:3825487/101‑1 (MQ=14)
                   gtccgtgggatcggtcgggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggat                                                                                   >  1:4862697‑M2/1‑96 (MQ=255)
                    tccgtgggatgggtagggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatc                                                                                  >  1:3925406‑M2/1‑95 (MQ=255)
                    gccgtggggtcgttgggttccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                  <  1:2456107/100‑1 (MQ=14)
                     tcgtgggatcggtcgggtccgt.cggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                                                                 >  2:3767626‑M2/2‑94 (MQ=255)
                     ccgtggggtccgtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgggggagcg                                                                                 >  2:3687003/1‑93 (MQ=11)
                     ccgtgggatcggtcgggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                                                                 >  1:1873846‑M2/1‑94 (MQ=255)
                     ccgtgggatcggtcgcgtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                                                                 >  2:4473664‑M2/1‑94 (MQ=255)
                       gtgggatcggtcgggcccgt.gggatcgggcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGTTCGGTCGGGACGGTGGGAtccgt                                                                               >  1:3329344/1‑101 (MQ=11)
                       gtgggatcggtcggatccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggt                                                                               <  1:3650842‑M2/101‑10 (MQ=255)
                       gtgggatcggtcggatccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggt                                                                               <  1:1001112‑M2/101‑10 (MQ=255)
                       gtgggatcggtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGATCCGAGGGATCGGTCGGATCCGTAGGATCGGTCGCGTCCGTAGGAtcggt                                                                               >  1:4242424/1‑101 (MQ=11)
                        tgggatcggtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtc                                                                              <  1:2264172‑M2/101‑11 (MQ=255)
                        tgggatcggtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGTTCCGTGGGATCGGTCGGGACCGTGGGCtccgtc                                                                              >  1:4355637/1‑101 (MQ=11)
                         gggatcggtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGCCGGCtcggtcg                                                                             >  1:866401/1‑101 (MQ=11)
                            atcggtcgggtccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggg                                                                          >  1:4103212/1‑100 (MQ=11)
                            atcggtcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggt                                                                          <  1:2702138‑M2/101‑15 (MQ=255)
                             tctgtcggggccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                                                         <  2:900145/98‑1 (MQ=11)
                              gggtggggtccgtcggggtc.gtgggATCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc                                                                        <  1:998666/100‑1 (MQ=11)
                              cggtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtcggtcgggtcc                                                                        <  1:1329977‑M2/101‑17 (MQ=255)
                              cggtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatcc                                                                        >  1:547995‑M2/1‑85 (MQ=255)
                              cgggcgggtccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatcc                                                                        >  1:2890114‑M2/1‑85 (MQ=255)
                               ggtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccg                                                                       <  2:1913155‑M2/101‑18 (MQ=255)
                                gtcgggtccgt.gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgt                                                                      <  1:2537986‑M2/101‑19 (MQ=255)
                                 tcggatccgt.gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtc                                                                     <  2:2781320/101‑1 (MQ=255)
                                 gcggatccgt.aggagcggtctgGTCCGTAGGATCGGTCGGGTCCGTCGGATCTGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtc                                                                     <  1:2172192/100‑1 (MQ=11)
                                  cgggtccgt.cggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                                    >  1:4503252/1‑99 (MQ=11)
                                      tccgt.gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCTGTCGGGTCCGTgggatcggtcgggtccgtgggatc                                                                <  1:625403‑M2/101‑25 (MQ=255)
                                       ccgt.gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcg                                                               <  2:1679913‑M2/101‑26 (MQ=255)
                                         gt.ggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggt                                                             <  1:1547391‑M2/101‑28 (MQ=255)
                                          t.gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggtc                                                            <  2:2054555/101‑2 (MQ=11)
                                            ggggtccgtcggATCCGTGGGATCGGCCGGATCCGTGGGATCGGTGGGGTCCGTGGGGTCGGTGGGGTCCGTCGGCTCGgtcggctccgtggggcccccag                                                           >  2:2334648/1‑97 (MQ=2)
                                                 cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcg                                                      <  2:2124913/101‑1 (MQ=14)
                                                  ggtcggNTCCGTGGGATNGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccg                                                     >  2:384088‑M2/1‑66 (MQ=255)
                                                  ggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccg                                                     <  1:3408830‑M2/101‑36 (MQ=255)
                                                    tcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGCTCCGTGGGGtccgtg                                                   >  1:1146235/1‑100 (MQ=11)
                                                      ggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgtggg                                                 >  1:4787231/1‑101 (MQ=11)
                                                      ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtggg                                                 >  1:5006197/1‑101 (MQ=11)
                                                       ggtCCTTGGGATCTGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcggg                                                <  2:4034581/99‑2 (MQ=11)
                                                        aTCGGTCGGATCCGTGGGGCCCGTGGGATCGGTCGGGCCCGTGGGATAGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCGGTGGGTTCGGACGGac                                               >  2:1165207/1‑100 (MQ=14)
                                                         tCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATc                                              <  2:2609418/101‑1 (MQ=11)
                                                           cGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGCGGGCTCCGTGGGagcgg                                            >  2:1378799/1‑97 (MQ=11)
                                                             tGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc                                          <  1:4433456‑M2/101‑47 (MQ=255)
                                                               ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgg                                        <  1:4173701‑M2/101‑49 (MQ=255)
                                                                gATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCggcggg                                       >  2:1554230/1‑96 (MQ=17)
                                                                  tgggtcggGTCCGTGGGAGCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                     <  1:4921210/99‑1 (MQ=17)
                                                                   ccgtcggatCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatcc                                    >  1:2089273‑M2/3‑49 (MQ=255)
                                                                     gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcggatccgg                                  <  2:732996‑M2/101‑55 (MQ=255)
                                                                     gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgt                                  <  1:2356988‑M2/101‑55 (MQ=255)
                                                                     gtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgt                                  >  1:1249286‑M2/1‑47 (MQ=255)
                                                                     gtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgt                                  >  2:2701953‑M2/1‑47 (MQ=255)
                                                                       cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                <  1:920244‑M2/101‑57 (MQ=255)
                                                                         ggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtggggtcggtcggatccgtgggg                              <  1:26112‑M2/101‑59 (MQ=255)
                                                                         ggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggg                              <  1:4011809‑M2/101‑59 (MQ=255)
                                                                           tCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGACCGGTCGGGTCCGTCGGATCGGTCGGGTc                            >  2:1438281/1‑101 (MQ=14)
                                                                            ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcg                           <  2:3202267‑M2/101‑62 (MQ=255)
                                                                             cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcgg                          <  2:4862080‑M2/101‑63 (MQ=255)
                                                                             cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                          >  1:3213100/1‑101 (MQ=255)
                                                                                  gATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcgga                     >  2:1798682‑M2/1‑34 (MQ=255)
                                                                                  gATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcggg                     <  2:2978477‑M2/101‑68 (MQ=255)
                                                                                   aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                    >  1:3739933/1‑99 (MQ=18)
                                                                                      ggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATGCGTGGGATCGCTCGGGTCCGTGGGGCtgg                 >  2:837070/1‑101 (MQ=255)
                                                                                         cggATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg              <  2:3947552‑M2/101‑75 (MQ=255)
                                                                                          ggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcgg             >  1:1664345/1‑92 (MQ=14)
                                                                                           gATCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcggg            <  1:4619012/101‑11 (MQ=11)
                                                                                            gtcggtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat           >  1:3562003‑M2/2‑24 (MQ=255)
                                                                                             tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatc          >  1:2791101‑M2/1‑23 (MQ=255)
                                                                                             tccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc          >  1:1759190‑M2/1‑23 (MQ=255)
                                                                                              ccgggggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcc         >  2:1067565‑M2/1‑22 (MQ=255)
                                                                                                gtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt       <  2:3320413‑M2/101‑82 (MQ=255)
                                                                                                  gggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtgggatccgggggatcggtcgggtccgtgggatcggtcg     <  2:494256‑M2/101‑84 (MQ=255)
                                                                                                   ggATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcgg    >  2:3573102/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
GCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAC  >  NC_014034/221948‑222148

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 40 ≤  ATCG