Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,051 G→C 1.4% T202T (ACC→ACG RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0510G→C 1.4% 23.5 ‑1.3 ‑0.51279T202T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/14):  ref base (545/711):  total (552/727)
Fisher's exact test strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.23e-01

Only 100 of 1924 total aligned reads displayed.

CGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGC  >  NC_014034/221951‑222151
                                                                                                    |                                                                                                    
 cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCGGGAtcgg                                                                                                     >  2:198955/5‑101 (MQ=11)
 cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                                                                                     <  1:3649362‑M1/97‑1 (MQ=255)
  gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTCGGATCGGGCGGGTCCGTGGGGTCGGGCGGATCCGTGGGAtcggt                                                                                                    >  2:1418832/4‑101 (MQ=14)
   tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                                                                   >  1:441898‑M1/3‑101 (MQ=255)
    gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGGCGGTCGGGTCCGTGGGGtcggtcg                                                                                                  >  2:4073333/2‑101 (MQ=11)
    gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                                  <  1:1177001‑M1/100‑1 (MQ=255)
    gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                                  >  2:3762018‑M1/2‑101 (MQ=255)
     ggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGAtcggacgg                                                                                                 >  2:3393771/1‑101 (MQ=14)
     ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                                                                                                 <  2:2690769/101‑1 (MQ=11)
     ggaacgggcgggtgcgtggggacggtaggatccgggGGGTCGGCCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgg                                                                                                 <  1:2089415/101‑1 (MQ=2)
      ggtcggtcgggtccgtgggatcggtcgggtccgtcggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCggg                                                                                                <  1:1599715/99‑1 (MQ=14)
       atcggtcggagccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcggat                                                                                               <  1:3906785/101‑3 (MQ=14)
        tccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTc                                                                                              >  2:533930/4‑101 (MQ=14)
         ccgtccgccccgtgcgatcggtcgggtccgtggcatcgggcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                                             <  1:4026774/92‑1 (MQ=11)
          ggtcgggtccgtgggatcggtcgggtccgggggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg                                                                                            <  1:3095396/101‑1 (MQ=14)
            tcggatccgtgggatcggtcgggtcggtgggatcggTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtg                                                                                          <  2:2064051/101‑2 (MQ=11)
            tcggatccgtgggatcggtcggatccgtgggatcggTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGCCGGGTCCGTg                                                                                          <  1:259414‑M2/101‑2 (MQ=255)
             cgggtccgtgggctcggtctggcgcgtcggctcggtCGGGTCCGTGTGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtgg                                                                                         <  2:4148905/101‑3 (MQ=11)
             cgggtccgtgggatcggtcgggtccgtcggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                                         >  1:3955710‑M2/1‑99 (MQ=255)
              gggttcgttggaatggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                                                        <  2:1595632/101‑1 (MQ=11)
               ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                                                                       <  1:763490‑M2/101‑5 (MQ=255)
               ggtccgtgggatcggtcgggtccgtgggatcggtggGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                                                                                       >  2:2452348‑M2/1‑97 (MQ=255)
                atccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgggt                                                                                      >  1:4160709/2‑99 (MQ=11)
                atccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTcccggggat                                                                                      >  1:93907/2‑94 (MQ=11)
                 tccgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatc                                                                                     >  2:130470‑M2/1‑95 (MQ=255)
                 tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                                                                                     >  2:1230328‑M2/1‑95 (MQ=255)
                 tccgggggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                     <  2:836456/101‑1 (MQ=11)
                  ccgtggggtgcgtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                    <  1:1771018/101‑1 (MQ=11)
                  ccgtgggatcggtcgggtccgtgggctccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcg                                                                                    >  2:3660328‑M2/1‑94 (MQ=255)
                  ccgggggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcg                                                                                    <  2:4308974/101‑1 (MQ=11)
                   cgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtccg                                                                                   <  2:2827918‑M2/101‑9 (MQ=255)
                     ttgggtcggtcggagccgtgggatcggtcggCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtcggtc                                                                                 <  2:4574783/96‑1 (MQ=11)
                     tgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtc                                                                                 <  1:1317153‑M2/101‑11 (MQ=255)
                     tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                 >  1:2667794/1‑101 (MQ=14)
                     tgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGGGGGAtgcggg                                                                                 >  2:2982680/1‑96 (MQ=11)
                      gggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGATCCGGGGGATCGGTCGGATCGGTGGGAtcggtcg                                                                                >  2:2149933/1‑101 (MQ=11)
                      gggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                                                                                <  1:4097315‑M2/101‑12 (MQ=255)
                      gggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                                                                                <  1:3816678‑M2/101‑12 (MQ=255)
                       ggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGAGGGATCGGACGGATCCGTTGGCtcggtggg                                                                               >  1:49218/1‑101 (MQ=11)
                             gtcgggtccgtgggttctgttgaATCCGTGGGATCCGTCAGATCCGTGCGATCGTTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                         <  2:3043192/101‑1 (MQ=11)
                              tcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtg                                                                        <  1:2243103‑M2/101‑20 (MQ=255)
                              gcgggtccgtcggatcggtgggATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                        <  1:2418842/100‑2 (MQ=11)
                               cgggtccgtgggatcgggcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcg                                                                       <  1:3283602/101‑1 (MQ=11)
                               cgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGggcgggtccgtcg                                                                       >  2:2539762/1‑101 (MQ=11)
                                gggtccgtgggatcggtcggGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggg                                                                      >  1:3202987‑M2/1‑80 (MQ=255)
                                ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGgtcgggtccgtcgg                                                                      >  1:2346802/4‑101 (MQ=11)
                                ggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtggggtccggcgg                                                                      >  2:152768/4‑101 (MQ=1)
                                 ggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggg                                                                     <  2:4212928‑M2/101‑23 (MQ=255)
                                 gatccgtaggatgggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcggg                                                                     <  1:1003866/99‑1 (MQ=11)
                                    ccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGgacgggtccgtaggatcg                                                                  >  2:2364124/1‑100 (MQ=11)
                                        gggctcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTgg                                                              <  1:201320/101‑1 (MQ=11)
                                         gggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTggg                                                             >  1:3554725/7‑101 (MQ=14)
                                          ggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGa                                                            <  2:2896088/93‑1 (MQ=14)
                                          ggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGCGGGa                                                            >  1:2499852/6‑101 (MQ=11)
                                             cggtcggGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc                                                         >  1:688268‑M2/1‑67 (MQ=255)
                                             cggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtcc                                                         <  2:3133904/101‑2 (MQ=11)
                                              ggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccg                                                        >  1:2860732/1‑99 (MQ=11)
                                               gtcggATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggcatcggtcgggtccgtgggatcggtcgggtccgt                                                       <  1:2796018‑M2/101‑37 (MQ=255)
                                               gtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgt                                                       <  1:124666‑M2/101‑37 (MQ=255)
                                                 cggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggacggggccgtgg                                                     <  2:1032007‑M2/101‑39 (MQ=255)
                                                 cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgtgg                                                     >  1:2933759/1‑99 (MQ=11)
                                                  ggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGGtcggccgg                                                    >  2:4001580/1‑101 (MQ=14)
                                                  ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtggg                                                    <  2:3203192‑M2/101‑40 (MQ=255)
                                                   ggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcggg                                                   >  1:553197/3‑100 (MQ=11)
                                                   gctccGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtggga                                                   >  2:3304456‑M2/3‑61 (MQ=255)
                                                     tCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATc                                                 <  1:1704658/101‑1 (MQ=11)
                                                        gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcggatccgtgggatcggt                                              <  1:4277303‑M2/101‑46 (MQ=255)
                                                         tcggATCCGTGGGGTCCGTGTGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                             <  2:2755818/99‑1 (MQ=18)
                                                             aTCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                         <  1:1929982/101‑1 (MQ=14)
                                                              gcggtcggGTCCGGCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtc                                        <  2:1910660/100‑1 (MQ=14)
                                                               cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                       >  1:217958‑M2/1‑49 (MQ=255)
                                                                ggtcggGTCCGTCGGATCGGTCGGATCCGTGGGCTCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGGGGGGTCGGTGCGGTCCGTCGGAgcgg                                      >  2:1674568/1‑101 (MQ=14)
                                                                ggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtcggtcgggtccgtgggatcggtcgggtccg                                      <  1:2725954‑M2/101‑54 (MQ=255)
                                                                cgtcggGTCCGTGGGCTCGGTCGGATCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                      <  2:39269/100‑1 (MQ=17)
                                                                  tcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                    >  1:85166/1‑101 (MQ=17)
                                                                  tcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                    >  1:3364341/1‑101 (MQ=17)
                                                                    gggTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGCCCGCCGGAtcggtccg                                  >  1:1075131/1‑99 (MQ=14)
                                                                    ggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtggg                                  <  2:73771‑M2/101‑58 (MQ=255)
                                                                      atcCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatccgtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                                >  2:2621567‑M2/2‑42 (MQ=255)
                                                                       tCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                               <  1:4356118/101‑1 (MQ=17)
                                                                       tCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatc                               >  1:3468240‑M2/1‑41 (MQ=255)
                                                                         cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCg                             <  1:1428376/101‑1 (MQ=17)
                                                                         cGGGGGTTCGGTCGGATCCGTGGGATCTGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                             >  2:613623‑M2/1‑39 (MQ=255)
                                                                          gtcggATCGGTCGGATCCGTGGGATCGGTGGGGTCCGTCGGATCCGTGGGGTCCGTCCGGTCCGGGGGATGCGTCGGGGCGGTCGGGTCGGTCGGGTcctg                            >  2:3860010/4‑99 (MQ=9)
                                                                           tGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtcggtcggatccgtgggatcggtcgggtccgtgggatcggtg                           <  1:4522559‑M2/101‑65 (MQ=255)
                                                                             ggATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                         >  2:4396997‑M2/1‑35 (MQ=255)
                                                                               aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggat                       >  1:1525413/1‑99 (MQ=17)
                                                                                tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTGGGGGGCGGgggggc                      >  2:3472029/1‑99 (MQ=11)
                                                                                 tggtcggaTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtcc                     >  2:2215950‑M2/2‑31 (MQ=255)
                                                                                       gggccGTGGGGTCGGGCGGGGCCGTGGGATCGGTCGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa               <  1:4595249/98‑1 (MQ=18)
                                                                                        gtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggat              >  2:3576817‑M2/2‑24 (MQ=255)
                                                                                           cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcgg           <  2:476560‑M2/101‑81 (MQ=255)
                                                                                           cGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcgg           <  2:2307456/101‑1 (MQ=255)
                                                                                              gggGTCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcg        <  1:3010780/101‑1 (MQ=255)
                                                                                               ggatcggtcgtgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgg       >  2:924451‑M2/1‑17 (MQ=255)
                                                                                               ggATCGGTCGGGTCCGTGGGATCCGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtcgg       <  1:3722096/101‑1 (MQ=255)
                                                                                                 atcgggtgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggt     >  2:4483297‑M2/1‑15 (MQ=255)
                                                                                                  tcggtcgggtccgtggggtcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtc    >  2:1775029‑M2/1‑14 (MQ=255)
                                                                                                   cggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcggatccgtgggatcgggcgggtcc   <  2:1142765‑M2/101‑89 (MQ=255)
                                                                                                    ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg  <  1:610824‑M2/101‑90 (MQ=255)
                                                                                                    |                                                                                                    
CGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGC  >  NC_014034/221951‑222151

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