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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0360C→G 3.9% 115.1 ‑0.8 ‑0.01788P207P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (29/43):  ref base (772/941):  total (804/984)
Fisher's exact test strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.77e-01

Only 100 of 2473 total aligned reads displayed.

TCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGG  >  NC_014034/221936‑222136
                                                                                                    |                                                                                                    
  ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                                    >  2:1314129/1‑101 (MQ=255)
       ggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                                               <  1:5093056‑M1/88‑1 (MQ=255)
        gtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGgatcggt                                                                                              <  1:5295200‑M1/89‑4 (MQ=255)
         tgggagcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc                                                                                             <  2:3228432/95‑2 (MQ=11)
         tGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                                             >  1:1100130/1‑101 (MQ=255)
          tggatcggttggGTCCGTGGGATCGGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGtcg                                                                                            <  2:1535841/97‑3 (MQ=11)
          cgggtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGCCGGGGCCGTgg                                                                                            >  1:3256760/8‑101 (MQ=11)
           gggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                                                           <  1:732308‑M1/92‑1 (MQ=255)
            gatcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCggcggg                                                                                          >  2:699247/3‑96 (MQ=11)
              tcggtcggGTCGGTGGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                                        <  2:2812526/101‑1 (MQ=11)
               gagttgactcggtcgggtgcgtgggatcggtcgagtcagtgggaTCGGTCGGGTACGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                                                                       <  2:396892/93‑1 (MQ=11)
                ggtaggATCGGTCGGATCCTTTGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                                                                      <  1:5705475/101‑1 (MQ=11)
                cgtgGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                      >  1:5600343‑M1/5‑101 (MQ=255)
                cgggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAtccgtggggtccg                                                                                      >  1:2643013‑M1/5‑96 (MQ=255)
                    ggatcggtcgggtccgtggggtccgtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCGGCGGGAccgggcgg                                                                                  >  2:4537644/1‑101 (MQ=11)
                    agatcagtcgggtccgtgggatcggtcgggtccgtggGGTCGGTCGGGTCTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                                  <  2:5587046/100‑1 (MQ=11)
                       gcggtgggggccgtgggagcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATc                                                                               <  1:2529775/100‑1 (MQ=11)
                        ggggcggggccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                              <  1:4298196/97‑1 (MQ=11)
                        cggtcgggtcggggggagcggtcgggccggtgggatCGGTCGGCCCGGTGGGAGCGGTCGGTCCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                              <  2:4979660/101‑1 (MQ=2)
                        cggtcgggtccgtgggatcggtcgggtccgtcggatCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                                                                              >  1:465800/1‑101 (MQ=11)
                         ggtcggcgccgtggggtcggtcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                             <  2:3450887/101‑1 (MQ=11)
                         ggtcggatccgtgggatcggtcgggtccgtgggatcGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCCGTCGGGGCCGTGGGATCTGTCGGGTACg                                                                             >  2:594113/1‑101 (MQ=14)
                          gtgggttccgtgggttcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                            <  2:1241918/95‑1 (MQ=11)
                          ggcgggtccgggggagcggtcgggtccgtgggatcggtGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                            <  1:827932/99‑1 (MQ=11)
                           tgggatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGGCCGtg                                                                           >  1:1464306/6‑100 (MQ=11)
                           tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtg                                                                           >  2:4753847‑M2/1‑94 (MQ=255)
                            cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtgg                                                                          >  1:3069372/1‑99 (MQ=11)
                             gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGACGGATCCGTGGGATCGGTCGGATCCGTGGGATGGGTCGGGTCCGTGgg                                                                         >  2:3872353/1‑101 (MQ=14)
                             ggatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTGGGATCCGTGGGAGCGGTGGGATcgggggg                                                                         >  1:4107291/4‑95 (MQ=11)
                                tccgtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                      <  2:4877059/101‑1 (MQ=14)
                                tccgtgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtaggatc                                                                      >  2:1636533‑M2/1‑89 (MQ=255)
                                 ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcg                                                                     <  1:401196‑M2/101‑14 (MQ=255)
                                   gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGTGTCCGTGGGATCGGTCGGGCCCggggggcccgg                                                                   >  1:4987656/1‑91 (MQ=11)
                                       gatcggtcgggtccgtgggatcggtcggATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATGGGTCGGGTCGGTGGGATCGgtcggg                                                               >  1:3493261/1‑101 (MQ=11)
                                         tcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGgtaggggc                                                             >  2:3125699/1‑99 (MQ=14)
                                          cggtcgggtccgtgggatcggtcggATCCGTG.........GGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgggtcc                                                   >  1:922593/1‑101 (MQ=11)
                                          cggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtcc                                                            <  2:5481152‑M2/101‑23 (MQ=255)
                                            gtcggctccggaggatcggtcggGTCCGTAGGATCGGTCGGGTCCCTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgt                                                          <  1:2718118/101‑1 (MQ=11)
                                             tggggtcggtggggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGTTCGgtcgggtccgtc                                                         <  1:4790845/99‑1 (MQ=11)
                                              gggatcggtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcg                                                        <  1:4892356/94‑1 (MQ=14)
                                              cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtgg                                                        <  1:3730649/101‑3 (MQ=11)
                                              cgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGGGCCGTGGGATCGgtcgggtccgtgg                                                        >  2:3747103/1‑99 (MQ=11)
                                               ggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                       >  2:1408207/4‑101 (MQ=255)
                                                ggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtggga                                                      <  1:5195530‑M2/101‑29 (MQ=255)
                                                  tccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtc                                                    <  1:1760653/101‑1 (MQ=11)
                                                  tccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcggatc                                                    <  2:1010228/101‑1 (MQ=11)
                                                   ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                   <  1:100380/101‑1 (MQ=14)
                                                   ccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcggatcg                                                   <  2:539255/101‑2 (MQ=255)
                                                    gggggggtccgtgggGTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGtcgggtccgtcgggtccg                                                  <  1:1666029/98‑1 (MQ=11)
                                                       gggatcggtcggGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcg                                               >  2:1939776‑M2/1‑66 (MQ=255)
                                                       gggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGatccgtgggatcggtcgggtccgtgggatcggtcg                                               <  2:2326511‑M2/101‑36 (MQ=255)
                                                       aggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGGCCGTCGGGTCCGTgg                                               >  2:4635132/2‑101 (MQ=14)
                                                        gggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTggg                                              <  2:5180967/92‑1 (MQ=14)
                                                         gatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGa                                             <  1:5485271/101‑1 (MQ=11)
                                                           tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtc                                           >  1:4562948/1‑101 (MQ=14)
                                                            cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                          <  1:3607972‑M2/101‑41 (MQ=255)
                                                             cgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcgg                                         >  2:3596699/5‑101 (MQ=14)
                                                                cggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcg                                      <  1:3837058/101‑1 (MQ=14)
                                                                 gggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                     >  2:4063857/4‑101 (MQ=11)
                                                                 ggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggg                                     <  1:4210960‑M2/101‑46 (MQ=255)
                                                                 ggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                     >  2:112852/1‑101 (MQ=17)
                                                                 ggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgcggg                                     <  1:1862950‑M2/101‑46 (MQ=255)
                                                                 ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg                                     <  1:1302496/101‑1 (MQ=11)
                                                                   gtCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgggg                                   >  1:2378664/2‑99 (MQ=14)
                                                                     ccGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                 <  2:806086‑M2/101‑50 (MQ=255)
                                                                      cGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCg                                >  1:1572276/1‑101 (MQ=14)
                                                                       gTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                               <  2:5064903‑M2/101‑52 (MQ=255)
                                                                       ctgGGATCGGTCGTGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggt                               >  2:4407851‑M2/2‑50 (MQ=255)
                                                                         gggATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcggtcg                             <  2:5291407‑M2/101‑54 (MQ=255)
                                                                         gggATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                             <  1:3157222‑M2/101‑54 (MQ=255)
                                                                            aTCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAGCGGTCGGATCCGTGGGAt                          >  2:212497/1‑101 (MQ=18)
                                                                              cggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                        <  2:4158521/101‑1 (MQ=14)
                                                                              ccgtcggatCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatcc                        <  1:393593‑M2/99‑59 (MQ=255)
                                                                               ggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                       <  2:3315289‑M2/101‑60 (MQ=255)
                                                                               ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGGTCCGTGGGATCGGCCGGGGCCGGGGGGtcgg                       >  2:2018775/1‑101 (MQ=14)
                                                                               cgtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGCCCGTGGGATCGGTCGGATCCGTGGGAtcgg                       >  2:2476782/5‑101 (MQ=255)
                                                                                 tcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc                     <  1:2033326/101‑1 (MQ=14)
                                                                                  cggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgg                    <  1:4354631‑M2/101‑63 (MQ=255)
                                                                                  cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                    >  1:4470111/1‑101 (MQ=17)
                                                                                  cggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcg                    <  2:425986/101‑1 (MQ=17)
                                                                                   ggatCCGTAGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                   >  2:3065465/4‑101 (MQ=14)
                                                                                     gTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggt                 >  2:182818‑M2/1‑36 (MQ=255)
                                                                                     gTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                 >  2:885214‑M2/1‑36 (MQ=255)
                                                                                      gccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                >  2:4325240‑M2/2‑35 (MQ=255)
                                                                                       ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcg               >  1:1221677‑M2/1‑34 (MQ=255)
                                                                                       ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcg               >  2:5684889‑M2/1‑34 (MQ=255)
                                                                                       ccGGGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc               <  1:1965922/101‑1 (MQ=17)
                                                                                           gggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgggggatcggtcg           <  2:2954158‑M2/101‑72 (MQ=255)
                                                                                           gggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtggggtcggtgg           <  1:2246451‑M2/101‑72 (MQ=255)
                                                                                             ggtcggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggg         <  1:2939982/99‑1 (MQ=255)
                                                                                             gatcggtcggatccgtgggatcggtcggatccgtaggatcggtcgggtccgtaggatcggtcgggtccgtcggatcggtcggatccgtgggatcggtcggg         <  2:2708879‑M2/101‑74 (MQ=255)
                                                                                             gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGTTCGGCCGGGTCCGTggga         >  1:2138483/1‑100 (MQ=17)
                                                                                              aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGGCGGGTCCGTGGGGTCGGGCGGGTCCggcggct        >  1:5389854/1‑95 (MQ=14)
                                                                                                ccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatcc      >  2:2606349‑M2/3‑25 (MQ=255)
                                                                                                ccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtcc      >  2:2225539‑M2/3‑25 (MQ=255)
                                                                                                 ggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccg     >  1:3370285‑M2/1‑24 (MQ=255)
                                                                                                  gtcggatccgtgggatcggtcggttccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgt    >  2:1716943‑M2/1‑23 (MQ=255)
                                                                                                  gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTaggatcggt    >  1:1362469/1‑95 (MQ=14)
                                                                                                   tcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg   <  1:285322/101‑1 (MQ=255)
                                                                                                   tcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc   >  1:603520/1‑94 (MQ=14)
                                                                                                    |                                                                                                    
TCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGG  >  NC_014034/221936‑222136

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