Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9640C→G 1.5% 26.5 ‑1.3 ‑2.21348P231P (CCG→CCCRCAP_rcc00185hypothetical protein
Rejected: Prediction has significantly lower quality scores supporting polymorphism compared to reference.
Reads supporting (aligned to +/- strand):  new base (13/8):  ref base (642/684):  total (655/693)
Fisher's exact test strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.85e-03

Only 100 of 3091 total aligned reads displayed.

TGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAG  >  NC_014034/221864‑222064
                                                                                                    |                                                                                                    
 gggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                     >  1:4416835‑M1/92‑101 (MQ=255)
 gggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                     <  1:1996212‑M1/10‑1 (MQ=255)
  ggatgggtcggatccgtggggtcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                                                                                                    >  2:2593899‑M1/91‑101 (MQ=255)
  ggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                                                                                                    <  2:377164‑M1/11‑1 (MQ=255)
     gACACCGTCGCCGAGGCTTCGGTGTTGGTGGGCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtc                                                                                                 <  2:678497/101‑1 (MQ=255)
          cgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                                            <  1:5250155‑M1/19‑1 (MQ=255)
            ggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggg                                                                                          <  1:4792462‑M1/21‑2 (MQ=255)
            gatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                          <  1:7121482‑M1/21‑1 (MQ=255)
             cGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggat                                                                                         <  2:3963544/101‑1 (MQ=255)
                cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccg                                                                                      >  2:2786039‑M1/77‑96 (MQ=255)
                 gacgCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCACGGTCggatcggtcgggtccgtgggatcggt                                                                                     >  1:5869275/4‑101 (MQ=255)
                    ggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgg                                                                                  <  1:6207657‑M1/29‑1 (MQ=255)
                      ttCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggt                                                                                <  2:633877/101‑1 (MQ=255)
                        tggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                                                              >  1:5389757‑M1/69‑101 (MQ=255)
                        cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                              <  2:5360299/101‑1 (MQ=255)
                         ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                                                                             <  1:243592‑M1/34‑1 (MQ=255)
                         ggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccg                                                                             <  2:9055818‑M1/34‑1 (MQ=255)
                           tgtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtccgtg                                                                           >  2:5532364/1‑101 (MQ=255)
                                  cgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcgg                                                                    <  1:6852172‑M1/43‑1 (MQ=255)
                                     gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                                                                 <  1:4079700‑M1/46‑1 (MQ=255)
                                     cggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                 <  1:4531407‑M1/46‑1 (MQ=255)
                                        gtcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAt                                                              >  2:5844549‑M1/53‑101 (MQ=255)
                                        atcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGAt                                                              <  1:2787113‑M1/49‑1 (MQ=255)
                                         tcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATc                                                             <  1:9039897‑M1/50‑1 (MQ=255)
                                          cggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGCTcc                                                            <  1:2122663‑M1/51‑1 (MQ=255)
                                           ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGATCGGATCCg                                                           <  2:5967596‑M1/52‑1 (MQ=255)
                                           ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCg                                                           >  1:2300342‑M1/50‑101 (MQ=255)
                                           ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                           <  1:3330911‑M1/52‑1 (MQ=255)
                                           ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                           >  2:7293395‑M1/50‑101 (MQ=255)
                                               ggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCAGTggg                                                       <  1:7450249‑M1/56‑1 (MQ=255)
                                                      tgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                <  2:4502957‑M1/63‑1 (MQ=255)
                                                      tgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc                                                >  1:4695134‑M1/39‑101 (MQ=255)
                                                      gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtg                                                >  2:4930028/1‑100 (MQ=255)
                                                      gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                >  1:3665841/1‑101 (MQ=255)
                                                       gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                               >  2:4779203‑M1/38‑101 (MQ=255)
                                                       gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                               <  1:8879027‑M1/64‑1 (MQ=255)
                                                        ggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                              <  1:6032070‑M1/65‑1 (MQ=255)
                                                        ggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                              <  2:1892943‑M1/65‑1 (MQ=255)
                                                         cATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAGCGGTCggg                                             >  2:2424575/1‑101 (MQ=255)
                                                          atcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGgatcggtcggac                                            <  1:1354683‑M1/67‑3 (MQ=255)
                                                           tcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTc                                           <  1:3167436‑M1/68‑1 (MQ=255)
                                                            cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                          >  1:8126496‑M1/33‑101 (MQ=255)
                                                            cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                          >  1:3000604‑M1/33‑101 (MQ=255)
                                                            cggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGtcg                                          <  2:4315092‑M1/69‑2 (MQ=255)
                                                            cggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                          >  1:7947811‑M1/33‑101 (MQ=255)
                                                            cggtcggatccgtgggagcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                          >  2:8226089‑M1/33‑101 (MQ=255)
                                                             ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                         <  2:6497843‑M1/70‑1 (MQ=255)
                                                              aCGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCgc                                        >  1:5061624/1‑100 (MQ=255)
                                                                cggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                      <  1:3934579‑M1/73‑1 (MQ=255)
                                                                 ggatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                     >  2:6652716‑M1/28‑101 (MQ=255)
                                                                 gTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                     >  1:8552879/1‑101 (MQ=255)
                                                                  gatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                    >  2:8360518‑M1/27‑101 (MQ=255)
                                                                   gtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                                   >  2:7503517‑M1/26‑101 (MQ=255)
                                                                   ggccgtgggatccgtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                   >  2:4109857‑M1/26‑101 (MQ=255)
                                                                     ccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                                 >  2:6928556‑M1/24‑101 (MQ=255)
                                                                      cgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                <  1:7255796‑M1/79‑1 (MQ=255)
                                                                      cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                >  2:768692‑M1/23‑101 (MQ=255)
                                                                      cGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                <  1:778936/101‑1 (MQ=255)
                                                                       gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                               >  1:8267480/1‑101 (MQ=255)
                                                                        tgggagtggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                              >  1:961469‑M1/21‑101 (MQ=255)
                                                                        tgcgATGCGTCGCCTCCCTGGCATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                              <  2:8684869/98‑1 (MQ=11)
                                                                             tcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGggatcggtcgggtc                         <  1:1789627‑M1/86‑4 (MQ=255)
                                                                              cggtcgggtccgtgGGATCGGTCGGATCTGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                        >  1:6389639‑M1/15‑101 (MQ=255)
                                                                               ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                       >  1:7650089‑M1/14‑101 (MQ=255)
                                                                                 cgatcTCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtc                     <  2:6170414/97‑2 (MQ=11)
                                                                                   gggtcgggcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTcgg                   <  1:1000386/98‑1 (MQ=11)
                                                                                   gggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCGGTGGGGTCGGTCGGATCCGTggg                   <  2:1428570‑M1/92‑1 (MQ=255)
                                                                                   gggtccgggGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                   >  1:2612961‑M1/10‑101 (MQ=255)
                                                                                     gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                 <  2:3230299‑M1/94‑1 (MQ=255)
                                                                                     gtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGCCGGAt                 >  1:8371923/5‑101 (MQ=11)
                                                                                      tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGACGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGggacc                <  2:5641603‑M1/95‑3 (MQ=255)
                                                                                       ccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGATCCGTCGGATCGGGGGGGTCCGTGGGATCGGGCGGGTCCGTGGGAtcg               >  2:1046960/6‑101 (MQ=11)
                                                                                        ggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTCGGAtcgg              >  2:2246658/1‑101 (MQ=11)
                                                                                        ggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg              >  1:6966123/1‑101 (MQ=255)
                                                                                         gtggggtccgtcggatccgtgggatcggtcggatccgtgggatcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtcggt             >  2:465939/10‑101 (MQ=11)
                                                                                          tggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtc            <  2:1739039/93‑1 (MQ=14)
                                                                                           gggatcggtcgggtccgtcggatcggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCCTGGGAtcgggcg           >  2:8493712/2‑101 (MQ=11)
                                                                                           gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATGGGGCGGGTCCGGGGGAtgcgtcg           >  2:3057508/2‑101 (MQ=11)
                                                                                            ggatcggtcgggtccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTGGGGTCGGTCGGGTCCGTGGGAtcggtcgg          >  1:5541720/1‑101 (MQ=11)
                                                                                              gtcggtcgggtccgtgggatcggtcgggtccgtcggaTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGt        <  2:791780/100‑1 (MQ=11)
                                                                                              gtcggtcgggtccgggggagcggtcgggtccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggat        >  1:6678659‑M2/2‑88 (MQ=255)
                                                                                              gtccgtcgggtccgtgggatcggtcggatccgtgggctcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGCTCGGTCGGATCCGTGGGATCGGTCGGGt        <  1:6511973/97‑1 (MQ=14)
                                                                                              atcggtgggatccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTCGGATCGGCCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCgggc        >  1:1650421/1‑100 (MQ=11)
                                                                                              atcggtcgggtccgtgggatcggtcgggtccgtgggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggat        >  1:1510716‑M2/1‑88 (MQ=255)
                                                                                               gcggtcggagccgtgggatcgggcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc       <  1:3501890/100‑1 (MQ=11)
                                                                                               gccgtcggctccgtgggatcggtcggatccgtgcgatggGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGATCCGTGGGATCGGTCGGGTc       <  2:3622128/98‑1 (MQ=14)
                                                                                                tggtcgtgtcctcgggatcgggcggtgcctggagatcGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc      <  1:4242956/100‑1 (MQ=11)
                                                                                                gggtcgggttcgtggggtcagtcgggtccttgggttcCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc      <  1:3816674/100‑1 (MQ=14)
                                                                                                cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc      <  2:209718/101‑1 (MQ=11)
                                                                                                cgggcgggtccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATcc      <  1:871610/101‑1 (MQ=14)
                                                                                                ccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGCGCCGTCGGGCCGGTCGAGTcc      >  2:3440298/3‑101 (MQ=11)
                                                                                                 ggtcgcgtccgtggggtccgtcggatccgtgggatcGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg     <  2:1580533/101‑1 (MQ=11)
                                                                                                  gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTGGGGGCGGTCGGATcggc    >  1:4917920/1‑98 (MQ=11)
                                                                                                  gtcggggccgtggggtccgtcggatccgtgggatcgGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt    <  1:5650767/101‑1 (MQ=11)
                                                                                                  gtcggatccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGGAGGGATCCgg    >  2:8866350/1‑100 (MQ=11)
                                                                                                   gcgggtccgtgggatcggtcgggtccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtgggatcggtc   >  2:808803‑M2/2‑83 (MQ=255)
                                                                                                    cgggtccgtgggatcggtcgggtccgtggggtccgtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgg  <  1:1170498‑M2/101‑3 (MQ=255)
                                                                                                    cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGGGGATCGGACGCGACCgggg  >  2:1226627/1‑98 (MQ=11)
                                                                                                    cgggtccgtgggatcggtcgggtccgtgggatcggtCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg  <  2:5201140‑M2/101‑3 (MQ=255)
                                                                                                    cggatccgtgggatcggtcggatccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtgg  >  2:7130091/1‑99 (MQ=11)
                                                                                                    |                                                                                                    
TGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAG  >  NC_014034/221864‑222064

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