Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,036 C→G 4.1% P207P (CCG→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0360C→G 4.1% 163.2 ‑1.0 ‑0.02448P207P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (50/56):  ref base (997/1344):  total (1047/1401)
Fisher's exact test strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.81e-01

Only 100 of 3075 total aligned reads displayed.

TCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGG  >  NC_014034/221936‑222136
                                                                                                    |                                                                                                    
tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                                      <  2:1350910‑M1/81‑1 (MQ=255)
tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                                      >  1:1713355‑M1/21‑101 (MQ=255)
 gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                                     >  2:2917475‑M1/20‑101 (MQ=255)
  ggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                                    <  1:1241740‑M1/83‑1 (MQ=255)
  ggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                                                                                    >  2:3432664‑M1/19‑101 (MQ=255)
     tcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGggctcggtcgggtc                                                                                                 <  2:6395764‑M1/86‑4 (MQ=255)
      ccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                                <  1:6619995/100‑1 (MQ=1)
       ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                                               <  1:3125359‑M1/88‑1 (MQ=255)
        gtcgggtccgtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                                              >  2:1230682‑M1/13‑101 (MQ=255)
        gtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGTCCGt                                                                                              <  1:7101148‑M1/89‑1 (MQ=255)
          gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                                                            >  1:3666572/1‑99 (MQ=1)
            gccaggtgggATCGGTCGGGACCGGGGGCTCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGa                                                                                          <  2:3995474/97‑1 (MQ=11)
              tcggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc                                                                                        <  1:23257/101‑1 (MQ=255)
              tccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                        >  1:3239922‑M1/7‑101 (MQ=255)
                 gtggggttggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                     <  2:3060721/93‑1 (MQ=11)
                 gtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATAGGTCGGGACCGTGGGATCGGTCGGGACCGTGGGATCGGTCGGATCCGTGGGCtcggg                                                                                     >  1:1543636/4‑100 (MQ=11)
                 gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                     >  2:6796299‑M1/4‑101 (MQ=255)
                 gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                                                                                     <  2:2145006‑M1/98‑1 (MQ=255)
                 gtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGGCCGTGGGAtcggg                                                                                     >  2:1084712/1‑100 (MQ=11)
                  tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                    >  1:6163126‑M1/3‑101 (MQ=255)
                  tgGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGgatccgtggggtccgtg                                                                                    >  1:2772752‑M1/3‑94 (MQ=255)
                   gggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGGTCGGGCGGATCCGTGGGGTCGGGCGGGTCCGTGGGAtcggtcg                                                                                   >  2:4889884/2‑101 (MQ=11)
                    ggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                                                                  <  2:328398/101‑1 (MQ=11)
                      gtcggttgcgtgcgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcggat                                                                                <  2:708691/100‑3 (MQ=11)
                      gtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTGCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGTTCCGTGGCATCGGTCCGGt                                                                                >  1:4024203/5‑101 (MQ=11)
                      atcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGtccgtcggat                                                                                >  2:77425/1‑99 (MQ=11)
                      atcggtcgggtccgtgggatcggtcgggtccgtgggAACGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGTCCGTGGGGtcggtcggat                                                                                >  1:6321174/1‑99 (MQ=11)
                       tcggtcggttccgtgggatcggtcgggcccgtcggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                                                                               <  2:6543585/101‑1 (MQ=11)
                       tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATc                                                                               <  1:672426/101‑1 (MQ=11)
                        cggtctggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                              <  2:6253264/101‑1 (MQ=11)
                         ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGAGGGATCGGTCGGGTCCGTGGGGGCGGTCGGATCCg                                                                             >  2:2458152/1‑101 (MQ=11)
                         ggtcgggtccgtgggatcggtcgggtccgtcggatcGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGCGGAGCGGTCGGGTCCg                                                                             >  2:5461976/1‑101 (MQ=11)
                         ggtcggatccgtgggatcggtcgggtccgtgggatcGGTCGGATCCGTGGGATCGGCCGGACCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                             >  1:6944085/1‑101 (MQ=14)
                           gcgggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                                                           >  1:4182661‑M2/2‑100 (MQ=255)
                            cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGACCGCGGGATCGGTCGGGAGCGTGGGGTTGTTCGGGGCCGCGGGATCGGTCgggacagagg                                                                          >  2:5325254/1‑94 (MQ=2)
                             gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCggggg                                                                         >  1:1812361/1‑97 (MQ=11)
                             ggatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGgg                                                                         >  1:4709201/4‑101 (MQ=11)
                             ggatccgggggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCgg                                                                         <  1:2170725/98‑1 (MQ=11)
                              ggtccgtaagagcggtcgggtccgtggcatccgtcgGGCCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                                                                        <  2:517356/101‑1 (MQ=11)
                               ggcggtgggatcggtcggggccgtgggatcggtcggatccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggat                                                                       >  1:81807‑M2/5‑96 (MQ=255)
                                tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtc                                                                      >  1:2997548‑M2/1‑95 (MQ=255)
                                gccgtgggatcgggcggggccgggggatcggtcggATCCGTGGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                      <  1:3557759/100‑1 (MQ=14)
                                 ccgtgggttcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                                     <  2:4615009/101‑1 (MQ=11)
                                 ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATGGGTCGGGTCCGTGGGGtcg                                                                     >  2:4364882/1‑101 (MQ=11)
                                  cgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgggcggg                                                                    >  1:1000254/1‑96 (MQ=11)
                                  cggggggtcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                    <  1:1101313/98‑1 (MQ=11)
                                    tgtgacctgtctggcccgtggcgtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                  <  2:6929227/93‑1 (MQ=11)
                                    tgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtc                                                                  <  2:4272499‑M2/101‑11 (MQ=255)
                                    gggggtcgggcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                  <  1:6050515/100‑1 (MQ=11)
                                     gggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                                                                 <  2:5878310‑M2/101‑12 (MQ=255)
                                          cgggcggatccgggggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtcc                                                            >  1:1928472‑M2/1‑85 (MQ=255)
                                               gggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggg                                                       >  1:5126912‑M2/1‑80 (MQ=255)
                                               gggtccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgg                                                       <  1:4430781/101‑1 (MQ=14)
                                               ggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggg                                                       <  1:2438714‑M2/101‑22 (MQ=255)
                                               ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATGCGTGGGATCGGACGGGTTCGTGGGATCGgtcgggtccgtggg                                                       >  2:892953/4‑101 (MQ=14)
                                               ggatccgtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAGCGggcggggcggtcgg                                                       >  2:4722717/4‑101 (MQ=11)
                                               ggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgg                                                       >  1:2877763/4‑101 (MQ=11)
                                                   ccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcg                                                   <  2:6662700‑M2/101‑26 (MQ=255)
                                                   ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                   >  1:3912160/1‑101 (MQ=14)
                                                     gtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggt                                                 <  1:5001711‑M2/101‑28 (MQ=255)
                                                     gtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATACGTgggatcggtcggatccgtgggatcggt                                                 >  1:2659151‑M2/1‑74 (MQ=255)
                                                     gtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTcgggtccgtcggatccgt                                                 >  2:5755699/1‑101 (MQ=11)
                                                      tcggatcggtgggatCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtc                                                >  2:165045‑M2/3‑73 (MQ=255)
                                                         gatcggtcggGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGCGATCGGTGGGCTCGGTGGCGGCCGTGCGa                                             >  2:317666/1‑101 (MQ=11)
                                                         gatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGCGGCTCCGTgggg                                             >  1:1369259/1‑100 (MQ=11)
                                                           tcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGtc                                           <  1:953903/101‑1 (MQ=11)
                                                             ggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccg                                         <  1:5400807‑M2/101‑36 (MQ=255)
                                                              gtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgt                                        <  1:2048707‑M2/101‑37 (MQ=255)
                                                              gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTGGGCTCGGGGGGAtcggc                                        >  2:363943/1‑100 (MQ=11)
                                                              ggcgggggcGTGGGATCGGGCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtccgt                                        <  2:5017818/93‑1 (MQ=11)
                                                               tgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtc                                       <  1:5638645/99‑1 (MQ=17)
                                                                cggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtccgtgg                                      >  1:3815827/1‑99 (MQ=11)
                                                                  gATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGATCCGGAGGATCGGTAGGGTCCGTGGGATCGGTgggggcaggcgga                                    >  1:3328902/1‑91 (MQ=11)
                                                                     ccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGCGTTGCGTCGGAtcggaggagcgg                                 >  2:2695519/1‑93 (MQ=11)
                                                                       ttgGGATCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggt                               >  2:3895664‑M2/2‑56 (MQ=255)
                                                                       gTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGt                               <  2:4482800/101‑1 (MQ=14)
                                                                         aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGGTCGGTCGGGTCCGtcg                             >  2:1900233/2‑99 (MQ=14)
                                                                         aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGGCCGTCGGGTCCGTGGGAACGGTCGGGGCCGtcg                             >  2:3552135/2‑99 (MQ=14)
                                                                          ggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGGCCGTGGGGTCGGTCGGGTcgggcgg                            >  2:3780849/1‑95 (MQ=14)
                                                                              gggtcggGGCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                        <  1:5271922/100‑1 (MQ=17)
                                                                                gtcggGCCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgt                      >  2:2818270‑M2/1‑47 (MQ=255)
                                                                                  tgggtCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgg                    >  1:4863444‑M2/2‑45 (MQ=255)
                                                                                  cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                    >  2:4138050/1‑101 (MQ=17)
                                                                                  cggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgggg                    <  2:2589036‑M2/101‑57 (MQ=255)
                                                                                      tCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                >  2:2030335‑M2/1‑41 (MQ=255)
                                                                                       gcGTGGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc               <  2:558871/100‑1 (MQ=17)
                                                                                       ctgTGGGATGGGTCGGGTGCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTcc               <  1:5524831/99‑1 (MQ=11)
                                                                                       ccttGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCTTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc               <  2:5972404/98‑1 (MQ=17)
                                                                                       ccGTGGGATCGGTCGGGGCCGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc               <  1:6467928/101‑1 (MQ=11)
                                                                                       ccGTGGGAGCGGTGGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc               <  1:5771351/101‑1 (MQ=11)
                                                                                        cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCCGTCGGGTCCg              >  2:1803920/1‑101 (MQ=17)
                                                                                         gTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtcggtggggtcggg             <  2:6123759‑M2/101‑64 (MQ=255)
                                                                                              aTCGGTCGGATCCGTGGGATCGTTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat        <  2:1491287/101‑3 (MQ=14)
                                                                                               tcggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatc       <  1:4562104/101‑4 (MQ=14)
                                                                                                cggtcggGTCCGTGGGGTGGGAGGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcg      <  1:345504/101‑5 (MQ=14)
                                                                                                cggtcggGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatcc      >  1:2536474‑M2/1‑31 (MQ=255)
                                                                                                cggtcggGTCCGTGGGATCGGTCGGGTCCGTcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcc      >  1:6821819‑M2/1‑31 (MQ=255)
                                                                                                 tgtcgggTCCGTGGGGTCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccg     >  2:2328266‑M2/2‑30 (MQ=255)
                                                                                                  gtcggGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgt    <  1:4303218‑M2/101‑73 (MQ=255)
                                                                                                  gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggt    >  1:5901606/1‑95 (MQ=14)
                                                                                                    |                                                                                                    
TCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGG  >  NC_014034/221936‑222136

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