Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,973 G→A 1.4% P228P (CCC→CCT)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9730G→A 1.4% 11.2 ‑2.2 ‑0.0721P228P (CCC→CCT) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/8):  ref base (337/372):  total (339/382)
Fisher's exact test strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01

Only 100 of 1673 total aligned reads displayed.

CCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCG  >  NC_014034/221873‑222073
                                                                                                    |                                                                                                    
 cgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgg                                                                                                     <  1:3995759‑M1/7‑1 (MQ=255)
  gggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggg                                                                                                    <  1:2614091‑M1/8‑1 (MQ=255)
   gatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtggga                                                                                                   <  1:1392771‑M1/9‑1 (MQ=255)
     gCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatc                                                                                                 <  2:2702923/101‑1 (MQ=255)
       cgggggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcgg                                                                                               >  1:3254382‑M1/89‑101 (MQ=255)
         tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc                                                                                             >  2:2501852‑M1/87‑101 (MQ=255)
               cGGTGTTGGAGGTCAGCGTCGCCGGATCGCGCTCCTCGACGGTGGGCGTCGGATGCGTGGGCTCGGCCggatcggtcgggtccgtgggatcggtggggtcc                                                                                       >  1:380878/1‑101 (MQ=39)
                ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCg                                                                                      <  1:1979013‑M1/22‑1 (MQ=255)
                 gtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGt                                                                                     <  1:2209533‑M1/23‑1 (MQ=255)
                  gggggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTg                                                                                    >  1:2845022‑M1/78‑101 (MQ=255)
                   cggatccgtgggatcggtcggatccgtgggagcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTgg                                                                                   >  1:1273481‑M1/77‑101 (MQ=255)
                     gatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGa                                                                                 >  1:2039032‑M1/75‑101 (MQ=255)
                       tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTCGGAtc                                                                               <  1:3892472‑M1/29‑1 (MQ=255)
                       tccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                               >  2:3599541‑M1/73‑101 (MQ=255)
                          gtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggaTCCGTAGGATCGGTCGGGTCCGTAGGAtcggt                                                                            >  1:4545452‑M1/70‑101 (MQ=255)
                           tgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCggatccgtgggatcggtcggatccgtgggatcggtc                                                                           <  2:1002993‑M1/39‑1 (MQ=255)
                            gggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggaTCCGTAGGATCGGTCGGGTCCGTAGGAtcggtcg                                                                          <  1:944330‑M1/34‑1 (MQ=255)
                              gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtcgGTCGggtccgtgggatcggtcgggtccgtggggtcggtcggg                                                                        <  1:3003386‑M1/42‑2 (MQ=255)
                               atcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGgatcggtcgggt                                                                       >  1:1917072‑M1/65‑99 (MQ=255)
                                 tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcgcccac                                                                     <  2:856123/101‑6 (MQ=255)
                                 cggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGtgggatcggtcgagatc                                                                     <  1:2543599‑M1/39‑6 (MQ=255)
                                  ggtcgtatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCg                                                                    >  1:1708213‑M1/62‑101 (MQ=255)
                                  ggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCg                                                                    >  1:2498753‑M1/56‑101 (MQ=255)
                                    tcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTg                                                                  <  2:3694509‑M1/42‑1 (MQ=255)
                                       ggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                                               <  1:3183780‑M1/45‑1 (MQ=255)
                                       gatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtgggc                                                               <  1:1002062‑M1/51‑2 (MQ=255)
                                          tgggcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                            <  2:645868/98‑1 (MQ=255)
                                          cgcgCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                            >  2:3201915/1‑101 (MQ=255)
                                            cgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                                          >  2:4286196/1‑101 (MQ=255)
                                             tgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                         >  1:3082084‑M1/51‑101 (MQ=255)
                                             tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                                         >  1:1506196‑M1/51‑101 (MQ=255)
                                              gggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                        >  2:4508488‑M1/50‑101 (MQ=255)
                                                 atcggtcgggtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                                                     >  1:4477508‑M1/41‑101 (MQ=255)
                                                   cgggcgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                   >  1:135069‑M1/45‑101 (MQ=255)
                                                    gACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTNCg                                                  >  1:3629149/1‑101 (MQ=255)
                                                      tcgggtccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCgta                                                >  2:3406057‑M1/42‑100 (MQ=255)
                                                      cGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                                <  1:3359561/101‑1 (MQ=255)
                                                       cgggtccgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                               <  1:1063949‑M1/61‑1 (MQ=255)
                                                       cgggtccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCcgtag                                               <  1:1130312‑M1/61‑3 (MQ=255)
                                                       cggatccgtgggatcggtcggatccgtgggatcggtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCcgtag                                               <  1:2707897‑M1/61‑3 (MQ=255)
                                                          gtccgtgggatcggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                            <  2:2987936‑M1/70‑1 (MQ=255)
                                                           gccgtcggatcggtcgggtccgtgggatccGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtc                                           >  2:278931‑M1/31‑101 (MQ=255)
                                                             cgtgggatcggtcggatccgtgggatcggtcgggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                         <  1:3528731‑M1/67‑1 (MQ=255)
                                                                    ctggcgggGTCGGTGGGGTCGGCCGGGGCCGTGGGATAGGTGGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc                                  <  1:344523/96‑1 (MQ=11)
                                                                     cggtcgggtccgtggggtccGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtcg                                 <  1:856580‑M1/81‑5 (MQ=255)
                                                                       gggggCTCGGTCGGATCGGTCGGGTCCGTGGGAGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGTATCCGt                               <  2:4600350/99‑1 (MQ=255)
                                                                          gggtcggtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                            <  1:2375130/98‑1 (MQ=11)
                                                                          ggCGCGATCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGATCCGTGgg                            <  2:1905128‑M2/101‑3 (MQ=255)
                                                                            gttcgtgtgatgggtcgggtccttgggatcggtcgggtccgtggGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                          <  2:1773873/93‑1 (MQ=11)
                                                                            gtccgtgggatcggtcggaTGCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAt                          >  2:1980251‑M1/20‑101 (MQ=255)
                                                                            atccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGTTCCGTGGGATCGGTCGGTTCCCTGGGATCGGACGGATACGTggatt                          >  1:547116/8‑98 (MQ=11)
                                                                             tgcttgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtc                         <  1:625094/95‑1 (MQ=14)
                                                                              cggtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcg                        <  1:4147902/101‑1 (MQ=255)
                                                                               ggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcgg                       <  2:4356612/101‑1 (MQ=11)
                                                                               cgtgggatcgttcgggtccgtgggagcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                       <  2:419704/97‑1 (MQ=11)
                                                                               cgtgggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                       >  2:1591789‑M1/17‑101 (MQ=255)
                                                                               cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGtccgagccg                       >  2:3946207/2‑93 (MQ=11)
                                                                                gtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt                      <  1:2091410/101‑1 (MQ=11)
                                                                                gtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGTtcggg                      >  2:3089813/1‑100 (MQ=14)
                                                                                 tgggatcggtcgggtccgtgggatcggtcgggtccgtgggATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGgggtccgtg                     >  2:403493‑M2/3‑92 (MQ=255)
                                                                                  ggggtccgtcggaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                    >  1:2157053‑M1/14‑101 (MQ=255)
                                                                                  gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGGTCCGTGGGAtcgggcg                    >  2:3266976/2‑101 (MQ=11)
                                                                                  gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAgcggtcg                    >  2:941643/2‑101 (MQ=14)
                                                                                  gggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTCGGGCCCGTGGGAgcgggcg                    >  1:1615298/2‑101 (MQ=11)
                                                                                  gggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCGGGGGGAtcggtcg                    <  2:4707864‑M1/88‑1 (MQ=255)
                                                                                   gggtcggtcgggtccgtgggatcggtcggtgccgtcggaTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                   <  2:2103814/98‑1 (MQ=11)
                                                                                   gggtcgggcgggtccgtgggatcggtcgggtccgtcggaTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggtcgg                   <  1:1819547/98‑1 (MQ=11)
                                                                                   gggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtcgg                   <  2:2329151/95‑1 (MQ=14)
                                                                                   ggatcggtcggatccgtgggatcggtcggatccgtaGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGAtcggtcgg                   >  1:4702159/1‑101 (MQ=11)
                                                                                   ggatcggtcggaTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGggatcggtggc                   <  2:2237110‑M1/89‑4 (MQ=255)
                                                                                    ggacgttcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgga                  <  2:1218758/95‑2 (MQ=14)
                                                                                     gtccggcgggtccgtgggatcggtcgggtccgtcggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGt                 <  1:4396683/95‑1 (MQ=14)
                                                                                       tggtcggatctgtgggatcggtcggacccgaaggctcGGTCGGGTCCGTAAGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc               <  2:734405/100‑1 (MQ=11)
                                                                                       cggtcgggtccgtgggatcggtcgggtccgtgggatCGGTCGTATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGAcc               >  2:2418838/1‑101 (MQ=14)
                                                                                        gttcgggtccgtgggattggtcgggtccgtggggtctgTCGGTGCCGGGGGATCGGTCGGCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg              <  2:1328160/99‑1 (MQ=14)
                                                                                         gtcgggtccgtgggatcggtcgggtccgtcggatcgGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGCCCGTCGGATCGGTCGGGGCCgg             >  2:1141118/1‑100 (MQ=11)
                                                                                         gtcgggtccgggggatcgggcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGt             <  1:230512/101‑1 (MQ=14)
                                                                                         gtcggaTCCGTGGGATCGGTCCGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt             >  2:413754‑M1/7‑101 (MQ=255)
                                                                                          tctgatccgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtc            <  2:2039467/96‑2 (MQ=11)
                                                                                          tcgggtccgtgggatcggtcgggtccgtgggagcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGtg            <  1:4328687/101‑2 (MQ=14)
                                                                                          tcgggtccgtgggatcggtcgggtccgtcggatcggTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATGGGTCGGGTTCGtg            >  2:3825883/1‑100 (MQ=11)
                                                                                           cgggtccgtcggatcggtcgggtccgtcggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGtgg           <  2:3492340/101‑3 (MQ=14)
                                                                                            gggtccgtgggatcggtcggatccgtgggatcggtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtagg          <  1:2898212‑M2/101‑10 (MQ=255)
                                                                                            gggtccgtgggatcgctcgggtccgtgggatcggtcGGGGCCGTGGGATCCGTCGGATCCGTGGGATCGGTCGGATCCCTGGGAACGGGCGGGTccaggag          >  1:762052/1‑96 (MQ=11)
                                                                                            ggatccgtggggtccgtcgggtccgtgggatcggtcggGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGTCCGTGgg          >  1:4630226/4‑101 (MQ=11)
                                                                                            ggatccgtgggatcggtcgggtccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtggg          <  2:1404932‑M2/101‑10 (MQ=255)
                                                                                            ggatccgtgggatcggtcgggtccgtgggatcggtcGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtggg          >  2:4410290‑M2/1‑92 (MQ=255)
                                                                                             gaTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGa         >  1:3717328‑M1/3‑101 (MQ=255)
                                                                                              gtccgtgggatcggtctggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggat        >  1:682982‑M2/1‑90 (MQ=255)
                                                                                              gtccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGatccgtgggat        <  1:428134‑M2/101‑12 (MQ=255)
                                                                                              gtccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggat        <  2:2591720‑M2/101‑12 (MQ=255)
                                                                                               tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGCGCGGATCCGTGGGATCGGTCGGATCCGTGGGAtc       >  1:4073858/1‑101 (MQ=11)
                                                                                               tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatc       >  2:1697116‑M2/1‑89 (MQ=255)
                                                                                               tccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGatccgtgggatc       >  2:1060059‑M2/1‑89 (MQ=255)
                                                                                               gccgtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGAtc       <  1:1190850/100‑1 (MQ=14)
                                                                                                ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGCGGGAtcg      >  1:419748/1‑101 (MQ=11)
                                                                                                ccgtgggatcggtcggatccgtaggatcgggcggGTCCGTAGGATCGGTCGGGTCGGTCGGATCGGGAGCATCCGCGGGGTGGTGCGGGTGGGTAGGAtcg      >  2:1475871/1‑101 (MQ=2)
                                                                                                  gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGgtccgtgggatcggt    <  2:2898035‑M2/101‑16 (MQ=255)
                                                                                                   tgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGCtcgggc   >  2:2160008/1‑99 (MQ=11)
                                                                                                   tgggatcggtcggggccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc   <  2:561544/101‑1 (MQ=14)
                                                                                                    |                                                                                                    
CCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCG  >  NC_014034/221873‑222073

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