Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,051 G→C 1.8% 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.8% 35.7 ‑3.2 ‑0.61481T202T (ACC→ACGRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/22):  ref base (660/788):  total (670/811)
Fisher's exact test strand distribution p-value = 4.23e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.54e-01

Only 100 of 2386 total aligned reads displayed.

CGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGC  >  NC_014034/221951‑222151
                                                                                                    |                                                                                                    
cggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcg                                                                                                      <  2:3275539/101‑1 (MQ=11)
ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                                                                      >  1:3840404‑M1/6‑101 (MQ=255)
ccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTGCGTGGGGTCGGCCGGGGCAGTGGGAtcg                                                                                                      >  2:317867/3‑101 (MQ=14)
 cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGAtcgg                                                                                                     <  2:1674742‑M1/97‑1 (MQ=255)
 cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                                                                                     <  1:796675‑M1/97‑1 (MQ=255)
  gtggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggt                                                                                                    >  2:5393539/10‑101 (MQ=11)
  gtgGGATCGGTCGGTTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                                    >  2:4890336‑M1/4‑101 (MQ=255)
  gtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                                    >  2:3098315‑M1/4‑101 (MQ=255)
  gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                                                                    >  2:5478054/1‑101 (MQ=255)
   tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                                                                   >  1:5652124‑M1/3‑101 (MQ=255)
   tgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                                                   <  1:1487915‑M1/99‑1 (MQ=255)
   ggGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                                                                   >  1:4108048‑M1/3‑101 (MQ=255)
    cggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtgg                                                                                                  >  2:1770851/1‑99 (MQ=255)
     ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGGGGGGTCGGTCGGGTCCGTGGGGtcggtcgg                                                                                                 >  2:2047502/1‑101 (MQ=11)
     ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGGGGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGCGGGTCGGTCGGAtcggttcg                                                                                                 >  2:415871/1‑98 (MQ=11)
       atcggtcgggtccgtggcatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcggat                                                                                               <  2:5195244/101‑3 (MQ=14)
        tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtagatcg                                                                                              >  2:2655326/1‑95 (MQ=14)
        ggagtgcggtccgtgggagcggtcggggccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                                              <  1:2647263/94‑1 (MQ=14)
         ctggcgtatccgtgggatcggtcgggtccgttggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                                                             <  1:1305871/93‑1 (MQ=11)
         cggtgggctccgcgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                             <  2:2147795/101‑1 (MQ=11)
         cggtcgggtccgtgggatcggtcgggtccgtcggatCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                                                                                             >  1:465800/1‑101 (MQ=11)
         cgggcgggtccgtgggatcgggcgggtccgtggggtCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                             <  1:1674841/101‑1 (MQ=11)
          ggtgggatccgtgggatcggtcgggtccgtgggatcGGTCGGGCCCGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGCCCGCTGGATCGGTAGGGTCCg                                                                                            >  1:2941113/1‑101 (MQ=11)
          ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGGTCGGACGGGTCCGTGGGGTCCGTCGGCTCCGTGGGGTCGGTCGGGTCCg                                                                                            >  1:5339952/1‑101 (MQ=11)
          ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcggagccg                                                                                            >  2:5094191/1‑96 (MQ=11)
          ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccg                                                                                            >  2:4069832‑M2/1‑96 (MQ=255)
          ggtcggatccgtgggatcggtcgggtccgtgggatcGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCCGTCGGGGCCGTGGGATCTGTCGGGTACg                                                                                            >  2:594113/1‑101 (MQ=14)
          ggtcggatccgtgggatcggtcggatccgtgggatcGGTCGGGTCCGCGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg                                                                                            >  2:251845/1‑101 (MQ=11)
           gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgt                                                                                           >  1:3262635‑M2/1‑95 (MQ=255)
            tcgggtccgtggggtcggtcgggtccgtcggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGtg                                                                                          <  1:3739548/101‑2 (MQ=14)
            tcgggtccgtgggatcggtcggatccgtgggatcggTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtg                                                                                          <  2:3788433‑M2/101‑8 (MQ=255)
             cgggtccgtgggatcggtcgggtccgtggggtccgtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtgg                                                                                         >  1:5757283/1‑99 (MQ=14)
             cggatccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTGGGGTCCGTGGGGTCGGTGGGGTCCGtgg                                                                                         >  2:3515628/1‑99 (MQ=11)
              gggtccgtggggtcggtcgggtccgtcggatcggtcGGGTCCGTGGGATCGCTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgg                                                                                        <  1:4894875/101‑1 (MQ=11)
              gggtccgtgggatcggtcgggtccgtcggatcggtgGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGCCGGATCGGTCGGGTCCGTGgg                                                                                        >  2:3027202/1‑101 (MQ=11)
               ggtccgtgggatcggtcgggtccgtgggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGCGATCGGTCGTGCCCGTGCGATCCGTGggggctgaggga                                                                                       >  1:2272670/1‑92 (MQ=11)
               ggtccgtgggatcggtcgggtccgtcggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGCCCGTGGGa                                                                                       >  1:2147917/1‑101 (MQ=14)
               gatccgtgggatcggtcgggtccgttggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                                                       >  1:3943457/3‑101 (MQ=11)
                  ccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcg                                                                                    >  1:1105026‑M2/1‑88 (MQ=255)
                   cgtggggtcggtcgggtccgtggggtccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                   <  1:5005178/101‑1 (MQ=11)
                   cgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcgg                                                                                   >  1:3146404‑M2/1‑87 (MQ=255)
                    gtgggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                                  <  1:1223071/101‑1 (MQ=11)
                     ggggatcgggcggatccgtgggatcgggcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGAtcggtc                                                                                 <  1:196957/100‑1 (MQ=11)
                       ggatcgggcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                               <  1:3080770/101‑1 (MQ=11)
                           cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtcc                                                                           >  1:55105‑M2/1‑79 (MQ=255)
                           cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGCCGgtcggggcc                                                                           >  1:970844/1‑101 (MQ=11)
                            ggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGATCGggcgggtccg                                                                          >  1:5109042/1‑101 (MQ=11)
                             gtcggatccgtgggatcngtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgt                                                                         <  2:2395371/101‑1 (MQ=255)
                             gtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgg                                                                         <  2:3668044‑M2/101‑25 (MQ=255)
                                ttatccgtgggatcgggcggTTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                                                                      <  2:1059168/98‑1 (MQ=14)
                                gggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTGGGatccgtgggatcggtcggatccgtggg                                                                      <  1:3388560‑M2/101‑28 (MQ=255)
                                 ggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGGCCGTGGGATCGgtcgggtccgtggga                                                                     >  2:4548153/1‑100 (MQ=11)
                                 gatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTCGGGCGGgcggggtccgtcgga                                                                     >  2:2887640/3‑100 (MQ=11)
                                   tccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtc                                                                   <  1:1760653/101‑1 (MQ=11)
                                    ccgtgggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                  <  1:100380/101‑1 (MQ=14)
                                    ccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                  <  2:1165415/101‑1 (MQ=11)
                                     cgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTCCGTGGGATCGgtcgggtcggtggggcccc                                                                 >  2:1866017/1‑100 (MQ=11)
                                     cgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCCGtggggtccgtcgggtccg                                                                 >  1:2042297/1‑101 (MQ=11)
                                       tgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtc                                                               >  2:4330856‑M2/1‑67 (MQ=255)
                                       tgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTg                                                               >  1:4350263/1‑101 (MQ=11)
                                            tcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtc                                                          >  1:4562948/1‑101 (MQ=14)
                                            tcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGtc                                                          <  1:3505051/101‑1 (MQ=11)
                                                  ggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgcggg                                                    <  1:1862950‑M2/101‑46 (MQ=255)
                                                  ggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGGCAGGTCCGTGGGATAGGACGGGTCCGTCGTGTTCGTCGGGtacgtcgg                                                    >  2:2265183/1‑101 (MQ=2)
                                                   ggtccGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                   >  1:4482020‑M2/3‑55 (MQ=255)
                                                    aTCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAt                                                  >  1:4282987/1‑101 (MQ=11)
                                                     tCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                                                 >  2:70385‑M2/1‑53 (MQ=255)
                                                      ccGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                                <  2:806086‑M2/101‑50 (MQ=255)
                                                      ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTcc                                                >  2:1583958/1‑101 (MQ=14)
                                                         tGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                             <  1:4796944‑M2/101‑53 (MQ=255)
                                                         tGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGTGGGGTCCGTGGGAtcggtgggggccggg                                             >  2:5586102/1‑94 (MQ=11)
                                                          cgtgtcttgagggTCGGTCGGAGCGGTGGGAGCGGTCGGGTCGGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                            <  2:5663917/91‑1 (MQ=11)
                                                          cggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                            >  2:3906945/2‑101 (MQ=255)
                                                            gATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                          >  2:4989381/1‑101 (MQ=17)
                                                                ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcggggccgtgggatcggtcgggtacg                                      <  2:1337536‑M2/101‑60 (MQ=255)
                                                                 gtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgg                                     <  2:2754140‑M2/101‑61 (MQ=255)
                                                                  gcgggTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtg                                    >  2:4643506‑M2/2‑40 (MQ=255)
                                                                   cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtag                                   >  2:4419066‑M2/1‑39 (MQ=255)
                                                                    gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGAtcggtcgg                                  >  2:1275622/1‑101 (MQ=17)
                                                                     ggTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                 <  1:885290‑M2/101‑65 (MQ=255)
                                                                      gTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggat                                >  2:885214‑M2/1‑36 (MQ=255)
                                                                       tccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatc                               >  1:4207384‑M2/1‑35 (MQ=255)
                                                                       tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatc                               >  1:4013512‑M2/1‑35 (MQ=255)
                                                                       tCCGTGGGATCGGTCTGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                               <  1:3639250/101‑1 (MQ=17)
                                                                        ccgtgggatcggtcggatccgtgggatcggtcggatccgtaggatcggtcgggtccgtaggatcggtcgggtccgtcggatcggtcggatccgtgggatcg                              <  1:5227592‑M2/101‑68 (MQ=255)
                                                                        ccGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTcc                              <  1:1584052/101‑1 (MQ=11)
                                                                         cGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCAGTCGGGTCCGTGGGATCGGACGGGTACGTAGGAGCGGTCGGGGCCg                             >  2:4763767/1‑101 (MQ=2)
                                                                          gtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggt                            <  1:2533077‑M2/101‑70 (MQ=255)
                                                                           tggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtc                           >  2:3385031‑M2/1‑31 (MQ=255)
                                                                           tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc                           <  1:2395505‑M2/101‑71 (MQ=255)
                                                                           tgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtc                           <  1:672467‑M2/101‑71 (MQ=255)
                                                                           tgggatcggtcggatccgtgggatcggtcggatccgtaggatcggtcgggtccgtaggatcggtcgggtccgtcggatcggtcggatccgtgggatcggtc                           <  1:2742156‑M2/101‑71 (MQ=255)
                                                                              gATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggga                        <  2:3992019/101‑2 (MQ=17)
                                                                                tcgggcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCt                      <  2:5294910/101‑1 (MQ=255)
                                                                                   gtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTaggatcggt                   >  1:1362469/1‑95 (MQ=14)
                                                                                    tcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtg                  <  1:285322/101‑1 (MQ=255)
                                                                                     cggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGACCGGTCGGGTCCGTGGGGCtggtgg                 >  2:4260396/1‑101 (MQ=255)
                                                                                         gccGTGGGGTCGGTCGGGTCCGTAGGGTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtc             <  1:5412195/100‑1 (MQ=255)
                                                                                          ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtagggtccgtgggatcggtcgggtccgtcggatcg            <  2:1972480‑M2/101‑86 (MQ=255)
                                                                                            gTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt          >  2:1474427/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGC  >  NC_014034/221951‑222151

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