Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,997 G→C 14.7% T220T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9970G→C14.7% 240.5 0.0 0.0751T220T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (60/52):  ref base (340/296):  total (401/350)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 1760 total aligned reads displayed.

TGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGG  >  NC_014034/221897‑222097
                                                                                                    |                                                                                                    
ccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                                      <  1:4231042‑M1/42‑1 (MQ=255)
   ttaGCGGCGCCGGATCTCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                   <  2:910496/99‑1 (MQ=255)
   tcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtc                                                                                                   <  2:1909496‑M1/45‑1 (MQ=255)
   tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                   <  1:2458955/101‑1 (MQ=255)
    gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                  >  2:2933530‑M1/56‑101 (MQ=255)
        tcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtc                                                                                              <  2:824487‑M1/50‑4 (MQ=255)
         tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                                             >  1:3354438/1‑101 (MQ=255)
            gccggtccgtgggatcggtcgggtccgggggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                                                                          >  2:1239693‑M1/48‑101 (MQ=255)
              ggATCGCGCTCATCGACGGTGGTCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                                                        <  2:2223724/101‑1 (MQ=255)
                      gggatcggttggttccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                >  2:3959772‑M1/38‑101 (MQ=255)
                         atcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGggtccgtcggat                                                                             <  1:1070100‑M1/67‑3 (MQ=255)
                            ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                          <  1:3726309‑M1/70‑1 (MQ=255)
                                gggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                                      <  1:3105929‑M1/74‑1 (MQ=255)
                                 gatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTAAGTGGGa                                                                     <  2:3713965‑M1/75‑1 (MQ=255)
                                  gttcgtgggatcggtcgggtccgtgGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                    >  2:3289672‑M1/26‑101 (MQ=255)
                                   tccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                   >  1:3545069‑M1/25‑101 (MQ=255)
                                    ccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                  >  1:1229286‑M1/24‑101 (MQ=255)
                                    ccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                  >  1:59160‑M1/24‑101 (MQ=255)
                                           atcggtcgggtccggtGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGgatcggtcgggt                                                           >  1:3853090‑M1/17‑99 (MQ=255)
                                           atccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGAtcggtcgggg                                                           >  2:4063746/4‑99 (MQ=11)
                                            tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                          <  1:489648/101‑1 (MQ=255)
                                              ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGATCCg                                                        <  1:187516‑M1/88‑1 (MQ=255)
                                              ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCg                                                        <  1:538758‑M1/88‑1 (MQ=255)
                                              ggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                        <  1:2143941‑M1/88‑1 (MQ=255)
                                                tcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCgtc                                                      >  2:346192‑M1/12‑100 (MQ=255)
                                                 cgggtcggtgggGGCGGTCGGTTGCGTGGTGGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTgg                                                     <  1:2946009/97‑1 (MQ=11)
                                                 cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTgg                                                     <  1:2315‑M1/91‑1 (MQ=255)
                                                 cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGtccgggg                                                     <  2:335712‑M1/91‑4 (MQ=255)
                                                  tggggctggggatcggtcgggtccgtggggtcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTggg                                                    <  1:3796938/92‑1 (MQ=11)
                                                  gggtccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATGGGCCGGGTCCGGGGGATGAGGCGGGACCGCGGGGGCCGCCGGGTCCTTggg                                                    >  2:4300111/7‑101 (MQ=2)
                                                   ggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCAGAGCGGTCCGGCGGAGCCGTGGGa                                                   >  2:4640263/9‑101 (MQ=1)
                                                    ggccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                                  <  2:3019469/97‑1 (MQ=11)
                                                    cTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                  <  1:4140831/101‑1 (MQ=255)
                                                     tcggacggATCGGTCGGGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                 >  2:152885/1‑101 (MQ=255)
                                                     tccgtgggatcggtcgggtccgtgggatcggtcgggaccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGACCGTGGGGTCCGTCGGACCCTTGGGAtg                                                 >  2:1765501/7‑100 (MQ=11)
                                                      ccgggggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcg                                                <  2:3103215/96‑1 (MQ=11)
                                                        ggggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggt                                              <  2:3554467/92‑1 (MQ=14)
                                                         tgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                             <  2:1320868/99‑1 (MQ=11)
                                                         tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                             <  1:1308284‑M1/99‑1 (MQ=255)
                                                         tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                             >  1:441898‑M1/3‑101 (MQ=255)
                                                         tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                             >  2:10420‑M1/3‑101 (MQ=255)
                                                         tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                             >  1:199930‑M1/3‑101 (MQ=255)
                                                         tgGGATCGGTCGGGTCCCTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGgatccgtggggtccgtg                                             >  2:750522‑M1/3‑94 (MQ=255)
                                                         tgGGATCGGTCGGAGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                             >  2:3808875‑M1/3‑101 (MQ=255)
                                                         tagagtcggtcggatcgttgggatcggtcgtgttcgtggtgTCGGGCTGATCCTTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgttc                                             <  2:584171/96‑1 (MQ=2)
                                                         gggggttggtggggtccgtgggatcggtcgggtccgggggctcGGTCGGGTCTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                             <  1:1690243/94‑1 (MQ=14)
                                                         ggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                             >  1:2225514‑M1/3‑101 (MQ=255)
                                                          gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTAGGATCGGGCGGGTCGCGGGTAtcggtag                                            >  2:208437/2‑99 (MQ=14)
                                                          ggcagctgtcgtgtccgctggttcggtcgagtccgtggggtccGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                            <  2:3896714/94‑1 (MQ=1)
                                                          ggagtcggtcggtttcgtggggtccgtcgggtccgggggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                            <  2:1431435/97‑1 (MQ=14)
                                                          gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                            >  2:2843590‑M1/2‑101 (MQ=255)
                                                          gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                            >  1:2066170‑M1/2‑101 (MQ=255)
                                                               gggtcggggccgagggatcggtcgggtcggtggggtcGGTCGGGCCCGTGGGACCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                       <  1:175895/100‑1 (MQ=14)
                                                                gggcgggtccgtgggctcggtcgggtccgtgggatcggtCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                      <  1:2714668/98‑1 (MQ=11)
                                                                cgtcggatccgtgggatcggtcggatccgtgggatcgGTCGGATCCGTGGGAGCGGTCGGTTCCGTGGGTTCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                      >  1:1530813/2‑101 (MQ=11)
                                                                 gtcgggaccgtgaggtcggtccgatccgtcggatcgGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGt                                     <  2:2316954/101‑1 (MQ=11)
                                                                  tcgggtccgtcggatcggtcgggtccgtgggatcggTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtg                                    >  2:660527/1‑100 (MQ=11)
                                                                  gcgggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                    >  1:4371261‑M2/2‑100 (MQ=255)
                                                                   cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg                                   >  1:4502417‑M2/1‑99 (MQ=255)
                                                                    gggtccgtgggatcggtcgggcccgctgggtccgtcGGCGCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCCGTGCGGTGGGGTGGGCCGGGTGCGTGgg                                  >  2:2372894/1‑101 (MQ=2)
                                                                    ggatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                  >  1:2761660/4‑101 (MQ=11)
                                                                    ggatccgtgggatcggtcggatccgtgggatcggtcggTTCCGTGGGATCGGCCGGGTCCGTGCGGTCCGTCGGACCCTTGGGACCGGTCGGATCCGTGgg                                  >  1:794657/4‑101 (MQ=11)
                                                                    ggatccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                  <  2:2003864/98‑1 (MQ=11)
                                                                      atccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                <  2:456129/100‑1 (MQ=11)
                                                                        ccgtggggtcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                              <  1:3028825/101‑1 (MQ=11)
                                                                        ccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcg                              >  1:11739‑M2/1‑94 (MQ=255)
                                                                        ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtcg                              <  1:960087‑M2/101‑8 (MQ=255)
                                                                         cgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcgg                             <  1:3925535‑M2/101‑9 (MQ=255)
                                                                           tgggatcggtcgggtccgtggggtcgggcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtc                           >  2:3191617‑M2/1‑91 (MQ=255)
                                                                           tgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTGGGGTCCGTgggatcggtg                           <  1:2578626‑M2/101‑11 (MQ=255)
                                                                           tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtc                           <  2:2456202‑M2/101‑11 (MQ=255)
                                                                           ggggatcggtcggatcggggggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtcggtc                           <  1:1619032/100‑1 (MQ=11)
                                                                            gggatcggtcggttgcgggggctcggtcggTCCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTTCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                          <  1:320308/101‑1 (MQ=11)
                                                                            gggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGATCCGGGGGATCGGTCGGATCGGTGGGAtcggtcg                          >  2:2149933/1‑101 (MQ=11)
                                                                            gggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                          <  1:3465335‑M2/101‑12 (MQ=255)
                                                                             ggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgg                         <  1:4622361‑M2/101‑13 (MQ=255)
                                                                               atcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGCGGGATCGGTCGGATCCGTGGGACCGCGCGGGACGGGGGGAtcgggccggc                       >  2:1690132/1‑95 (MQ=2)
                                                                               atcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggt                       <  2:1599610‑M2/101‑15 (MQ=255)
                                                                                tcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatc                      <  2:2572114/101‑1 (MQ=11)
                                                                                tcggtcgggtccgtggaatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtc                      <  1:4075399‑M2/101‑16 (MQ=255)
                                                                                tccgtggggtcggtggggtccgtcggATCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                      <  2:2707158/95‑1 (MQ=11)
                                                                                   gtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGTGGGgtcgcgtccgt                   >  1:864002/1‑101 (MQ=11)
                                                                                     cgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgg                 <  2:2647590‑M2/101‑21 (MQ=255)
                                                                                     cgagtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTCGGCTCGGCCGGGTCCGTGGGCTCGgtcgggtccgtgg                 >  1:2655853/4‑99 (MQ=11)
                                                                                      gggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGGTAGGTCGGGCCCGTGGGATCGgtcgggtccgtcgg                >  1:3971322/1‑101 (MQ=11)
                                                                                      ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                >  2:3658674/4‑101 (MQ=255)
                                                                                       ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggg               <  2:4212928‑M2/101‑23 (MQ=255)
                                                                                       gatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATGGGTCGGGTCGGGGGGATCGgtcggggccggggga               >  2:2774787/3‑96 (MQ=11)
                                                                                       gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg               <  1:1003866/99‑1 (MQ=11)
                                                                                          ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGGCCGTCGGAACCggggggtccgtcgggtcc            >  2:1149379/1‑101 (MQ=11)
                                                                                           cgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtccg           <  2:518522‑M2/101‑27 (MQ=255)
                                                                                           cgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcgg           <  1:3294750‑M2/101‑27 (MQ=255)
                                                                                           cgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcgg           <  2:4128738‑M2/101‑27 (MQ=255)
                                                                                             tgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtc         <  1:237280‑M2/101‑29 (MQ=255)
                                                                                               ggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcgg       <  2:4382103‑M2/101‑31 (MQ=255)
                                                                                                   cggtcggGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc   >  1:688268‑M2/1‑67 (MQ=255)
                                                                                                   cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgcgggatcggtcgggtca   <  1:1431530‑M2/101‑35 (MQ=255)
                                                                                                    ggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccg  >  1:3658942‑M2/1‑66 (MQ=255)
                                                                                                    ggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTACGGCGGGTCGGTagggttcg  >  2:4449267/1‑96 (MQ=11)
                                                                                                    ggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccg  >  1:2860732/1‑99 (MQ=11)
                                                                                                    |                                                                                                    
TGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGG  >  NC_014034/221897‑222097

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