Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9580A→G 3.6% 37.8 ‑6.1 ‑3.9612D233D (GAT→GACRCAP_rcc00185hypothetical protein
Rejected: Prediction has significantly lower quality scores supporting polymorphism compared to reference.
Reads supporting (aligned to +/- strand):  new base (18/4):  ref base (287/302):  total (306/306)
Fisher's exact test strand distribution p-value = 2.16e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.31e-04

Only 100 of 1272 total aligned reads displayed.

CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058
                                                                                                    |                                                                                                    
ggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                                      >  2:306980‑M1/99‑101 (MQ=255)
gatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                                      >  1:5049721‑M1/99‑101 (MQ=255)
 atccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcgggggggtccgtgggat                                                                                                     <  2:3672313‑M1/4‑1 (MQ=255)
   ccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                                                                                   <  2:4117932‑M1/6‑1 (MQ=255)
   cGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcg                                                                                                   <  2:1843220/101‑1 (MQ=255)
    ggTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgg                                                                                                  >  2:4720356/1‑101 (MQ=255)
     gtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                                 <  1:3572914‑M1/8‑1 (MQ=255)
     gtgggatcggtcgggtccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                                 <  2:1942198‑M1/8‑1 (MQ=255)
     gtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                                 <  1:2530673‑M1/8‑1 (MQ=255)
     gtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                                                 <  1:1001112‑M1/8‑1 (MQ=255)
     gTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggt                                                                                                 >  2:265517/1‑101 (MQ=255)
      tGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtc                                                                                                <  2:541756/101‑1 (MQ=255)
       gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                               <  1:2014520‑M1/10‑1 (MQ=255)
        ggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtggggtccgtgggatcgggcgggtccgtgggatcggtcgg                                                                                              <  2:2276556‑M1/11‑1 (MQ=255)
        cccGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgg                                                                                              <  1:3743926/101‑1 (MQ=255)
          atcggtcgggtccgtggggtccgtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggt                                                                                            <  1:1368692‑M1/13‑1 (MQ=255)
           tcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtc                                                                                           >  1:3430224‑M1/88‑101 (MQ=255)
             ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                                                                                         <  1:1203845‑M1/16‑1 (MQ=255)
               tcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtg                                                                                       <  2:2560681‑M1/18‑1 (MQ=255)
                cggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtaggatcggtcgggtccgtgggatcggtcgggtccgtcg                                                                                      <  2:4286105‑M1/19‑3 (MQ=255)
                    tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtc                                                                                  <  1:495672‑M1/23‑4 (MQ=255)
                       gtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggt                                                                               <  1:1547391‑M1/26‑1 (MQ=255)
                           gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtagat                                                                           <  1:3017300‑M1/30‑5 (MQ=255)
                           gatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggg                                                                           <  1:665671‑M1/30‑1 (MQ=255)
                             tCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcgggcgggtccgtggggtcggtcgggtc                                                                         >  1:927272/1‑101 (MQ=255)
                              cggtcgggtccgtggggtccgtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtcc                                                                        >  2:3879459‑M1/69‑101 (MQ=255)
                              cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                        <  2:4503379/101‑1 (MQ=255)
                               ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccg                                                                       <  1:2693706‑M1/34‑1 (MQ=255)
                                  gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGACGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                    >  2:2082855/1‑101 (MQ=255)
                                   ggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggg                                                                   <  1:1687410‑M1/38‑1 (MQ=255)
                                    gatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggga                                                                  <  2:1014478‑M1/39‑1 (MQ=255)
                                      tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatc                                                                >  2:1429218‑M1/61‑101 (MQ=255)
                                      gtggtCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                >  1:4697473/1‑101 (MQ=255)
                                        cgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggggccgtgggatcgg                                                              <  2:903777‑M1/43‑1 (MQ=255)
                                         gtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcggctccgtgggatcggtcgggtccgtgggatcggc                                                             <  2:356502‑M1/44‑2 (MQ=255)
                                         gtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggt                                                             >  2:4876776‑M1/58‑101 (MQ=255)
                                          tgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtc                                                            >  1:1433877‑M1/57‑101 (MQ=255)
                                            ggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcgg                                                          >  2:495607‑M1/55‑101 (MQ=255)
                                            ggatcggtcgggaccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgg                                                          <  1:2740576‑M1/47‑1 (MQ=255)
                                             gatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggaccgtgggctcggtcgggtccgtggggtcggtcggg                                                         <  1:4852966‑M1/48‑2 (MQ=255)
                                              gtcggtcgggtccgtgggatcggtcggggccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                                        >  1:2208617‑M1/53‑101 (MQ=255)
                                               nccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGNCGGATc                                                       <  1:384141‑M1/50‑1 (MQ=255)
                                                cggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                                      <  2:47982‑M1/51‑1 (MQ=255)
                                                    cgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                  <  1:3954818‑M1/55‑1 (MQ=255)
                                                    cgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                  <  1:4205932‑M1/55‑1 (MQ=255)
                                                     gggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTggg                                                 <  1:2282338‑M1/56‑1 (MQ=255)
                                                     ggATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                                 >  1:754079/1‑101 (MQ=255)
                                                         cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                             <  1:4410823/101‑1 (MQ=255)
                                                           gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                           <  2:4597043‑M1/62‑1 (MQ=255)
                                                            tgggatcggtccgatccgtgggatcggtcgtgtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                          >  2:841143‑M1/39‑101 (MQ=255)
                                                             gggatcggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                         <  1:214664‑M1/64‑1 (MQ=255)
                                                                  ctgtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                    >  1:2605150‑M1/33‑101 (MQ=255)
                                                                  cggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                    <  1:3090480‑M1/69‑1 (MQ=255)
                                                                  cggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                                    <  2:1252936‑M1/69‑1 (MQ=255)
                                                                       gggtccgtgtgatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               >  2:1741296‑M1/28‑101 (MQ=255)
                                                                       gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               >  2:1655832‑M1/28‑101 (MQ=255)
                                                                       gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               >  1:4914979‑M1/28‑101 (MQ=255)
                                                                       gggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                               >  1:125813‑M1/28‑101 (MQ=255)
                                                                       ggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                               <  1:191629‑M1/74‑1 (MQ=255)
                                                                        ggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGa                              <  2:3230917‑M1/75‑1 (MQ=255)
                                                                        gatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                              <  1:3665753‑M1/75‑1 (MQ=255)
                                                                         gtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                             <  1:2893923‑M1/76‑1 (MQ=255)
                                                                          tccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc                            >  2:4441012‑M1/25‑101 (MQ=255)
                                                                           gCGTCGGATGCGTGGGCTGGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                           <  2:3410818/101‑1 (MQ=255)
                                                                             gtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                         <  1:3368471‑M1/80‑1 (MQ=255)
                                                                             gcgtgccccctGGGCTAGGTGGGGTCGGTGGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggt                         <  2:1088987/91‑1 (MQ=11)
                                                                              tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGggtccgtggggtccgtg                        <  1:458538‑M1/81‑8 (MQ=255)
                                                                              tgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                        >  1:2817461‑M1/21‑101 (MQ=255)
                                                                              ggggatcggtcgggtgcgtgGGGGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                        >  1:405103‑M1/21‑101 (MQ=255)
                                                                              gcGGATCGGTCGGGTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                        <  2:3175577/100‑1 (MQ=11)
                                                                               gggatcggtcgggtccgtgGGATCGGGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                       >  1:871405‑M1/20‑101 (MQ=255)
                                                                                   tccGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc                   <  1:3242147/99‑1 (MQ=11)
                                                                                    gCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                  <  2:2840668/101‑1 (MQ=255)
                                                                                    cggtcgttccggtgggATCGGTCGGGCCCTTGGGATCGGTCGTATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                  <  2:748397/92‑1 (MQ=11)
                                                                                    cggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc                  >  1:1776158‑M1/15‑101 (MQ=255)
                                                                                    cggtcggCTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGTTCGGGTcc                  <  2:3994722/96‑1 (MQ=11)
                                                                                     ggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCg                 <  2:2376339‑M1/88‑1 (MQ=255)
                                                                                     gggggggtccgtgGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCg                 >  1:1069629‑M1/14‑101 (MQ=255)
                                                                                     cGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGGtcggtcggggccg                 >  2:4365585/1‑96 (MQ=11)
                                                                                      gtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGt                <  2:2739069‑M1/89‑1 (MQ=255)
                                                                                      gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGGGGGATCGGTCGGGTCCgg                >  1:413556/1‑100 (MQ=255)
                                                                                        cgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg              >  1:3504324‑M1/11‑101 (MQ=255)
                                                                                         gtgcccgtcggATCGGTCGGGTCCGTGGGATCGGTCGCGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg             <  1:4320586/95‑1 (MQ=11)
                                                                                         ggatcggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTGCGTcgg             >  1:1844097/4‑101 (MQ=1)
                                                                                           gtccgggggatcgttcgagtccgtggcatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt           <  2:636168/94‑1 (MQ=11)
                                                                                           gctcgggggatctgtcggctccgtgggatcggtcgggcccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt           <  2:2460813/94‑1 (MQ=11)
                                                                                             ccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTGGGGTCCGTGGGTTCGGTCGGGCCGGTGGGAtcg         >  1:4452698/3‑101 (MQ=11)
                                                                                             ccgggggatgggtcgggtccgtggggtccgtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg         <  1:4769152/96‑1 (MQ=11)
                                                                                                tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTGGGGCTCGTGGGCTCGGTCGGCTTCGTGGGAgcgggt      >  1:552213/3‑99 (MQ=11)
                                                                                                tgggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtc      <  1:2092250/99‑1 (MQ=11)
                                                                                                gcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtc      <  2:3769562/100‑1 (MQ=14)
                                                                                                gcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtc      <  2:1491822/100‑1 (MQ=11)
                                                                                                 tggatcggttgggcccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     <  2:2940034/100‑1 (MQ=11)
                                                                                                 gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGCCGGGTCCGTGGGATCGGTCGGGCCCNTGGGATCGGCCGGGTCCGCGGGGccggtcg     >  1:626167/2‑101 (MQ=11)
                                                                                                 gggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcg     >  1:3179072/2‑101 (MQ=11)
                                                                                                 gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     <  1:239488‑M1/100‑1 (MQ=255)
                                                                                                 gGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg     <  2:1188674‑M1/100‑1 (MQ=255)
                                                                                                  ggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGAtcggtcgg    >  2:2184969/1‑101 (MQ=14)
                                                                                                  ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGGGAGACCGGGCGGGCCCGTGGGATCGGTCGAGTCCGTGGGAgcggtcgg    >  2:4182753/1‑101 (MQ=11)
                                                                                                    atcggtcggatccgtgggatcggtcgggtccgtgggATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGGGGGAtcggtcggat  >  2:3495704/1‑99 (MQ=11)
                                                                                                    |                                                                                                    
CGCCGGTGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGT  >  NC_014034/221858‑222058

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