Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,081 C→G 24.5% P192P (CCG→CCC)  RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0810C→G24.5% 251.2 ‑2.8 ‑0.4397P192P (CCG→CCC) RCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (40/57):  ref base (156/143):  total (196/201)
Fisher's exact test strand distribution p-value = 6.29e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.86e-01

Only 100 of 589 total aligned reads displayed.

TGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGG  >  NC_014034/221972‑222181
                                                                                                             |                                                                                                    
          cggatccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtggggtccgtcg                                                                                                     >  2:4066988/1‑101 (MQ=14)
          cggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGGCGGGTCCGTCGGATCCgtggggtccgtcg                                                                                                     >  2:4005287/1‑101 (MQ=11)
             atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcggggccgtcggat                                                                                                  >  1:141785/2‑99 (MQ=14)
              tccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcggatc                                                                                                 <  1:2934253/101‑1 (MQ=11)
               tcgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCgtggggtccgtcgggtcc                                                                                                <  2:1732942/100‑1 (MQ=11)
               ccgtgggctcgggcggGTCCGTGGGATCGGTCGGGTCCGTGGTATCGGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtggggtcc                                                                                                <  1:1806957/101‑1 (MQ=11)
                 gtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTCCGTGGGATCGGTcgggtccgtggggtccgt                                                                                              >  2:253771/1‑101 (MQ=11)
                 gtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTGGGATCGGTCGCGTCCGTGGGATCCgtcgcgccggtcggcgcggt                                                                                              >  2:2846246/1‑95 (MQ=11)
                  tgggttcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcgggtccgtggggtccgtc                                                                                             <  2:3340822/101‑2 (MQ=11)
                  gcgggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTg                                                                                             <  1:3820588/93‑1 (MQ=11)
                   gggatgggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTgg                                                                                            <  2:2332464/101‑1 (MQ=11)
                   gggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTgg                                                                                            >  2:1191368/1‑101 (MQ=255)
                   gggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCTTGGGCTCGgtcggagccgtgggatgggccg                                                                                            >  2:3483165/1‑93 (MQ=11)
                   cgggtccgtcggATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTgg                                                                                            >  1:4144948/8‑101 (MQ=11)
                    ggttccgttggATCCGTGGTGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGTGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTggg                                                                                           <  2:363214/92‑1 (MQ=11)
                     gatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGGTCCGTGGGa                                                                                          <  1:1191462/101‑1 (MQ=255)
                      gtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAt                                                                                         >  1:645994/8‑101 (MQ=14)
                      atcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAt                                                                                         <  1:2476616/101‑1 (MQ=11)
                       tcggtgggATCGGTGGGATCGGTCGGATCCGTGGGATCGGTTGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtc                                                                                        <  2:1876789/101‑1 (MQ=11)
                       tcggtcggATCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtc                                                                                        <  1:1242535/101‑1 (MQ=11)
                       tcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCTGATCGGTCGGAtc                                                                                        <  2:395231/101‑1 (MQ=11)
                        gaggggggcccGTTGGAGAGGTCGGGTGCGTGGGATCGGTGGGGGCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtcc                                                                                       <  1:1072221/92‑2 (MQ=11)
                        ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg                                                                                       <  2:3060736/96‑1 (MQ=14)
                        ccgtcggATCTGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGAtcc                                                                                       <  2:1298320/99‑2 (MQ=11)
                        ccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGAtcc                                                                                       >  1:2144037/3‑100 (MQ=11)
                         ggtgtgctCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccg                                                                                      <  1:1010002/94‑3 (MQ=11)
                          gtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCAGTGGGGTGCGTCGGGTCCGTGGGGtcggt                                                                                     >  2:2988000/1‑101 (MQ=14)
                           tcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGGGCCGTCGGGTCGGTGGGAtcggtg                                                                                    >  2:382816/1‑100 (MQ=14)
                            ggggtcggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                                   <  2:2659704/91‑1 (MQ=17)
                            gggagcggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcg                                                                                   <  1:3036078/91‑1 (MQ=17)
                             gggttcGTTGGACCGGGTGGGGCCGTGGGGGTTGTCGGATCCGTGGGATGGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                                                                  <  1:405848/96‑1 (MQ=2)
                             gggtcggTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtcgg                                                                                  <  1:3734451/95‑1 (MQ=11)
                             ggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgg                                                                                  <  2:154198/101‑1 (MQ=17)
                             ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTGGGAtccgcggg                                                                                  >  1:1127882/1‑95 (MQ=14)
                              ggtCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtgggtcggg                                                                                 >  2:2229011/3‑100 (MQ=11)
                              gATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcggtcgag                                                                                 >  1:3969635/1‑99 (MQ=17)
                              gATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGa                                                                                 >  2:2758112/1‑101 (MQ=11)
                               gtCCTTGGGGTTGTTCCGGTCAGTGGGGTTCGTCCGGTCCGTGGGACCGGTCGGTTCCGTGGGATCCGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                                                <  2:1129180/100‑1 (MQ=9)
                               gtCCGTAGGACCGGCGGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcgggt                                                                                <  2:1597731/100‑3 (MQ=14)
                                gccGTGGGATCGGGCGGGTCCGTCGGATCGGGCGGATCCGTTGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTc                                                                               <  1:4005547/100‑1 (MQ=11)
                                 ccGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGAtcggtcggggcg                                                                              >  2:44523/1‑97 (MQ=14)
                                 ccGTAGGATCGGTCGGGTCCGTCGGAGCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                              <  1:3335024/101‑1 (MQ=14)
                                  ggtcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCg                                                                             >  2:3409463/5‑101 (MQ=17)
                                   gtaggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTGCGt                                                                            >  1:3272533/4‑101 (MQ=14)
                                   gTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGCTCGGTCGGGTCCTTCTGGTCCGTCGGAtccgtgggggcggt                                                                            >  2:1895616/1‑95 (MQ=14)
                                     gggctctctggggTCCGGCCGATCGGTGGGATCGGTAGGGTCCGTGGGATCGGGAGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                                                          <  2:787418/91‑1 (MQ=1)
                                     aggATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGACGGGTCCGTGGGGTCGGTGAGGGCCGTCGGAGCCGTGGGCGCGGTCGGTTCAGtcg                                                                          >  2:2371413/2‑99 (MQ=2)
                                          ccgtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGCCCGTGGGATCGGTCGGATCCGTGGGAtcg                                                                     >  2:2492135/6‑101 (MQ=255)
                                           ggtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                    >  1:595810/1‑101 (MQ=17)
                                           ggtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                                                    >  1:843627/1‑101 (MQ=14)
                                            gtggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGAtcggc                                                                   >  2:869231/4‑100 (MQ=255)
                                            gtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt                                                                   >  1:2312432/1‑101 (MQ=17)
                                            gtcggGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcggt                                                                   <  2:1091579/101‑1 (MQ=14)
                                            gtcggGGCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                                   <  1:1376482/101‑1 (MQ=17)
                                             tggggTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGGCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                  >  2:4160559/3‑101 (MQ=255)
                                               gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                >  2:3230250/1‑101 (MQ=17)
                                               gggTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcgg                                                                >  1:3886261/1‑101 (MQ=17)
                                                   gcGTGGGCTCGGTCGGGTCCGTGGGATCGGTCGGGTCGGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTcc                                                            <  2:1521221/100‑1 (MQ=255)
                                                   cggTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                            <  1:97338/99‑1 (MQ=17)
                                                   ccGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCGGTCGGGTcc                                                            >  2:2920537/1‑101 (MQ=255)
                                                   ccGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc                                                            <  2:570414/101‑1 (MQ=17)
                                                    cGTGGGTGCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                           <  2:3563930/101‑1 (MQ=255)
                                                     gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGt                                                          <  1:1878848/101‑1 (MQ=255)
                                                     gTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGCATCGGTCGGGTCCGt                                                          >  2:1791296/1‑101 (MQ=11)
                                                      tGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTg                                                         >  2:1236862/1‑101 (MQ=255)
                                                      tGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGta                                                         <  2:949484/101‑2 (MQ=17)
                                                       cggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgg                                                        <  1:2643110/100‑1 (MQ=17)
                                                        ggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggg                                                       >  1:2516790/1‑101 (MQ=255)
                                                          aTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggat                                                     <  1:1957878/101‑3 (MQ=14)
                                                           tcggtcggATCGGTGGGACCGGTCGGGTCCGTGGGGTCGGTCGGACCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatc                                                    <  2:512645/101‑4 (MQ=14)
                                                            ttcgcggAGCGGTGGGCTCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcg                                                   <  2:2055583/97‑5 (MQ=1)
                                                            cggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtg                                                   >  1:2053752/1‑101 (MQ=255)
                                                             ggtcggGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtgg                                                  >  1:3508511/1‑101 (MQ=255)
                                                             ggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcgg                                                  <  2:2896757/101‑6 (MQ=14)
                                                             ggtcggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcgg                                                  >  2:1679120/1‑96 (MQ=17)
                                                               tcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                                >  1:516569/1‑94 (MQ=14)
                                                               tcggATCCGGGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtc                                                <  1:2898056/101‑8 (MQ=11)
                                                               ggggATCCGTGGGGGCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTaggatcggtc                                                <  1:3666006/99‑8 (MQ=14)
                                                                tgggtttgTGTGATCTGTCGGGTCCGTAGGGTCGGTCGGTTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCtggtgg                                               <  2:3786050/94‑1 (MQ=18)
                                                                cggATCCGTGGGATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTggggtcggtcg                                               >  1:1680346/1‑94 (MQ=14)
                                                                   gcccGTGAGATGGGTGGGGTCCGGAGGATCCGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAt                                            <  2:1038449/99‑1 (MQ=21)
                                                                     ccGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGGGGGAtcg                                          >  1:3091657/1‑101 (MQ=255)
                                                                       gTGGGATCGGTCGGGTCCGTAAGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCAGTGGGGCTGGTGGGAtcggt                                        >  1:2398914/1‑101 (MQ=255)
                                                                       gTGGGATCGGGCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGAtcggt                                        <  1:1478024/101‑1 (MQ=255)
                                                                            tctggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGAc                                   <  2:2168234/98‑1 (MQ=255)
                                                                              cggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCg                                 >  2:1974625/1‑101 (MQ=255)
                                                                               ggtcggGTCCGTGGGATCTGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGc                                <  1:1541068/101‑1 (MQ=255)
                                                                                   gggTCGGTCGGTTCGGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTc                            <  2:365251/101‑1 (MQ=17)
                                                                                    ggTCCGGAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcg                           <  2:637886/101‑1 (MQ=255)
                                                                                       ccGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccg                        <  1:3409708/101‑1 (MQ=255)
                                                                                       ccGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccg                        >  1:748146/1‑101 (MQ=255)
                                                                                        cGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTcgccgc                       >  2:1415958/1‑101 (MQ=255)
                                                                                              aTCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCgg                 >  1:1621329/1‑101 (MQ=255)
                                                                                                 ggtcggGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACa              >  1:432231/1‑101 (MQ=255)
                                                                                                  gtggggtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAg             <  1:2853410/98‑1 (MQ=255)
                                                                                                    cgggtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAc           >  2:661518/1‑101 (MQ=255)
                                                                                                     gggtccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGAcc          >  2:3178192/1‑101 (MQ=255)
                                                                                                         ccgtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAgg      >  2:3861830/1‑101 (MQ=255)
                                                                                                           gtcgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAggcc    >  1:3562599/1‑101 (MQ=255)
                                                                                                             cgggtccgtGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAggccgg  >  2:3701871/1‑101 (MQ=255)
                                                                                                             |                                                                                                    
TGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCTGGTGGGATCGGTCGGACCGCTGTCGCCGCCGCCGGACAGACCCAGGCCGG  >  NC_014034/221972‑222181

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