Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,063 A→G 100% P198P (CCT→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0630A→G100.0% inf 837P198P (CCT→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (171/258):  ref base (123/97):  total (383/454)
Fisher's exact test strand distribution p-value = 1.23e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.96e-02

Only 100 of 1072 total aligned reads displayed.

CGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGC  >  NC_014034/221964‑222163
                                                                                                   |                                                                                                     
 gggtccgtgggatcggttgggtccgtggggtccgtcGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgg                                                                                                     <  1:4804519/101‑1 (MQ=11)
 ggatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCGGTGGGGTCCGTTgg                                                                                                     >  2:3651599/4‑101 (MQ=11)
  ggtccgtgggatcggtcgggtccgtcggatcggtggGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGa                                                                                                    >  2:4804714/1‑101 (MQ=14)
   tctcgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                   <  2:4665892/98‑1 (MQ=11)
   gaccgtgtgatccgtcgggtccgtgggaccggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                                   <  2:5394203/99‑1 (MQ=11)
     ccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGGtcg                                                                                                 >  2:1591751/1‑101 (MQ=11)
     ccgtgggatcggtcggatccgtgggagcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGAtcg                                                                                                 <  2:899862/101‑1 (MQ=11)
       gtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                               <  2:1976276/101‑1 (MQ=11)
        tgcggtcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                              <  2:1682355/96‑1 (MQ=11)
         gggatcggttggatccgggggatcggtcggGTCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                             <  1:5237514/101‑1 (MQ=11)
         gggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                             >  2:2714051/1‑101 (MQ=14)
         gggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                                                             <  1:2474641/101‑1 (MQ=11)
             tcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                                                                         <  2:4926860/101‑1 (MQ=11)
              cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGGCCGTGGGATCGgtcgggtcc                                                                                        >  1:1125334/1‑101 (MQ=11)
              cgggcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatcc                                                                                        <  1:1362309/101‑1 (MQ=11)
               ggtggggtccgggggatcgggcggAGCCGTGGGATCGGTCGGGTCCGTGGGATTGGTCGGATCCGTGGGAGCGGTCGGATCCGTGGGATCGgtcggatccg                                                                                       <  2:2076184/101‑1 (MQ=11)
               ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGTTCGGTCGGGTCCGTGGGCTCGGTCGGGTCCGTAGGCtcggtcggtaccg                                                                                       >  1:654560/1‑96 (MQ=11)
               ggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTGGGATCGGTGGGGTCCGTGGGGTCCgtcgggtcgg                                                                                       >  2:3066645/1‑99 (MQ=11)
               cgtcggttccgtgggatcggtcggATCCGTTGTATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCgtcggatccg                                                                                       <  2:1217750/100‑1 (MQ=11)
                gtcgggtccgtggggtccgtcgtATCCGTGGGACCGGTCGGATCTGTGGGATCGGTCGGATCCTTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgt                                                                                      <  1:40657/101‑1 (MQ=11)
                 tcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtg                                                                                     >  1:76083/1‑100 (MQ=11)
                 tcgggtccgtcggatcggtgggATCCGTGGGATCGGTAGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtccggg                                                                                     >  2:3125710/1‑99 (MQ=11)
                 tcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGTGACGGTGGGATGGgtcggatccgtg                                                                                     >  1:1703413/1‑100 (MQ=11)
                  cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                                                    <  1:3730649/101‑3 (MQ=11)
                  cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                                                    >  2:1283548/1‑99 (MQ=11)
                  cggatccgggggatcgggcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcg                                                                                    <  1:303885/101‑1 (MQ=14)
                  ccggtcgttgggatgtgtcggCTCCGTGTGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                                                    <  2:2611623/93‑3 (MQ=1)
                   ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                   >  2:3810101/4‑101 (MQ=255)
                   ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                                   <  2:5029662/98‑1 (MQ=255)
                    ggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGGCCGTGGGATCGgtcgggtccgtggga                                                                                  >  2:4548153/1‑100 (MQ=11)
                    ggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGGGGGATCGggcgggtccgtggga                                                                                  >  2:836836/1‑100 (MQ=11)
                    gatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTCGGGCGGgcggggtccgtcgga                                                                                  >  2:2887640/3‑100 (MQ=11)
                       ccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGgtcggatccgggggatcg                                                                               >  2:4198239/1‑93 (MQ=11)
                       ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                               <  1:100380/101‑1 (MQ=14)
                       ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTCGGATCCgcggggtccgtcgggtcc                                                                               >  2:5237618/1‑101 (MQ=14)
                        cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGgtcgggtccgtgggagatc                                                                              >  2:1865217/1‑96 (MQ=11)
                        cgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGACCGGTCGGGCCCGTGGGGTCGggggggaccgtggggtcgg                                                                              >  1:4334162/1‑99 (MQ=11)
                        cgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcgg                                                                              >  1:1610355/1‑99 (MQ=255)
                         gtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTcgggtccgtcggatcggt                                                                             >  1:1229353/1‑101 (MQ=11)
                           gggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGT.gg                                                                          <  1:3227381/101‑1 (MQ=11)
                           gggatcgggcggATCCGTAGGATCGGTCGGGTCCGTACGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGGCCGTCGGATCCGT.gg                                                                          <  1:4884471/101‑1 (MQ=11)
                           cggttcggttggGTTCGTGGGATCGGTCTGGTCCGTCGGCTCGGGGGGATCCGTGGGATCTGTCGGGTCAGTGGGATCGgtcgggtccgtcggatcgggcg                                                                           <  1:3304972/97‑6 (MQ=2)
                            ggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGGGGGGACCGT.cgg                                                                         >  2:35237/1‑101 (MQ=11)
                            ggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGT.ggg                                                                         <  2:5336190/101‑1 (MQ=11)
                             gatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGGCGGATCCGT.gggg                                                                        >  1:4479992/1‑100 (MQ=11)
                              gtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGT.CGGac                                                                       >  1:5756475/5‑100 (MQ=11)
                               tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGAtc                                                                      >  1:4562948/1‑101 (MQ=14)
                               tcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGT.GGGGtc                                                                      <  1:4245902/101‑1 (MQ=11)
                               tccgtgggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGT.GGGAtc                                                                      <  1:5516266/95‑1 (MQ=11)
                               tccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGT.GGGAtc                                                                      <  1:972276/98‑1 (MQ=11)
                                ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGAtcg                                                                     >  1:1928200/6‑101 (MQ=14)
                                ccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGT.GGGAtcc                                                                     <  1:2106369/99‑2 (MQ=11)
                                 ctccggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGT.GGGAtccg                                                                    <  2:2378530/98‑3 (MQ=11)
                                 cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGAtcgg                                                                    <  2:2482203/97‑1 (MQ=14)
                                  gtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGT.GGGGtccgg                                                                   >  2:135196/1‑98 (MQ=11)
                                   tggggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCCGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGGTCCGT.GGGAtcggtc                                                                  <  2:337651/96‑1 (MQ=11)
                                    cggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCCGTGGGGTCCGTCGGGTCCGT.GGGAtcgggcg                                                                 >  2:2139626/1‑101 (MQ=14)
                                     gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGAtcggtcgg                                                                <  2:4355420/98‑1 (MQ=11)
                                     gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGAtcggtcgg                                                                >  2:4063857/4‑101 (MQ=11)
                                     gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGAtcggtcgg                                                                >  1:2635436/4‑101 (MQ=11)
                                     gggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGAtcggtcgg                                                                >  1:5404887/4‑101 (MQ=14)
                                     ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGT.CGGAtccgtggg                                                                <  1:1302496/101‑1 (MQ=11)
                                      ggtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTCGGGTCCGTCGGGTCCGT.GGGGtcggtcggg                                                               >  1:3972460/3‑100 (MQ=14)
                                       gtCCGTGCGTTCGGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGAtcggtcgggt                                                              <  2:5157472/100‑3 (MQ=11)
                                         ccGTGGGATCGGTCGGATCCGTGGGATCGGTAGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGAGAtcggaagagcg                                                             >  2:4293822/1‑94 (MQ=11)
                                            taggATCTGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGtc                                                         <  2:1816195/99‑2 (MQ=14)
                                             cgtgtcttgagggTCGGTCGGAGCGGTGGGAGCGGTCGGGTCGGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGT.GGGATCGGTCGGGTCCGTgg                                                        <  2:5663917/91‑1 (MQ=11)
                                             cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTgg                                                        >  2:3906945/2‑101 (MQ=255)
                                             aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGGTCGGTCGGGTCCGtcg                                                        >  1:3876583/2‑99 (MQ=14)
                                              ggATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGATCGGTCGGGTCCGTggg                                                       >  1:2906832/1‑101 (MQ=14)
                                               gATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGa                                                      <  2:2239603/101‑1 (MQ=255)
                                                gtcggtcggGTCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGATCGGTCGGGTCCGTGGGAt                                                     <  1:4668293/100‑1 (MQ=17)
                                                gtcggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGAt                                                     <  1:149410/100‑1 (MQ=14)
                                                aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGAt                                                     >  2:560327/1‑101 (MQ=18)
                                                 tcggtgggGTGGGGGGGGTCGGTAGGGTCGGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGT.CGGGTCCGTGGGGCTGGTGGGAtc                                                    <  1:180903/101‑1 (MQ=14)
                                                   ggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGAtcgg                                                  >  2:1131421/1‑101 (MQ=17)
                                                     tggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGAtcggtc                                                <  1:309885/99‑1 (MQ=255)
                                                          tCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                           >  2:3240537/1‑101 (MQ=255)
                                                            cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGCTCGGTCGGGTCCGTCGGATCGGGCGGGTCCg                                         >  2:5515974/1‑101 (MQ=17)
                                                            cGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGGTCGGTCGGGTCCg                                         >  2:3020036/1‑101 (MQ=17)
                                                               cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGGTCGGTCGGGTCCGTgg                                      >  2:577092/2‑101 (MQ=17)
                                                               cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgg                                      <  1:21695/100‑1 (MQ=17)
                                                                ggATCGGTCGGATCCGTGGGATCGGTCGGGTCCTGCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggg                                     <  1:4979063/101‑1 (MQ=17)
                                                                 gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                                    <  2:3992019/101‑2 (MQ=17)
                                                                  gtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggat                                   <  2:76072/97‑3 (MQ=17)
                                                                     ggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGGTCGGTCGGGTCCGTTGGATCGCGCGGGTCCGTggcatcgg                                >  1:3420117/1‑95 (MQ=14)
                                                                     ggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcgg                                >  1:5664212/1‑96 (MQ=17)
                                                                     ggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCAGTCGGGTCGGTCGGGTCCGTgggatcgg                                >  2:4481748/1‑96 (MQ=17)
                                                                       tcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg                              <  1:285322/101‑1 (MQ=255)
                                                                       tcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtc                              <  2:1046973/101‑8 (MQ=11)
                                                                        cggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGT.GGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcg                             >  1:4142886/1‑93 (MQ=11)
                                                                         gggtCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTggg                            >  2:697725/4‑101 (MQ=255)
                                                                          ggaccGTGGGATGGGGGGGGGCGGTAGGATCTGTCGGGTCCGCCGGGTCCGT.GGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGa                           <  1:3539023/98‑1 (MQ=18)
                                                                          gATCCGTGGGATCGGTGGTGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGT.GGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtcggg                           <  1:2285848/101‑11 (MQ=11)
                                                                            tCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                         <  2:1480128/101‑1 (MQ=255)
                                                                            tCCGTCGGAACGGTCGGGTCGGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc                         <  1:536096/101‑1 (MQ=255)
                                                                                tGGGATCGGTCGGGTCGGTGGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggtc                     <  1:198782/101‑1 (MQ=255)
                                                                                      cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCg               >  1:92421/1‑101 (MQ=255)
                                                                                       ggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGc              >  1:2742879/1‑101 (MQ=255)
                                                                                               ccGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccg      >  1:5083230/1‑101 (MQ=255)
                                                                                                   |                                                                                                     
CGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGT.GGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGC  >  NC_014034/221964‑222163

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