Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,951,674	+G	95.5%	intergenic (‑81/‑111)	ssb1 ← / → RCAP_rcc01806	single‑stranded DNA‑binding protein/transglycosylase, Slt family

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,951,674	1	.→G	95.5%	2.4 0.0 ‑0.6	67	intergenic (‑81/‑111)	ssb1/RCAP_rcc01806	single‑stranded DNA‑binding protein/transglycosylase, Slt family
Reads supporting (aligned to +/- strand):  new base (32/32):  ref base (2/1):  total (34/33)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.68e-01

TTGTTCACCGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.GG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAAGAGGGC  >  NC_014034/1951574‑1951768                                                                                                      |                                                                                                 ttcttCGCCGCTCCTGCCATGAATCTCTCCGTAGTGTCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTc                                                                                                   <  1:1984642/98‑1 (MQ=255) ttGTTCCCCGACCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTc                                                                                                   <  1:4040023/101‑1 (MQ=255) ttGTTCACCGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTc                                                                                                   <  1:2851636/101‑1 (MQ=255) ttGTTCACCGACCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTc                                                                                                   <  1:1514522/101‑1 (MQ=255)  tGTTCCCCGATCCCGCCATGAAACCCTCCGTAGGGGCCGGAGGGACCCGGCGATTTTCCGAAATGCCCCCGGCCGGGGCCGCGGCAAGGCGGCCCCGGTC.g                                                                                                 <  2:461060/101‑1 (MQ=255)   tttCACCGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg                                                                                                <  2:1666584/100‑1 (MQ=255)   gTTCACCGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg                                                                                                >  1:1124002/1‑101 (MQ=255)     tctccGATCCTGCCATTGAACTATCCGTATTGGCCGGAGCGACCCGGCGATTGTCCGAACTGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gg                                                                                             <  2:780578/98‑1 (MQ=255)     tccccGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gg                                                                                             <  1:2779307/98‑1 (MQ=255)     tCACCGACCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gg                                                                                             <  2:4213293/101‑1 (MQ=255)       ccggacccTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gggc                                                                                           <  2:1102262/95‑3 (MQ=255)       cccGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gggc                                                                                           <  1:357697/100‑3 (MQ=255)        ccGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gggcg                                                                                          <  2:2205894/101‑4 (MQ=255)         cGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gggcgg                                                                                         >  2:2093789/1‑97 (MQ=255)          gATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.gg.gggcgtc                                                                                        >  1:2809372/1‑96 (MQ=255)                  cATGAACCTCTCCGTAGTGGCCGGAGCGCCCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATAc                                                                                <  2:3186818/101‑1 (MQ=255)                  cATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATAc                                                                                >  2:3533282/1‑101 (MQ=255)                      aaaCTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGGCGGCAATACTTTg                                                                            >  2:3169877/1‑101 (MQ=255)                      aaaCTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCCGTCGGG.GG.CGGCGGCaacagcttg                                                                            >  2:4543614/1‑94 (MQ=255)                       aaCTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGCCGGG.GG.CGCCGGCAATGCATTGc                                                                           >  1:17985/1‑101 (MQ=255)                        aCTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGcc                                                                          >  1:4308538/1‑101 (MQ=255)                         ctctCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGGCGGCAATACTTTGccc                                                                         >  2:1096952/1‑101 (MQ=255)                           ctcCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCAt                                                                       >  1:284124/1‑101 (MQ=255)                             cggTAGTGGCCGAAGAGACACGTCGATTGGCCGAACTCCCCCCTACCGGGGCCGCGGCACGGCGTCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATaa                                                                     <  2:776588/99‑1 (MQ=255)                              cGTGGTGGCCGGGGCTGCCCGGCGATTGTCCGAGTGGCCCCCGTCGGGGGCCGCGGCAAGGCGGCCGCGGTC.GG.GGCCGTCGGCAATACTCTGcacgcaaa                                                                    <  2:141585/101‑8 (MQ=255)                               gTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAAt                                                                   <  2:296667/101‑1 (MQ=255)                               gTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGGCGGG.GG.CGTCGGCAATACTTTGCCCATAAAt                                                                   >  2:1227823/1‑101 (MQ=255)                                tAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGGCGGG.GG.CGGCGGCAATACTTTGCCCATAAATc                                                                  >  1:3156688/1‑101 (MQ=255)                                   tGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAGGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTGGCCCATAAATAggg                                                               >  2:243222/1‑101 (MQ=255)                                   tGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCGATACTTTGCCCATAAATCggg                                                               >  2:1868422/1‑101 (MQ=255)                                      ccGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCcg                                                            <  1:2670804/101‑1 (MQ=255)                                      ccGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCcg                                                            >  2:335864/1‑101 (MQ=255)                                        ggAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCcgcg                                                          <  1:4421587/101‑1 (MQ=255)                                          aGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCgcggc                                                        >  1:4473884/1‑101 (MQ=255)                                            cGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGgcgc                                                      <  2:4002659/101‑1 (MQ=255)                                               cccGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAg                                                   >  2:628249/1‑101 (MQ=255)                                               cccGGCGATTGCCCGAAATGCCCCCGCCGGGGCCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAg                                                   <  1:2821478/101‑1 (MQ=255)                                                 cGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGAc                                                 <  1:4419349/101‑1 (MQ=255)                                                 cGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGAc                                                 <  2:359459/101‑1 (MQ=255)                                                 cGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGAc                                                 >  1:3035180/1‑101 (MQ=255)                                                 cGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGAc                                                 <  1:2093929/101‑1 (MQ=255)                                                   gCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGc                                               >  2:1673991/1‑101 (MQ=255)                                                   gCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGc                                               >  2:2760727/1‑101 (MQ=255)                                                     gATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTa                                             >  2:3391988/1‑101 (MQ=255)                                                     gATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTa                                             <  1:4618465/101‑1 (MQ=255)                                                      aTTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTag                                            <  2:1687061/101‑1 (MQ=255)                                                       ttGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.GGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCCAGACGCTagg                                           >  1:162396/1‑101 (MQ=255)                                                        tGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTaggc                                          >  1:1354032/1‑101 (MQ=255)                                                         gTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCTGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTaggca                                         <  2:3598987/101‑1 (MQ=255)                                                         gTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTaggca                                         >  2:742803/1‑101 (MQ=255)                                                           ccGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTaggcagg                                       <  1:3821952/101‑1 (MQ=255)                                                             gAAATGCCCCCGACCTGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCg                                     >  1:1042761/1‑101 (MQ=255)                                                              aaaTGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCgg                                    >  1:2560215/1‑101 (MQ=255)                                                              aaaTGCCCCCGACCGGGGCCGCGGCAAGGCCGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCgg                                    <  2:3862327/101‑1 (MQ=255)                                                               aaggCCCCGGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGa                                   <  2:980587/98‑1 (MQ=255)                                                                aTGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGAc                                  >  1:2835025/1‑101 (MQ=255)                                                                  gCCCCCGCCCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCg                                <  2:3266477/101‑1 (MQ=255)                                                                   cccccGCCCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGa                               <  1:4380877/101‑1 (MQ=255)                                                                    ccccGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGAGGCGCAAGACGCTAGGCAGGCGGACCGAt                              >  2:1131802/1‑101 (MQ=255)                                                                     cccGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGAtt                             >  2:2375540/1‑101 (MQ=255)                                                                     cccGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGAtt                             <  1:3533540/101‑1 (MQ=255)                                                                      cggccggggcCCGGGGCAGGGCGCCCCGGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTg                            <  1:3197951/95‑1 (MQ=255)                                                                      ccGACCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTg                            <  2:2563155/101‑1 (MQ=255)                                                                         aCCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTgcgc                         >  2:3810000/1‑101 (MQ=255)                                                                         aCCGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAAACGGGCCGCGGCGCAAGACGCGAGGCAGGCGGACCGAGTgcgc                         >  2:1834836/1‑101 (MQ=255)                                                                           cGGGGCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCt                       >  1:260585/1‑101 (MQ=255)                                                                            ggggCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCta                      >  1:1564150/1‑101 (MQ=255)                                                                            ggggCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCta                      >  1:1699401/1‑101 (MQ=255)                                                                            ggggCCGCGGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCta                      >  2:4272986/1‑101 (MQ=255)                                                                                 cgcgGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGtt                 <  1:3392241/101‑1 (MQ=255)                                                                                  gcgGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTc                >  1:2740621/1‑101 (MQ=255)                                                                                  gcgGCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTc                <  2:4642210/101‑1 (MQ=255)                                                                                    ggCAAGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCg              <  1:793261/101‑1 (MQ=255)                                                                                       aaGGCGGCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGcc           <  1:2962200/101‑1 (MQ=255)                                                                                       aaGGCGCCCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGcc           <  2:3156444/101‑1 (MQ=255)                                                                                        aGGCGGCCCCGGTC.GGTGG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCg          >  1:2096971/1‑101 (MQ=255)                                                                                         ggcgccCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGa         <  2:17980/101‑1 (MQ=255)                                                                                          gcggcCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGaa        <  2:888548/101‑1 (MQ=255)                                                                                           cggcCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAag       <  1:1868548/101‑1 (MQ=255)                                                                                            ggcCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAaga      >  1:3016026/1‑101 (MQ=255)                                                                                             gcCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAagag     <  2:2809156/101‑1 (MQ=255)                                                                                             gcCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAagag     <  2:1172455/101‑1 (MQ=255)                                                                                             gcCCCGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAagag     >  1:1701209/1‑101 (MQ=255)                                                                                             gcCCCGGTC.GGAGG.CGTCGGCAATACTTTGCGCATAAATCGGGCGGCGGCGCAAGACGCTAGGCAGGGGGCCCGATTTCGCCTATCGGTCCCGCCGGagag     >  2:4576111/1‑101 (MQ=255)                                                                                               cccGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAAGAggg   >  2:3597112/1‑101 (MQ=255)                                                                                                ccGGTCGGG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAAGAGGGc  <  1:3547628/101‑1 (MQ=255)                                                                                                      |                                                                                                 TTGTTCACCGATCCTGCCATGAAACTCTCCGTAGTGGCCGGAGCGACCCGGCGATTGTCCGAAATGCCCCCGACCGGGGCCGCGGCAAGGCGGCCCCGGTC.GG.GG.CGTCGGCAATACTTTGCCCATAAATCGGGCCGCGGCGCAAGACGCTAGGCAGGCGGACCGATTGCGCCTATAGTTCCGGCCGAAGAGGGC  >  NC_014034/1951574‑1951768 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG