Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,976 A→G 3.0% D227D (GAT→GAC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9760A→G 3.0% 31.8 ‑0.2 ‑1.2712D227D (GAT→GACRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/11):  ref base (329/361):  total (340/372)
Fisher's exact test strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.66e-02

Only 100 of 1668 total aligned reads displayed.

TCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGT  >  NC_014034/221876‑222076
                                                                                                    |                                                                                                    
  tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatc                                                                                                    >  2:4551538‑M1/91‑101 (MQ=255)
  tccgtgggatcggtcgggtccgtgggatcggacggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtctgatccatgggata                                                                                                    <  1:82816‑M1/11‑2 (MQ=255)
   ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                                   >  1:2654881‑M1/90‑101 (MQ=255)
     gtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggt                                                                                                 <  2:3714980‑M1/14‑1 (MQ=255)
      aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                                                <  1:3113131/101‑1 (MQ=255)
         gatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggttgga                                                                                             <  2:475288‑M1/18‑5 (MQ=255)
          gtcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGt                                                                                            >  1:4668117‑M1/83‑101 (MQ=255)
             ggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGGTCGGTCGGGTCCg                                                                                         <  2:3861458‑M1/22‑1 (MQ=255)
              gtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGGTCGGTCGGGTCCGt                                                                                        <  1:1829014‑M1/23‑1 (MQ=255)
                cggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTgg                                                                                      <  2:1969702‑M1/25‑1 (MQ=255)
                cggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggaTCCGTAGGATCGGTCGGGTCcgtag                                                                                      <  1:2744327‑M1/25‑3 (MQ=255)
                   atccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGAt                                                                                   >  1:4372920‑M1/74‑101 (MQ=255)
                      cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                >  2:3183580‑M1/71‑101 (MQ=255)
                      cgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                                <  1:933368‑M1/31‑1 (MQ=255)
                        tCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                              <  2:653487/101‑1 (MQ=255)
                               ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcgggtccgtgggatcggtcgggaccg                                                                       <  2:1745798‑M1/46‑6 (MQ=255)
                                gtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCcgg                                                                      <  2:2013186‑M1/47‑2 (MQ=255)
                                  cggatccgtgggattggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                                    >  2:1382626‑M1/53‑101 (MQ=255)
                                  cggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                                    >  2:4072664‑M1/53‑101 (MQ=255)
                                   gggcgcgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                                   >  2:1050979‑M1/58‑101 (MQ=255)
                                    gatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtgggc                                                                  <  1:1002062‑M1/51‑2 (MQ=255)
                                            gggtccgtcgggtccgtgggatcggtcggatccgggggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                                          >  1:2228136‑M1/49‑101 (MQ=255)
                                              atcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCCTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTcgggc                                                        <  1:99564‑M1/61‑2 (MQ=255)
                                               tcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                       >  1:3201681‑M1/40‑101 (MQ=255)
                                                    cgggtccgtggggtccgtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                  <  1:3578048‑M1/61‑1 (MQ=255)
                                                    cgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                  <  1:1063949‑M1/61‑1 (MQ=255)
                                                    cgggtccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCcgtag                                                  >  1:584058‑M1/41‑99 (MQ=255)
                                                     gggtccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                                 >  2:4167511‑M1/40‑101 (MQ=255)
                                                         ccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                             >  2:3436902‑M1/36‑101 (MQ=255)
                                                           gtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                           <  1:7498‑M1/68‑1 (MQ=255)
                                                             cGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                         >  1:3308926/1‑101 (MQ=255)
                                                                atcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATGGGTCGGAt                                      <  2:568295‑M1/79‑1 (MQ=255)
                                                                 tcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGACcgtgggatcggctgggcg                                     <  1:2861935‑M1/74‑8 (MQ=255)
                                                                 tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGCGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                     <  1:4017661/101‑1 (MQ=255)
                                                                   ggTGGGGTCGGGGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                                   <  1:2322773/100‑1 (MQ=11)
                                                                    gTGGGGTCCGGCGGAGCGGGCGGTTCGGGGGGACCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                  <  1:173033/101‑1 (MQ=11)
                                                                    gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                                  >  1:653915/1‑101 (MQ=255)
                                                                     tcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTg                                 <  2:2091289‑M1/78‑1 (MQ=255)
                                                                      gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                >  1:2581653/1‑99 (MQ=1)
                                                                       gggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               >  1:1987064‑M1/16‑101 (MQ=255)
                                                                       gggtccgtcggAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                               <  2:2950781/95‑1 (MQ=11)
                                                                       ggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                               <  1:150902‑M1/80‑1 (MQ=255)
                                                                         tgccgttgggtcgggcggctccgtcggatcggtcgggtccgtcggtTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                             <  2:2390715/91‑1 (MQ=11)
                                                                         gttcgtgggatcggtcggctccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGGGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                             <  2:363319/94‑1 (MQ=11)
                                                                          tccttgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtc                            <  2:167100/95‑1 (MQ=11)
                                                                          tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGAGCCGTGGGATCGGTCGGATCCGTGGGCtc                            >  2:1273390/7‑101 (MQ=11)
                                                                             gtgggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTGGGAtcggg                         >  1:2516983/4‑100 (MQ=11)
                                                                             gtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                         >  1:3848245‑M1/16‑101 (MQ=255)
                                                                              tggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtc                        <  2:2869155/93‑1 (MQ=14)
                                                                              tgggatcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgggtccgtc                        >  2:1449741‑M2/3‑92 (MQ=255)
                                                                              tgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc                        <  2:4378221/99‑1 (MQ=11)
                                                                              ggggatcggtcgggttcgtgggagcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                        <  1:4595084/99‑1 (MQ=14)
                                                                              ggggagctggcgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtc                        <  1:3399978/91‑1 (MQ=11)
                                                                               ttgatgggtcgggttcgtggggtttgtcggatccgtgcGATCTGTCGGATCTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                       <  2:4659081/99‑1 (MQ=11)
                                                                               ggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtcg                       <  1:2963273/94‑1 (MQ=14)
                                                                               gggggccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtcg                       <  2:247712/94‑1 (MQ=14)
                                                                                ggatcggtcggaTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                      <  1:1409157‑M1/89‑1 (MQ=255)
                                                                                 gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGa                     <  1:4508738‑M2/101‑2 (MQ=255)
                                                                                  gttggtcgggtcggggggatggggcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt                    <  1:1354481/98‑1 (MQ=11)
                                                                                  atcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGATCCTTGGGATCGGTcgggc                    <  1:1203595‑M1/91‑2 (MQ=255)
                                                                                   tcggtgggatccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                   >  1:3028671/1‑101 (MQ=11)
                                                                                   tcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatc                   <  1:2172344‑M2/101‑4 (MQ=255)
                                                                                   tccgtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGAATCGGTCGGATCCGTGGGATCGGTCGGGTc                   <  1:4309850‑M1/92‑1 (MQ=255)
                                                                                    tggtcggatctgtgggatcggtcggacccgaaggctcGGTCGGGTCCGTAAGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                  <  2:734405/100‑1 (MQ=11)
                                                                                    gggtcgggtccgtggggnccgtcgggtccgtcggatcGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                  <  2:94140/100‑1 (MQ=11)
                                                                                    cggtgggacccgggggatgggtcgggtccgtgggatcggtcggGTCCGGCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  1:3946449/92‑1 (MQ=11)
                                                                                    cggtcgtgtccgtgggctccgtcggatccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  2:4468024/101‑1 (MQ=11)
                                                                                    cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatcc                  <  2:1891770‑M2/101‑5 (MQ=255)
                                                                                    cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGatcg                  <  2:2555714‑M2/101‑5 (MQ=255)
                                                                                    cggtcggcgccgtgggatcggtcgggtccgtcggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  2:849055/101‑1 (MQ=11)
                                                                                    cggtcggatccgtgggatcggtcggatccgtaggatCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                  >  2:3059258/1‑101 (MQ=11)
                                                                                    cggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                  >  2:3980244‑M1/9‑101 (MQ=255)
                                                                                    cggtccggtccgtgggatcgttcgggtccgtcggatCGGTCGCGTCCGTTGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  2:3327970/101‑1 (MQ=11)
                                                                                     ggtcgggtccgtgggatctgtcgggtccgtggtatcGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcgg                 <  1:1923270‑M2/101‑17 (MQ=255)
                                                                                     ggtcgggtccgtgggatcggtcgggtccgtggggtcCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGATCGGATCCGTCGGATCGGTCGGGTCCg                 >  1:204614/1‑101 (MQ=11)
                                                                                     ggtcgggtccgtgggatcggtcgggtccgtgggatcGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccg                 <  2:2228012‑M2/101‑6 (MQ=255)
                                                                                     ggtcggatccgttggatcggtcgtgtccgtgggatcGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcgg                 >  2:3548798‑M2/1‑85 (MQ=255)
                                                                                      gtgggatccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                >  1:4223082/7‑101 (MQ=11)
                                                                                      gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGTCGCTGGGTTCCGTCGGCTCGGCGGGCTCGCTCGGGTCCgg                >  2:1121807/1‑100 (MQ=11)
                                                                                      gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgt                <  1:1132289‑M2/101‑7 (MQ=255)
                                                                                       gcgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGGCCGGCGGATCCGTGGGATCGGTCGGATCCGtg               <  1:417595/100‑2 (MQ=14)
                                                                                        cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTccagat              >  2:422644/1‑97 (MQ=11)
                                                                                        cgggtccgtgggatcggtcgggtccgtgggatcggtCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatacgtgg              <  2:1202696‑M2/101‑9 (MQ=255)
                                                                                         gggttcgtgggtcctgtcgggtgcgttggatcggtcGGGTCCTTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgg             <  2:1557644/101‑1 (MQ=1)
                                                                                         gggtccgtggggtccgtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGGGGGATAGGTCGGGTCCGTGGGAtcggtcgcgagcgtgga             >  1:504362/1‑91 (MQ=11)
                                                                                         gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTGCGTGGGCTCCGTCGGTTCCGTGgg             >  2:3241498/1‑101 (MQ=11)
                                                                                         gggtccgtgggatcggtcgggtccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtggg             <  1:2412931‑M2/101‑10 (MQ=255)
                                                                                         gggtccgtgggatcggtcgggtccgtcggatcggtgGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGATCGGTCGGGCCCGTGgg             >  1:4439363/1‑101 (MQ=14)
                                                                                         ggatccgtgggatcggtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtggg             <  1:2105931‑M2/101‑10 (MQ=255)
                                                                                          gatccgtggggtccgtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCGGCGGGATCGCTAGGGAcccagcga            >  2:2654727/3‑95 (MQ=2)
                                                                                           gtccgtgggatcggtcgggtccgggggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATccttgtgat           >  2:2318584/1‑94 (MQ=14)
                                                                                           ggccgtgggatcgggcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAt           <  1:3745268/99‑1 (MQ=14)
                                                                                            gccgtgtgctcggtccggtccgtggggtccgtcggACCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc          <  2:431975/100‑1 (MQ=11)
                                                                                             ccgtggggtccgtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGGGGGATCCGTGGGAtcg         >  1:3479400/1‑101 (MQ=11)
                                                                                             ccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatcg         <  1:1630305‑M2/101‑14 (MQ=255)
                                                                                             ccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcg         <  1:3643766‑M2/101‑14 (MQ=255)
                                                                                               gtggggtccgtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGGGCGGGAGGGTCGGTGggggcggg       >  2:2138345/1‑95 (MQ=0)
                                                                                                tggggtccgtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtcggcc      >  1:412188/1‑99 (MQ=11)
                                                                                                tgggatcggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc      >  2:2752726/1‑101 (MQ=11)
                                                                                                  ggatcggtcgggtccgtgggatcgggcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg    <  1:1295110/101‑1 (MQ=17)
                                                                                                    |                                                                                                    
TCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGT  >  NC_014034/221876‑222076

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