Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9580A→G 4.9% 85.0 ‑4.9 ‑3.1942D233D (GAT→GACRCAP_rcc00185hypothetical protein
Rejected: Prediction has significantly lower quality scores supporting polymorphism compared to reference.
Reads supporting (aligned to +/- strand):  new base (34/14):  ref base (457/437):  total (491/451)
Fisher's exact test strand distribution p-value = 7.66e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.82e-04

Only 100 of 1900 total aligned reads displayed.

GCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221859‑222058
                                                                                                   |                                                                                                    
gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggat                                                                                                     >  1:281739‑M1/98‑101 (MQ=255)
     tgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtc                                                                                                >  1:3759561‑M1/93‑101 (MQ=255)
     ggggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc                                                                                                >  2:1075875‑M1/93‑101 (MQ=255)
       cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgg                                                                                              >  2:6200380/1‑101 (MQ=255)
          tcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatc                                                                                           >  1:6714087‑M1/88‑98 (MQ=255)
          tcggtcgggtccgtaggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtc                                                                                           <  1:520491‑M1/14‑1 (MQ=255)
          tccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtc                                                                                           <  2:6875531‑M1/14‑1 (MQ=255)
          gACGCCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                           <  1:2900653/101‑1 (MQ=255)
          gACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                           <  1:4949506/101‑1 (MQ=255)
           cggtcggatccgtgggatcggtcggggccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatcc                                                                                          >  1:4514654‑M1/87‑101 (MQ=255)
               cGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgg                                                                                      <  2:1172255/101‑1 (MQ=255)
                  gtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggct                                                                                   <  1:1526010‑M1/22‑3 (MQ=255)
                  atccgtgggatcggtcgggtccgtgggatcggtcggatccgtggggtcggtcgggtccgtggggtcggtcgggtccgtggggtcggtcgggtccgtgggat                                                                                   <  2:853279‑M1/22‑1 (MQ=255)
                    ccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                 <  1:6066156‑M1/24‑1 (MQ=255)
                    ccGAGGCTTCGGTGGTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcg                                                                                 <  1:4251252/101‑1 (MQ=255)
                      gtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgc                                                                               <  1:6195600‑M1/26‑7 (MQ=255)
                       tgggatcgggcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc                                                                              >  2:4619776‑M1/75‑101 (MQ=255)
                       ggggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc                                                                              >  1:4160266‑M1/75‑101 (MQ=255)
                          gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcgga                                                                           <  2:6062020‑M1/30‑2 (MQ=255)
                          cTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcggg                                                                           >  2:3536175/1‑101 (MQ=255)
                           atcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggg                                                                          <  2:325060‑M1/31‑2 (MQ=255)
                             gggtcggctccctgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                                                        >  2:2009260‑M1/69‑101 (MQ=255)
                              ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggggtcggtcgggtccg                                                                       <  1:4853728‑M1/34‑1 (MQ=255)
                               ggcgggtccgtgggatcggtccggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgt                                                                      >  1:1246028‑M1/67‑101 (MQ=255)
                                tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgggcgggtccgtg                                                                     >  1:2634251/1‑101 (MQ=255)
                                 cggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgggg                                                                    <  2:2589036‑M1/37‑4 (MQ=255)
                                  gggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggg                                                                   >  1:3438935‑M1/64‑101 (MQ=255)
                                  gggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggg                                                                   >  1:5887640‑M1/64‑101 (MQ=255)
                                  ggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggg                                                                   <  1:3913461‑M1/38‑1 (MQ=255)
                                   ggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                  <  1:2053619‑M1/39‑1 (MQ=255)
                                   gatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggga                                                                  >  2:300676‑M1/63‑101 (MQ=255)
                                      tggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGatcggtcgggtccgtggggtcggtcggggccgggggggcg                                                               >  1:5437386/1‑97 (MQ=38)
                                      ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcg                                                               >  2:3636082‑M1/60‑101 (MQ=255)
                                         tgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtg                                                            <  2:2856496‑M1/45‑2 (MQ=255)
                                          gggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtggggtccgtggggtcggtcg                                                           <  2:2133946‑M1/46‑1 (MQ=255)
                                           ggatcggtcgggtccgtgggttcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgg                                                          >  1:4388395‑M1/55‑101 (MQ=255)
                                            gatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatctgtcggg                                                         <  1:31251‑M1/48‑2 (MQ=255)
                                              gcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATc                                                       >  1:3239274‑M1/52‑101 (MQ=255)
                                                   cggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                                  >  2:5930173‑M1/47‑101 (MQ=255)
                                                   cGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                  >  1:3601882/1‑101 (MQ=255)
                                                    ggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGGTCCGAggg                                                 <  1:3318439‑M1/56‑1 (MQ=255)
                                                       tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                                              >  2:4260484‑M1/43‑101 (MQ=255)
                                                        gcgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                             >  1:6404096‑M1/42‑101 (MQ=255)
                                                        ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATccgtggggtcc                                             >  1:2240484‑M1/42‑97 (MQ=255)
                                                        ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATccgtggggtcc                                             <  1:4711948‑M1/60‑5 (MQ=255)
                                                        ccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                             <  1:2423695‑M1/60‑1 (MQ=255)
                                                         cgtgggatcggtcggggccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                            >  1:5728556‑M1/41‑101 (MQ=255)
                                                          gtcggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                           >  1:6586424‑M1/40‑101 (MQ=255)
                                                              gatcggtcggatccgtgggatcggtcggctccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatcggtcgga                                       >  2:4103107‑M1/36‑100 (MQ=255)
                                                              cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGACGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg                                       >  2:3452329/1‑101 (MQ=255)
                                                               atcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGgatcggtcggat                                      <  1:5960859‑M1/67‑3 (MQ=255)
                                                               atcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGggtccgtcggat                                      <  1:3019169‑M1/67‑3 (MQ=255)
                                                               agcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                                      >  1:6293182‑M1/35‑101 (MQ=255)
                                                                tcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                     >  1:1104265‑M1/34‑101 (MQ=255)
                                                                tcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                     >  2:3128864‑M1/34‑101 (MQ=255)
                                                                 cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc                                    >  2:3950424‑M1/33‑101 (MQ=255)
                                                                 cggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                    >  1:7109491‑M1/33‑101 (MQ=255)
                                                                   gtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCcgg                                  <  2:4350184‑M1/71‑2 (MQ=255)
                                                                    tcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTg                                 <  2:1752245‑M1/72‑1 (MQ=255)
                                                                      gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTggg                               >  1:4337120‑M1/28‑101 (MQ=255)
                                                                      gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                               <  1:5328799‑M1/74‑1 (MQ=255)
                                                                      gggtccgtgggatcggtcgggtccgtcGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               >  2:7109892‑M1/28‑101 (MQ=255)
                                                                       ggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                              <  2:1495042‑M1/75‑1 (MQ=255)
                                                                       gatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                              <  1:16912‑M1/75‑1 (MQ=255)
                                                                        gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                             >  2:5284276/1‑101 (MQ=255)
                                                                         tccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                            >  2:4914565‑M1/25‑101 (MQ=255)
                                                                         tccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                            >  1:4611217‑M1/25‑101 (MQ=255)
                                                                         ggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                            >  2:1331066/1‑101 (MQ=255)
                                                                          cggTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcggggcc                           >  1:5086008/3‑97 (MQ=11)
                                                                            gtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCAGTGGGatcggc                         <  2:228479‑M1/80‑2 (MQ=255)
                                                                            gtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcgga                         <  2:273993‑M1/80‑2 (MQ=255)
                                                                            gtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                         <  1:1805145‑M1/80‑1 (MQ=255)
                                                                            ggcGGATGCGGGGGCTCGGTCGGATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                         <  1:1915966/99‑1 (MQ=255)
                                                                            gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggg                         >  1:2721378/1‑100 (MQ=255)
                                                                             tgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGatcggtg                        <  2:1627964‑M1/81‑2 (MQ=255)
                                                                               ggATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGCCGGGTCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGGTCCGTGGGAtcggtcgg                      >  2:4743081/1‑101 (MQ=255)
                                                                                     gtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                <  1:954688‑M1/89‑1 (MQ=255)
                                                                                      tgggagcggtcggATACGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCATCGGGTCCGTGGGATCGGtc               >  1:6860707/7‑100 (MQ=14)
                                                                                      tcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTg               <  2:1867106‑M1/90‑1 (MQ=255)
                                                                                      gggggtccgtgGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg               >  1:6997326‑M1/12‑101 (MQ=255)
                                                                                        ggatcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGGCCGTGGGATCGGTCGGGTCCGTcgg             >  1:5552399/4‑101 (MQ=11)
                                                                                          gtccgtgtgatcggtcggctccgttggatcggtcggggccgtggGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt           <  2:3846593/93‑1 (MQ=11)
                                                                                          atcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt           <  1:2512084/100‑1 (MQ=11)
                                                                                           tcggtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtc          >  1:100103‑M1/7‑101 (MQ=255)
                                                                                           tccgtcggGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtc          <  2:6775141/98‑1 (MQ=11)
                                                                                           gcggtcggGCCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtc          <  2:566693/100‑1 (MQ=1)
                                                                                           gccgtaggatcggtcgggtccgtggggtccgtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc          <  2:5552109/95‑1 (MQ=11)
                                                                                            cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg         >  2:63362/1‑101 (MQ=255)
                                                                                            ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTGGGGTCCGTGGGATCGGTAGAATCCGTgggggca         >  1:3868550/6‑97 (MQ=11)
                                                                                             cgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGGAGGGTCGGCGGGAccgg        >  2:6656204/5‑101 (MQ=11)
                                                                                              gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt       >  1:3530342/1‑101 (MQ=14)
                                                                                               ttggatcgttcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc      <  2:5031177/99‑1 (MQ=14)
                                                                                               tgggattggtcgggtcggtgggatcgggcggatccgtgGTATATGTGGGATCGGTTGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc      <  2:5437114/99‑1 (MQ=11)
                                                                                               tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGACGGGTCCGTGGGATCCGTCGGACACGTggggaccgtg      >  1:85551/3‑94 (MQ=11)
                                                                                               tgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc      <  2:448857/99‑1 (MQ=11)
                                                                                                gggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgggcg     >  2:6537395/2‑101 (MQ=14)
                                                                                                gggaccggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     <  2:2438620/100‑1 (MQ=14)
                                                                                                 ggatcgggcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg    <  1:4013641/101‑1 (MQ=11)
                                                                                                  gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGGtccgtcgg  >  2:5792580/1‑100 (MQ=11)
                                                                                                   atcggtcgggtccgtgggatcggtcgggtcggtgggGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGt  <  1:1350961/101‑1 (MQ=11)
                                                                                                   |                                                                                                    
GCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221859‑222058

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