Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,979 G→C 4.5% T226T (ACC→ACG)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9790G→C 4.5% 50.0 ‑0.1 0.0702T226T (ACC→ACG) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (16/18):  ref base (327/341):  total (343/359)
Fisher's exact test strand distribution p-value = 8.62e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 1915 total aligned reads displayed.

CCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCC  >  NC_014034/221879‑222087
                                                                                                    |                                                                                                            
ccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                                              <  2:2696211‑M1/24‑1 (MQ=255)
ccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcg                                                                                                              >  1:3547190‑M1/78‑101 (MQ=255)
    gggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcg                                                                                                          >  1:1953531‑M1/74‑101 (MQ=255)
         cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                                                     <  2:1797298/101‑1 (MQ=255)
            tcgggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtg                                                                                                  <  1:120866‑M1/36‑1 (MQ=255)
              gggtccgtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGatcggtcgggtccgtgggatcggtcgggtccgtggg                                                                                                <  1:4210960‑M1/38‑1 (MQ=255)
               ggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGAtcggtcgggtccgtgggatcggtcgggtccgtggga                                                                                               >  1:4482020‑M1/63‑101 (MQ=255)
                 tccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtggggtc                                                                                             <  2:1587510‑M1/41‑4 (MQ=255)
                 tccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCggtcgggtccgtgggatcggtcgggtccgtgggatc                                                                                             >  1:932886‑M1/61‑101 (MQ=255)
                  ccgtgggatcggtcgggtccgtggggtccggcggatccgtgggatcggtcggatccgtgGGATCGgtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                                            >  2:5044442‑M1/60‑101 (MQ=255)
                   cgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGtcggatccgtgggatcggtcggatccgtgggatcgg                                                                                           <  1:5556768‑M1/43‑1 (MQ=255)
                    gtgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTcgggtccgtgggatcggtcgggtcggtgggatcggt                                                                                          <  2:4013265‑M1/44‑1 (MQ=255)
                     tgggatcggtcgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                         <  1:4796944‑M1/45‑1 (MQ=255)
                      gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGggtccgtgggatcggtcgggtccgtggggtcggtcg                                                                                        <  1:4353032‑M1/46‑1 (MQ=255)
                      gggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                        <  1:4894484‑M1/46‑1 (MQ=255)
                      cAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGggtccgtgggatcggtcgggtccgtgggatcggtcg                                                                                        <  1:2335417/101‑1 (MQ=255)
                       aGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGgtccgtgggatcggtcgggtccgtgggatcggtcgg                                                                                       >  2:4449078/1‑101 (MQ=255)
                         cGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                                                                     <  2:1165835/101‑1 (MQ=255)
                           tCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                                   >  2:5750087/1‑101 (MQ=255)
                            cgtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCg                                                                                  <  2:1987908‑M1/52‑1 (MQ=255)
                              tcgggtccgtgggatccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                                                                >  2:5296654‑M1/48‑101 (MQ=255)
                               cgggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                                               <  1:1853199‑M1/55‑1 (MQ=255)
                                 ggtccgtggggtccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                                                             >  1:3747331‑M1/45‑101 (MQ=255)
                                  ggccgtgtgatcggtcggaaccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                                            >  1:364748‑M1/44‑101 (MQ=255)
                                  aTCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAt                                                                            >  1:671945/1‑101 (MQ=255)
                                     cgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                                         >  2:4568525‑M1/41‑101 (MQ=255)
                                        gggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGggtcggtgg                                                                      <  1:2246451‑M1/64‑3 (MQ=255)
                                         ggatcggtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTTCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgg                                                                     <  1:4985288‑M1/65‑1 (MQ=255)
                                             ccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGggtac                                                                 <  2:4946405‑M1/69‑3 (MQ=255)
                                             ccgtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                                 >  2:4354365‑M1/33‑101 (MQ=255)
                                               gtcggatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCcgg                                                               <  2:3303793‑M1/71‑2 (MQ=255)
                                                 cggatccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                             <  1:4305897‑M1/73‑1 (MQ=255)
                                                    gggCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                                                          <  2:4766403/101‑1 (MQ=255)
                                                          gggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGgatcggtgg                                                    <  2:1362221‑M1/82‑3 (MQ=255)
                                                           ggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                   <  1:2024051‑M1/83‑1 (MQ=255)
                                                             gtGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGAt                                                 <  1:5781073/100‑1 (MQ=255)
                                                              tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                <  1:5371721/101‑1 (MQ=255)
                                                                   gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                           >  1:396471/1‑99 (MQ=11)
                                                                    ggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                                          <  2:259139‑M1/92‑1 (MQ=255)
                                                                    ggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                                          >  1:1008891‑M1/10‑101 (MQ=255)
                                                                     gatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                         >  1:4666411‑M1/9‑101 (MQ=255)
                                                                      gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGACCCGTGGGAt                                        >  1:1636694/8‑101 (MQ=11)
                                                                      gtccgtcggATCGGTCGGGTCCTTCCGATCGGTCGGTTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                        <  2:3201835/97‑1 (MQ=11)
                                                                      gtccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAt                                        <  2:773823/97‑1 (MQ=11)
                                                                       tccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGatg                                       <  2:3679055‑M1/95‑2 (MQ=255)
                                                                        ccgtcggATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGAtcg                                      >  2:5504153/3‑101 (MQ=11)
                                                                        ccgtccgATCGTGCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                      <  2:3929582/99‑1 (MQ=11)
                                                                         cgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGATCGGATCCGTGGAGTCGGTCGGATCCGCGGGATCAATCGGATCCGTGGGAtcgg                                     >  2:649835/5‑101 (MQ=11)
                                                                         cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcgg                                     <  1:796675‑M1/97‑1 (MQ=255)
                                                                          gtggggtcggtcgggtccgtgggatcggtcgggtccgtggggTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                    <  1:4575861/95‑1 (MQ=11)
                                                                          gtgggatcggtcgggtccgttggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                                    <  1:2587461/98‑1 (MQ=11)
                                                                          gtgGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggt                                    >  1:4035672‑M1/4‑101 (MQ=255)
                                                                          gtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggt                                    <  1:5459062‑M1/98‑1 (MQ=255)
                                                                            gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                  >  1:339601‑M1/2‑101 (MQ=255)
                                                                            gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                  >  2:4210700‑M1/2‑101 (MQ=255)
                                                                             ggctcgggcgggtccgttggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtcgg                                 <  2:4092211/98‑1 (MQ=11)
                                                                             ggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGGtccgtggg                                 >  2:3830562/1‑101 (MQ=11)
                                                                               atcggtcgggtccgtgggatcggtcgggtccgtgggGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGATCCGTGGCGtgggtcggaa                               >  2:2313829/1‑99 (MQ=14)
                                                                               atcggtcggatccgtgggatcggtcgggtccgtgggATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGat                               >  2:2212193‑M2/1‑99 (MQ=255)
                                                                                tgggtcgggtccgtgggatcggtcgggtccgtcggatcGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                              <  1:2805255/99‑1 (MQ=11)
                                                                                tccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTc                              <  1:3194822/98‑1 (MQ=14)
                                                                                 cttccgagctcgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATcc                             <  2:1699587/91‑1 (MQ=11)
                                                                                 cggtcgggtccgggggatcgggcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                             <  1:2692952/101‑1 (MQ=11)
                                                                                 cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTGGGGTCCGTGGGGTCGGCGGGGTcc                             >  2:2931728/1‑101 (MQ=11)
                                                                                 cggtcggatccgtgggatcggtcgggtccgtgggatCGGCCGGGTCCGTGGGGTCGGCCGGGTCCGCGGGGGCGGTCGGGTCGGTGGGGTAGGTCggggca                             >  2:853018/1‑98 (MQ=2)
                                                                                 cgggcgggtccgttggatctgtcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                             <  1:5657544/101‑1 (MQ=11)
                                                                                  ggtcgggtcggtgggatcggtcgggtccgtggggtcCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                            <  1:5065187/101‑1 (MQ=14)
                                                                                  cgtcgggtccgtgggatcggtcgggtccgtgggatcggTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccg                            <  2:542114‑M2/100‑6 (MQ=255)
                                                                                   gtcgggtccgtgggatcggtcgggtccgtgggatcgGTCGGGTCCGTGGGGTCCGTCGGATCCGGGGGATCGGGCGGAGCCGTGGGTCCGGACGGGACCgg                           >  2:3556323/1‑100 (MQ=11)
                                                                                    tggggtgcgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGGTCCATCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtg                          <  1:1545160/99‑2 (MQ=11)
                                                                                    tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGATCCGTGGGATCGGTCGGGTCCATGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtg                          >  1:527887‑M2/1‑94 (MQ=255)
                                                                                    tcgggtccgggggatctgtcgggtccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGtg                          <  2:1304805/101‑2 (MQ=11)
                                                                                    tcgggtccggcggatcggtcgtgtccgtggggtccgTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAACGGTCGGATCCGtg                          <  2:5647854/101‑2 (MQ=11)
                                                                                    tcggatccgtgggatcggtcgggtccgtgggatcggTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGGGGGATCGGTCGGATGCGtg                          >  2:303852/1‑100 (MQ=11)
                                                                                     cggggccgtgggatcggtcgggtccgtcggatcggtGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGtgg                         <  1:3866307/101‑3 (MQ=11)
                                                                                      gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCggggg                        >  1:3098320/1‑97 (MQ=11)
                                                                                       ggtccgtgggatcggtcgggtccgtggggtcagtcgGATCCGTGGGGTCGGTCGGATCCGTGGGCTCGGTAGGATCCGTGGGATCGGTAGGGTCCGTTGGa                       >  1:5157758/1‑101 (MQ=14)
                                                                                       ggtccgtgggatcggtcgggtccgtgggatcggtcgGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgagatc                       <  1:1719951‑M2/101‑11 (MQ=255)
                                                                                       gctctgtgggatcggtctggtcagtggggtccgtcggCTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                       <  2:1057113/96‑1 (MQ=11)
                                                                                        tcttccgatctcggtcggagccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                      <  1:4113566/91‑1 (MQ=14)
                                                                                        tctcgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAt                      <  2:4665892/98‑1 (MQ=11)
                                                                                        gtgttgggggtcggtggggtcggtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGgggt                      <  1:145524/95‑3 (MQ=1)
                                                                                         tccgtgggatcggtcgggtccgtgggatgggtcggGGCAGGTCGAGTGGGCGGGGCGGTGGGATCGCTGGGGTCGGTGGGTTCTGACGGGCCCGTGGGAtg                     >  2:486524/1‑100 (MQ=2)
                                                                                         tccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatc                     >  1:1299333‑M2/1‑89 (MQ=255)
                                                                                          ccgtggggttggtccgatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcg                    <  2:2802558/101‑1 (MQ=11)
                                                                                          ccgtgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcg                    >  1:1105026‑M2/1‑88 (MQ=255)
                                                                                           cgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTTGGGTCCGTGGGGtcag                   >  2:2086487/1‑99 (MQ=11)
                                                                                           cgtgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGGCGGATCGGTCGGGACGGTGGGAtagg                   >  2:3399807/1‑101 (MQ=11)
                                                                                             tgggatcggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGAACGTCGGATCGGCCGGGATCGTGGGAtcggtc                 >  1:2949978/1‑101 (MQ=11)
                                                                                              gggatcggtcgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                <  1:2474641/101‑1 (MQ=11)
                                                                                               ggatcgttcgggtccgtggaatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCCTGGGATCGGTCGGGTCCGTGGGAtcggtcgg               <  1:3263142/101‑1 (MQ=14)
                                                                                               ggatcggtcgggtccgggggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcggtcgg               <  1:2139748/101‑1 (MQ=11)
                                                                                                gatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGgtccgtgggatcggtcggg              <  1:863656‑M2/101‑20 (MQ=255)
                                                                                                 atcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggt             >  2:1502539/1‑101 (MQ=11)
                                                                                                  tcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGatccgtgggatcggtgggatc            <  2:3113019‑M2/101‑22 (MQ=255)
                                                                                                   gggtcgggtgcgtgggatcggtcggGTGCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtcc           <  2:2491112/100‑1 (MQ=11)
                                                                                                   cggtcggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtcc           <  2:5146147/101‑1 (MQ=255)
                                                                                                   cgggcgggtccggggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcggatcc           <  1:1362309/101‑1 (MQ=11)
                                                                                                    ggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTGCGTGGGGTCGGTGGGGTCCGTGGGATCGgtccggtccg          >  2:4297936/1‑101 (MQ=11)
                                                                                                    ggtcgggtccgtcggatcggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAGCGgtcgggtccg          >  2:1922789/1‑101 (MQ=11)
                                                                                                    |                                                                                                            
CCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCC  >  NC_014034/221879‑222087

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