Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 221,964 C→G 2.5% P231P (CCG→CCC RCAP_rcc00185 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034221,9640C→G 2.5% 15.3 ‑0.5 ‑0.2525P231P (CCG→CCCRCAP_rcc00185hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (239/271):  total (244/281)
Fisher's exact test strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.56e-01

Only 100 of 1265 total aligned reads displayed.

TGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAG  >  NC_014034/221864‑222064
                                                                                                    |                                                                                                    
tgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                                                      <  1:2166776‑M1/9‑1 (MQ=255)
tgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtc                                                                                                      >  2:3352654‑M1/93‑101 (MQ=255)
  ggatcggtcgggtccgtggggtcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgg                                                                                                    >  2:1369599‑M1/91‑101 (MQ=255)
   gatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggg                                                                                                   <  1:3385186‑M1/12‑1 (MQ=255)
         tcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtg                                                                                             <  2:1148960‑M1/18‑1 (MQ=255)
         ccGTCGCGGAGGTGTCGGTGGTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGGCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                                             <  2:881120/101‑1 (MQ=255)
         ccGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtg                                                                                             <  2:2954613/101‑1 (MQ=255)
          cgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                                            <  2:2408581‑M1/19‑1 (MQ=255)
              tccgtggggtcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatc                                                                                        >  1:2879443‑M1/79‑101 (MQ=255)
               ccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcg                                                                                       <  1:1124720‑M1/24‑1 (MQ=255)
                cGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcgg                                                                                      >  2:3874495/1‑101 (MQ=255)
                 gtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggt                                                                                     <  2:4324698‑M1/26‑1 (MQ=255)
                  aGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtc                                                                                    >  2:2633733/1‑101 (MQ=255)
                   ggCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcg                                                                                   <  2:1576670/101‑1 (MQ=255)
                      atcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccgtcggat                                                                                >  2:325018‑M1/71‑99 (MQ=255)
                        cGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCggatcggtcgggtccgtgggatcggtcgggtcc                                                                              <  2:1083251/101‑1 (MQ=255)
                         ggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcggatccg                                                                             <  2:1299155‑M1/34‑1 (MQ=255)
                            gtTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGgatcggtcgggtccgtgggatcggtcgggtccgtgg                                                                          >  1:2435846/1‑101 (MQ=255)
                            cggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcggatccgtgg                                                                          <  1:3659380‑M1/37‑1 (MQ=255)
                               atccgtgggatcggtcggatccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggat                                                                       >  1:4105033‑M1/62‑101 (MQ=255)
                                 ccgtcggatcggtgggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtcggatcggtcgggtccgtgggatcggtcgggtccgtgggatcg                                                                     >  1:1806083‑M1/60‑101 (MQ=255)
                                  cgtgggatcggtcgggtccgtgggatcggtcgggtccgggggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtggggtccg                                                                    >  1:1473353‑M1/59‑96 (MQ=255)
                                    tgggatcggtcgggtccgtgggatcggtcggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgggatcggtc                                                                  <  1:4286955‑M1/45‑1 (MQ=255)
                                             tcggatccgtgggatcggtcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTg                                                         <  2:576058‑M1/54‑1 (MQ=255)
                                              cggatccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTGGGATCGGTCGGAGCCGTgg                                                        <  2:4047660‑M1/55‑1 (MQ=255)
                                                    gcgcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                                  <  2:2550190/101‑1 (MQ=255)
                                                     gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                                                 >  1:231954‑M1/40‑101 (MQ=255)
                                                      tgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtc                                                >  1:1062679‑M1/39‑101 (MQ=255)
                                                      gcTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtc                                                >  2:1788020/1‑101 (MQ=255)
                                                       cTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                                               >  2:3080677/1‑101 (MQ=255)
                                                         gatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGgtccgtcgga                                             >  2:2322589‑M1/36‑100 (MQ=255)
                                                          atcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGt                                            >  2:3657695‑M1/35‑101 (MQ=255)
                                                           tcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGTGATCGGTCGGATCCGTGGGATCGGTCGGGTc                                           >  1:3981852‑M1/34‑101 (MQ=255)
                                                           tcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGggatcggtcggatc                                           >  1:3908648‑M1/34‑98 (MQ=255)
                                                               tcggatccgtgggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCGGGTccggg                                       <  2:1239703‑M1/72‑3 (MQ=255)
                                                               gcgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTg                                       >  1:1577101‑M1/30‑101 (MQ=255)
                                                                cgggtccgtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGCCGTgg                                      <  2:3160543‑M1/73‑1 (MQ=255)
                                                                cgggtcagtgggatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTgg                                      >  2:4361839‑M1/29‑101 (MQ=255)
                                                                  gatccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGa                                    >  1:562777‑M1/27‑101 (MQ=255)
                                                                      cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcgg                                <  1:4434320‑M1/79‑1 (MQ=255)
                                                                      cgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg                                <  2:945697‑M1/79‑1 (MQ=255)
                                                                       gtgggatcggtcggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                               <  1:12000‑M1/80‑1 (MQ=255)
                                                                       gTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggt                               <  2:2047280/101‑1 (MQ=255)
                                                                        ttcgttccGTTGGGTCCGTCGGAGCGGTCGGGTCCGTGGGATCGGTCGGGTCCTTTGCATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTCGGAtccgtg                              <  2:215016/94‑2 (MQ=18)
                                                                        tgggatgggtcggtgccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc                              >  1:3828647‑M1/21‑101 (MQ=255)
                                                                         ggggtccGTGGGGTCGGTGGGATCGGT.GGGTCCCTGGGGTTCGGCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGtccgtcgg                            <  1:950923/95‑1 (MQ=2)
                                                                         gggatcggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtcg                             <  1:1730168‑M1/82‑1 (MQ=255)
                                                                           gatcggtcgggtccgtcGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGatcggtcggg                           <  1:1355440‑M1/84‑2 (MQ=255)
                                                                            atcggtgggatccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGgatcggtcgggt                          <  2:3938319‑M1/85‑3 (MQ=255)
                                                                              gggtcgggtccgtgGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc                        >  1:4045859‑M1/15‑101 (MQ=255)
                                                                              cggtcgggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATcc                        <  1:3546489‑M1/87‑1 (MQ=255)
                                                                              ccgtcgtatcgggcggATCCGTGGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATcg                        <  2:732361/93‑2 (MQ=14)
                                                                                 tGGGCTCGGTCGGATCGGACGGGCCCGTGGGATCGGACGGGTCCGTGGGATCGGTCGGATCCGTGGGATCCGTGGGGTCCGTGGGATCGGTCGGATCCggg                     >  2:246925/1‑99 (MQ=17)
                                                                                  gggATCGGTCGGGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtcg                    >  1:4362110/1‑99 (MQ=1)
                                                                                   gggttcgtgggatgggtcgggtccgtgggatcggtcgggtccgtgGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTggg                   <  2:818666/92‑1 (MQ=11)
                                                                                   gggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                   <  1:787717‑M1/92‑1 (MQ=255)
                                                                                   ggatccgtgggatcggtcgggtccgtgggatcggtggggtcggtgGGAACGGCGGGGTGCGTGGGATCGGTCGGGGCCGTGGGATCGGTGCGATCTGTggg                   >  1:2103484/10‑101 (MQ=0)
                                                                                   ggCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTggg                   <  2:3822963/101‑1 (MQ=255)
                                                                                   ggCTCGGTCGGAGCCGTTGCGTCCGTGGGATCGGTCGCGTCCGTGGGATCGGTCGGGTCCGTAGGATCGGTCGGGTCCGTCGGGTCCGTGGGATCGGTcgg                   <  2:2844672/101‑1 (MQ=17)
                                                                                    ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                  <  1:88926‑M1/93‑1 (MQ=255)
                                                                                    ggtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGa                  >  1:2637212‑M1/9‑101 (MQ=255)
                                                                                     gtccgtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGGTGCGTGGGGTCCGTCGGGTCCGTGGGATCGGGCGGATCGGTGGGAt                 >  2:1957773/8‑101 (MQ=11)
                                                                                     gtccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCGGTCGGGTCCGTGGgag                 <  1:501404‑M1/94‑2 (MQ=255)
                                                                                     gtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAt                 <  1:3253278‑M1/94‑1 (MQ=255)
                                                                                      tcggtgggGTCCGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtc                <  1:3737507/101‑1 (MQ=14)
                                                                                      tccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGAtc                <  2:1860029‑M1/95‑1 (MQ=255)
                                                                                       ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcg               <  2:3285694‑M1/96‑1 (MQ=255)
                                                                                       ccgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGGCGGATCCGTGGGAtcg               <  1:2955516‑M1/96‑1 (MQ=255)
                                                                                       ccgtcggATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGCGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGAtag               >  2:2742604/3‑99 (MQ=11)
                                                                                        ggtcggGTCCGTCGGATCCGTGGGGTCCGTCGGGTCCGTGGGATCGGTCGGGGCCGGCGGGTCGGCAGGGCTCGTGGGATCGGACGGGTCCGTCGGGtcgg              >  1:2387418/1‑101 (MQ=2)
                                                                                        cgtggtatcggtcgggtccgtgggatcggtcgagtccgtgggGGCCGTCGGATGCGTTGGATCGGTCGGAGCCCTGGGATCAGTCCGATCCGCGGGAtcgt              >  1:307578/7‑100 (MQ=2)
                                                                                        cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg              >  1:131580‑M1/5‑101 (MQ=255)
                                                                                        cgtgGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcgg              <  1:2278998‑M1/97‑1 (MQ=255)
                                                                                         gtgggatcggtcgggtccgtgggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGGCCGGTCGGCTCCGGGGGGtccgg             >  1:670636/4‑98 (MQ=11)
                                                                                          tggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtc            <  2:3554516/93‑1 (MQ=14)
                                                                                          tgGGATCGGTCGGGTCCGTGGGAACGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtc            >  1:1912451‑M1/3‑101 (MQ=255)
                                                                                          tgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGAtcggtc            <  1:2721560‑M1/99‑1 (MQ=255)
                                                                                           ggggtccgtcgggtccgtgggatcggtcgggtccgtcggatcgGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGAtcggtcg           <  2:1656983/94‑1 (MQ=14)
                                                                                           gggatgggtcgggtccgggggatcggtcgggtccgtgGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtccgtgg           <  1:2015677/100‑3 (MQ=11)
                                                                                           ggaatcggtcgggtccgtggggtccgtcggatccgtgggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg           >  1:3307112/4‑101 (MQ=11)
                                                                                           gGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGAtcggtcg           >  2:646291‑M1/2‑101 (MQ=255)
                                                                                             gatcggtcgggtcggtgggggcggtcgggtccgtcgGAGCGGGGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTGGGGTCCGTCGGATCGGTCggg         <  1:3459921/101‑1 (MQ=11)
                                                                                             gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGGGCGGGGGATCGGTCggg         >  2:3458719/1‑101 (MQ=11)
                                                                                             gatcggtcgggtccgtgggatcggtcgggtccgtggGATCGGTCGGATCCGTGGGATCCGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtccgtggga         <  2:2271056/101‑2 (MQ=11)
                                                                                             gatcggtcgggtccgtgggatcggtcgggcccgtggGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCggg         >  1:3147835/1‑101 (MQ=14)
                                                                                              gtcggtgcggtgcgtgggatcggtggggtccgtgggaTCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGt        <  1:4181055/100‑1 (MQ=11)
                                                                                              gtcggtcgggtccgtgggatcggtcgggtccgtgggaTCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGGtccgtgggat        <  2:3988733/100‑3 (MQ=11)
                                                                                              atcggtcgggtccgtggggtccgtcggatccgtgggATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGGCgggg        >  2:1512829/1‑100 (MQ=11)
                                                                                               tcggtggggcccgtgggatcggtcgggtccgtgggaTCGGTCGGGTCCGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTc       <  2:2696931/101‑1 (MQ=11)
                                                                                                cggtcgggtccgtgggatcggtcgggtccgtcggatCGGTGGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTCGGGTcc      >  1:1657904/1‑101 (MQ=11)
                                                                                                cggtcgggtccgtgggacccgtcgggtccctcggacCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTcc      <  2:4103914/101‑1 (MQ=11)
                                                                                                cggtcggatccgtgggatcggtcgggtccgtgggatCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGATcc      <  2:454772/101‑1 (MQ=11)
                                                                                                cggtccggtccgtgggatcgttcgggtccgtgggatCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCCGATCCGTGGGATCGGTCGGATcc      <  2:1441491/101‑1 (MQ=14)
                                                                                                   tcgggtgcgtgggatcggtcgggtccgtgcgatcgcCCGGGTCCGTGGTGTCCCTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGtg   <  2:3721083/101‑2 (MQ=11)
                                                                                                   tcgggtccgtgggatcggtcgggtccgtggggtccgTCGGATCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtg   <  1:255270/101‑2 (MQ=11)
                                                                                                   tcgggtccgtgggatcggtcgggtccgtgggatcggTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTCGGATCCGTGGGATGGGTCGGATCCGtg   >  1:2772519/1‑100 (MQ=11)
                                                                                                    cgggtccgtgggatcggtcgggtccgtgggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGATCGGTGGGATCCGtgg  >  1:3202448/1‑99 (MQ=14)
                                                                                                    cgggtccgtgggatcggtcgggtccgtcggatcggtCGGGTCCGTGGGATCGGTCGGGTCCGTCGGATCGGTGGGATCCGTGGGGTGGGGCGGGTTCGtgg  >  2:1016294/1‑99 (MQ=14)
                                                                                                    cgggaccgtgggatgggtcgggtccgtgaggtcggtCCGATACGTGAGATCCGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGtgg  <  2:4243589/101‑3 (MQ=11)
                                                                                                    cggatccgtgggatcggtcggatccgtgggatcggtCGGGTCCGTGGGGTCGGTCGGGTCCGTGGGGTCCGTCGGATCCGTGGGGTCGGGGGGATCCGtgg  >  1:1004321/1‑99 (MQ=11)
                                                                                                    |                                                                                                    
TGCCCGACACCGTCGCCGAGGCTTCGGTGTTGGTGGTCAGCGTCGCCGGATCGCGCTCATCGACGGTGGGCGTCGGATGCGTGGGCTCGGTCGGATCGGTCGGGTCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTGGGATCGGTCGGATCCGTGGGATCGGTCGGGTCCGTAG  >  NC_014034/221864‑222064

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