Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 222,009 G→A 8.0% P216P (CCC→CCT RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034222,0090G→A 8.0% 0.0 ‑0.8 ‑0.02577P216P (CCC→CCTRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (95/102):  ref base (1142/1096):  total (1321/1256)
Fisher's exact test strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.47e-01

Only 100 of 2983 total aligned reads displayed.

CCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGG  >  NC_014034/221909‑222109
                                                                                                    |                                                                                                    
   gATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGa                                                                                                   >  1:4218367/1‑101 (MQ=255)
          gggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGggatcgggcg                                                                                            <  1:230318‑M1/64‑4 (MQ=255)
          cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                                                            >  1:4596071/1‑101 (MQ=255)
              tcggtcgggtccgtgggntcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                                        >  2:1312535‑M1/34‑101 (MQ=255)
              tCGACGGTGGGGGTCGGATGCGTGGGCTCGGTCGGGTCGGTCGGGTCCGGGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                                                                        <  2:4701319/101‑1 (MQ=255)
               cggtcgggtccgtgggagcggtcggggccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                                                       >  2:4348315‑M1/33‑101 (MQ=255)
               cggtcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATcc                                                                                       >  2:4874251‑M1/33‑101 (MQ=255)
                ggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCg                                                                                      >  1:2564162‑M1/32‑101 (MQ=255)
                 gtgggatcggTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGt                                                                                     <  1:2765005/92‑1 (MQ=11)
                     gatccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGa                                                                                 >  1:5413867‑M1/27‑101 (MQ=255)
                      gtgcgttggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                                                >  1:1495102‑M1/26‑101 (MQ=255)
                       tccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                               >  2:3891838‑M1/25‑101 (MQ=255)
                       gccgtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGAtc                                                                               >  2:2017713‑M1/25‑101 (MQ=255)
                           tgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                                           >  2:5797468‑M1/21‑101 (MQ=255)
                              gatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcggg                                                                        >  1:3902129‑M1/18‑100 (MQ=255)
                                tGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATc                                                                      <  1:6650191/101‑1 (MQ=255)
                                 ccGTCGGGTCGGTCGCAGCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATcg                                                                     <  2:1022248/100‑2 (MQ=14)
                                   gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGt                                                                   >  1:4890891/1‑101 (MQ=255)
                                   gTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGt                                                                   <  1:2669116/101‑1 (MQ=255)
                                     gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                                                                 >  1:2847987/1‑99 (MQ=1)
                                     cggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTgg                                                                 >  2:1879488‑M1/11‑101 (MQ=255)
                                       ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGa                                                               >  1:891737‑M1/9‑101 (MQ=255)
                                        gtgcgtgggatccgtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAt                                                              <  2:6869419/94‑1 (MQ=11)
                                        gtccgtgGGATCGGTCGGCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                                                              >  1:7025990‑M1/8‑101 (MQ=255)
                                        ggcggtcggCGCCGTCGCATCGGTGTGATCCGTCGGGTTCGTGGGATCTGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTcgggt                                                              <  2:2310696/99‑3 (MQ=14)
                                        atccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGATCGGTGGGAt                                                              <  1:3995647‑M1/94‑1 (MQ=255)
                                          ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtcg                                                            <  1:4475328‑M1/96‑1 (MQ=255)
                                           cgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtcgg                                                           >  1:6384297/2‑101 (MQ=14)
                                           cgtcggAGCGGTCGGGCCCGTGGGATCGTTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccg                                                           <  1:1781737/100‑3 (MQ=11)
                                             tgggatcggtcgggtccgtggggtccgtcggatccgtgGGATAGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtg                                                         >  2:3440685/3‑100 (MQ=11)
                                              gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg                                                        >  2:3625454‑M1/2‑101 (MQ=255)
                                              gGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                                        <  1:6538979‑M1/100‑1 (MQ=255)
                                               ggttcggtcggatccggcggatccgtcgggtccgtcggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcgg                                                       <  2:3535218/98‑1 (MQ=14)
                                                gatcggtcgggtccgtgggatcggtcggatcggtcgGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCggg                                                      <  1:734325/101‑1 (MQ=11)
                                                 attggtcggatccgttggatcggtcggatccgtaggatcGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCTGTCGGGt                                                     <  2:7043369/98‑1 (MQ=11)
                                                  tcggtcgggtccgtgggatcggtcgggtccgtcggaTCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTc                                                    <  1:4722957/101‑1 (MQ=11)
                                                  tcggtcggatccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATc                                                    <  2:4136769/101‑1 (MQ=11)
                                                   cggtcgggcctgtgggatccgtcgggtccgtgggatCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                                                   <  2:3051311/101‑1 (MQ=14)
                                                    ggtcgggtccgtgggatcggtcgggtccgtcggatcGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCg                                                  <  1:2155384/101‑1 (MQ=11)
                                                     gtcgggtccgtgggctcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGt                                                 <  2:6652116/101‑1 (MQ=14)
                                                     gtcgggtccgtgggatcggtcgggtacgtgggatcgGTCGGGTGCAGGGGATGGGTGGGGTTCTTGGGGTTCGCCGGATCCTTGGGATAGGTCGGGTCAGt                                                 >  2:2776181/1‑101 (MQ=9)
                                                      tcgcgtctgtgggatggggcgtatccgtgggatcggTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGtg                                                <  2:5137244/101‑2 (MQ=11)
                                                       ggggtccgtgggatccgtcgggtccgtgggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                               >  1:1880412‑M2/2‑99 (MQ=255)
                                                       cggatccgtgggatcggtcgggtccgtgggatcggtCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgg                                               <  1:324346‑M2/101‑3 (MQ=255)
                                                        gggtccgtgggatcggtcgggtccgtggggtccgtcGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGACCGGCCGGGTCCGTGgg                                              >  1:5955755/1‑101 (MQ=14)
                                                        gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGGTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTTCGTCGGGTccttggg                                              >  2:675069/1‑96 (MQ=11)
                                                        gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggg                                              >  2:2676385‑M2/1‑98 (MQ=255)
                                                        gggtccgtgggatcggtcgggtccgtgggatcggtcGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTggg                                              <  1:2959230‑M2/101‑4 (MQ=255)
                                                         ggtccgtgggatcggtcgggtccgtggggtccgtcgGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                             >  2:4251529‑M2/1‑97 (MQ=255)
                                                         ggtccgtgggatcggtcgggtccgtgggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                             <  2:177006‑M2/101‑5 (MQ=255)
                                                         ggtccgggggatcggtcgggtccgtgggatcggtcgGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTggga                                             >  1:6232247‑M2/1‑97 (MQ=255)
                                                          atccgtgggatcggtcggatccgtaggatcggtcggGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTCGGAt                                            <  2:4459887/100‑1 (MQ=11)
                                                           tccgtgggatcggtcggctccgtgggatcggtcggCTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatc                                           >  2:2431483‑M2/1‑95 (MQ=255)
                                                            ccgtgggctccgtgcggtccgtggggtccgtcggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcg                                          <  2:2120971/101‑1 (MQ=1)
                                                            ccgtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTggggtcc                                          >  1:6084984‑M2/1‑94 (MQ=255)
                                                            ccgtgggatcggtcgggtccgggggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcg                                          >  1:122434‑M2/1‑94 (MQ=255)
                                                            ccgtgggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg                                          >  1:5315881/1‑101 (MQ=11)
                                                             cgtgggatcggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGgggcggg                                         >  1:1000254/1‑96 (MQ=11)
                                                              gtgggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTCGGATCGGTGGGGtcggg                                        >  2:6063537/1‑100 (MQ=11)
                                                               tgggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtc                                       >  2:4988183‑M2/1‑91 (MQ=255)
                                                               tgggatcggtcgggtccgggggatcggtcggGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                                       <  1:2020875/101‑1 (MQ=14)
                                                                gggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                                      <  2:3159326‑M2/101‑12 (MQ=255)
                                                                gggatcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcg                                      <  2:891713‑M2/101‑12 (MQ=255)
                                                                gggatcggtcggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                      >  2:731374/1‑101 (MQ=14)
                                                                 ggatcggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgg                                     <  2:6200235‑M2/101‑13 (MQ=255)
                                                                   atcggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcggac                                   <  1:3471120‑M2/101‑15 (MQ=255)
                                                                    tgcgtcgggtccgtcgggtcggtgggATCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGgtcgggtc                                  <  2:5914519/98‑1 (MQ=11)
                                                                     cggtcgggtccgtgggatcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCGGTgggatcggtcggatcc                                 <  2:4561865‑M2/101‑17 (MQ=255)
                                                                      ggtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGTCGgtgggatccg                                >  1:4133468/1‑101 (MQ=14)
                                                                      ggtcgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGCCCGTGGCGTCGgtcgggtccg                                >  1:1097804/1‑101 (MQ=11)
                                                                      gggtgggtccgtgtgatccctaggATCCTTCGTATCGGCCCGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccg                                <  1:5694048/97‑1 (MQ=11)
                                                                      cctctgagctgtgggacccggcggAGCCGTGGGATCGGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGgtcgggtccg                                <  2:2885899/91‑1 (MQ=14)
                                                                       gtcgggtccgtggggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgt                               >  1:6768093‑M2/1‑83 (MQ=255)
                                                                         cgggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgg                             <  2:5521058‑M2/101‑21 (MQ=255)
                                                                         cgggtccgtcggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGGGGGATCGGTCGGGTCCGTGGGATCGgtcgggtccgtcg                             >  2:2569088/1‑101 (MQ=11)
                                                                          ggatccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggg                            <  1:1642353‑M2/101‑22 (MQ=255)
                                                                          ggatccgtgggatcggtcggATCCGTGGGATCGATCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcgggtccgtcgg                            >  1:5017954/4‑101 (MQ=14)
                                                                           ggtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgggggg                           <  2:3775451‑M2/101‑23 (MQ=255)
                                                                            gtccgtgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGACCGGTCGGACCCGTGGGATGGGTCGGGTCCGTGGGGCCGgtcggatcccggggcg                          >  2:1459004/1‑94 (MQ=11)
                                                                            atccgtgggatcggtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGgtcggggccgtcggat                          >  2:6338151/2‑99 (MQ=14)
                                                                                gtgggatcggtcggATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGATCCGTGGGATCGGtcgggtccgtcggatccgg                      >  1:370106/1‑100 (MQ=11)
                                                                                 tgggatcggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCAGTGGGGTGGgtggggtacgctgggtgggca                     >  2:5407913/1‑96 (MQ=11)
                                                                                   gggtccgtcggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTcgg                   >  2:3989943/7‑101 (MQ=11)
                                                                                    ggtccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGa                  >  1:5176349/9‑101 (MQ=11)
                                                                                    gggtggtccgACCCGTGCGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTcggg                  <  2:3270091/97‑2 (MQ=11)
                                                                                     atcggtcggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcggatccgtgggatcggtcggat                 <  2:56704‑M2/101‑33 (MQ=255)
                                                                                       cggtgggATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc               >  1:4348489‑M2/1‑67 (MQ=255)
                                                                                       cggtcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtcc               >  2:5276867‑M2/1‑67 (MQ=255)
                                                                                       ccgtgggATCGGTCGGATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGAtcg               >  1:5816942/6‑101 (MQ=14)
                                                                                         gtccgATCCGTTGTTTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGAtccgt             <  2:1613359/101‑1 (MQ=11)
                                                                                          tcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTAGGacccgtg            >  1:333670/1‑95 (MQ=11)
                                                                                          gcggGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGTCCGTGGGAtcggtc            <  1:3096678/100‑1 (MQ=11)
                                                                                           cggGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgg           >  2:1157863‑M2/1‑63 (MQ=255)
                                                                                            ggATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTAGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGATCGGTCGGAtccgtggg          <  1:381422/101‑1 (MQ=11)
                                                                                             ggtccGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtggggtcggtcgggtccgtggga         <  1:3574329‑M2/99‑41 (MQ=255)
                                                                                                ccGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg      >  1:243242‑M2/1‑58 (MQ=255)
                                                                                                  gTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggt    <  1:555917‑M2/101‑46 (MQ=255)
                                                                                                   tGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtg   <  2:6515980‑M2/101‑47 (MQ=255)
                                                                                                    gggATGGGTCGGGTCGGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg  <  1:6869773‑M2/101‑48 (MQ=255)
                                                                                                    gggATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcg  <  1:1333710‑M2/101‑48 (MQ=255)
                                                                                                    |                                                                                                    
CCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGG  >  NC_014034/221909‑222109

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