Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,843,369	+C	100%	intergenic (+31/+117)	ppiB → / ← tyrS	peptidyl‑prolyl cis‑trans isomerase B/tyrosyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,843,369	1	.→C	100.0%	288.7	89	intergenic (+31/+117)	ppiB/tyrS	peptidyl‑prolyl cis‑trans isomerase B/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  new base (17/71):  ref base (1/0):  total (18/71)

GCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCA.GG.G.G.GCTCTGCCCCCCCTCT.G.TGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCG.TCCGCT.GGGGCGCCCTTTCATCTTGCCACAAATATCCCGGG  >  NC_014034/1843276‑1843461                                                                                                     |                                                                                               gcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCa                                                                                          <  2:3436897/101‑1 (MQ=255) gcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCa                                                                                          <  1:436772/101‑1 (MQ=255) gcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCa                                                                                          <  1:648945/101‑1 (MQ=255)  cgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCAt                                                                                         <  1:3637528/101‑1 (MQ=255)  cgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCAt                                                                                         <  1:1550226/101‑1 (MQ=255)  cgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCAt                                                                                         <  1:1562843/101‑1 (MQ=255)   gcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCAtt                                                                                        <  1:4129530/101‑1 (MQ=255)   gcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCAtt                                                                                        <  2:1586335/101‑1 (MQ=255)   gcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCAtt                                                                                        <  1:1208755/101‑1 (MQ=255)    cccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTc                                                                                       <  1:235927/101‑1 (MQ=255)    cccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTc                                                                                       <  1:4210029/101‑1 (MQ=255)    cccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTc                                                                                       <  2:306737/101‑1 (MQ=255)     cttccgatctcGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCa                                                                                      <  1:2727848/91‑1 (MQ=255)       gAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCAcc                                                                                    <  1:646576/101‑1 (MQ=255)        aGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCAccc                                                                                   >  1:4432211/1‑101 (MQ=255)         gCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCAcccc                                                                                  <  1:594346/101‑1 (MQ=255)          cTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCAccccc                                                                                 <  1:68179/101‑1 (MQ=255)           tGAACGGCCTTGTCGAAGCCCCCGATTACATCGGCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCAcccccc                                                                                <  1:1943639/101‑1 (MQ=255)           tGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCAcccccc                                                                                <  2:768171/101‑1 (MQ=255)            gAACGGCCTTGTCGAAGCCCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCg                                                                               <  1:4519095/101‑1 (MQ=255)             aaCGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACTCCCCgg                                                                              <  2:269920/101‑1 (MQ=255)             aaCGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCgg                                                                              <  1:840716/101‑1 (MQ=255)             aaCGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCgg                                                                              <  1:4458030/101‑1 (MQ=255)              aCGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCggg                                                                             <  2:3767636/101‑1 (MQ=38)               cGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCaagatcggaag                                                                            >  2:2727650/1‑91 (MQ=255)                ggCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTGCCCCCCC.CG.C.CGCGCCTTCAccccccccggta                                                                          >  1:3063159/1‑95 (MQ=255)                ggCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGat                                                                           <  2:4583199/101‑1 (MQ=38)                ggCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCACTCT.GCTGCGCATTCACCCCCCGGGat                                                                           <  1:951270/101‑1 (MQ=38)                ggCATTGTAGAAGCCCACAATTCCCTCGTCCCGGCCGCGATCGAGTGACAC.GGGCGCA.GG.GCG.GCTCAG.CCCCCCTCT.GCTGCTCATTCACCCCCCGGGat                                                                           <  2:2847708/101‑1 (MQ=255)                   cTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATAtt                                                                        <  1:1385388/101‑1 (MQ=37)                        cgagggccccGATTACAACGTCAAGGCCACCAACGAGTGTAACTGGGGGAG.GG.G.G.GCTCTGCCCCCCC.CT.GCTGCGCATTCACCCCCCGGGATATTTTCGc                                                                   <  2:3847525/95‑1 (MQ=12)                          aaGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCa                                                                 >  2:1479056/1‑101 (MQ=37)                          aaGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCACTCACCCCCCGCGGTCTTTTCTCCa                                                                 >  2:3665715/1‑101 (MQ=255)                           aGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCaa                                                                <  2:3397469/101‑1 (MQ=37)                           aGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCaa                                                                <  1:3446771/101‑1 (MQ=37)                                 cGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaa                                                          <  2:3448336/101‑1 (MQ=25)                                 cGATTACATCGTCAAGGCCACGAGCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaa                                                          <  1:649610/101‑1 (MQ=25)                                  gATTACCTCGTCAAGGCCACGATCGAGTGACTC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGaaa                                                         <  1:1509661/101‑1 (MQ=21)                                  gATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaaa                                                         <  2:3823821/101‑1 (MQ=25)                                  gATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCA.GGAG.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaaa                                                         >  2:1515379/1‑101 (MQ=25)                                   aTTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAAt                                                        <  1:3362858/101‑1 (MQ=25)                                    ttACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATg                                                       <  2:366617/101‑1 (MQ=25)                                     tACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATgg                                                      <  2:79434/101‑1 (MQ=25)                                     tACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATgg                                                      <  1:423565/101‑1 (MQ=25)                                       cATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaa                                                    <  2:4061575/101‑1 (MQ=25)                                        aTCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaa                                                   <  2:3244147/101‑1 (MQ=25)                                         tCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaaa                                                  <  1:3297256/101‑1 (MQ=25)                                         tCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaaa                                                  <  2:3008766/101‑1 (MQ=25)                                         tCGTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaaa                                                  <  2:2422171/101‑1 (MQ=25)                                           gTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAgg                                                <  1:597427/101‑1 (MQ=25)                                           gTCAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAgg                                                <  1:2963271/101‑1 (MQ=25)                                             cAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GGTCTGCCCCCCCCTC.C.TGCGCCTTTACCCCCCCGGGTAATTTTGCCCAGGTGAAAAGGAAAAaggg                                              >  1:1469876/1‑100 (MQ=255)                                             cAAGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGc                                              <  1:545474/101‑1 (MQ=25)                                              aaGGCCGCGATCGAGTGACAC.GGGCGCA.GG.G.GCGCTTTG.GCCCCCTCTGG.TGCGGCTTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCg                                             <  2:3790504/101‑1 (MQ=25)                                               aGGCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgc                                            >  1:449729/1‑100 (MQ=25)                                                ggCCACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgct                                           <  1:2416102/101‑3 (MQ=25)                                                  ccACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctcc                                         <  2:4619031/101‑5 (MQ=25)                                                   cACGATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctccg                                        <  2:3430890/101‑6 (MQ=25)                                                      gATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.g                                    <  1:3637734/101‑1 (MQ=25)                                                      gATCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCACCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.g                                    >  2:1599778/1‑101 (MQ=21)                                                       aTCGAGGGACAC.GGGCGGAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gg                                   <  2:2547234/101‑1 (MQ=21)                                                        tCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.ggg                                  <  1:4039596/101‑1 (MQ=25)                                                        tCGAGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.ggg                                  <  2:969819/101‑1 (MQ=25)                                                           aGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggc                               <  2:3787784/101‑3 (MQ=21)                                                           aGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggc                               <  1:3708239/101‑3 (MQ=21)                                                           aGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggc                               <  1:3756269/101‑3 (MQ=21)                                                           aGTGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggc                               <  1:3682779/101‑3 (MQ=21)                                                             tGACAC.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggcgc                             <  2:4089810/101‑5 (MQ=21)                                                              gACAC.GGGCGCAGGG.G.G.GCTCTGCCCCCCC.CT.GCTGCGCATTCACCCCCCGGGAGATTTTCGCCCAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggcgcc                            >  2:2722341/1‑96 (MQ=21)                                                                 ac.GGGCGCAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCtt                          <  1:3290981/101‑1 (MQ=21)                                                                      gcgcAGGG.G.G.GCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCTCCCTTTCAt                      <  1:3163213/101‑1 (MQ=18)                                                                       cgcaggg.g.g.gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATc                     <  2:4546926/96‑1 (MQ=21)                                                                       cgcaggg.g.g.gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATc                     <  2:2496662/96‑1 (MQ=21)                                                                         caggg.g.g.gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCtt                   <  1:3470676/98‑1 (MQ=21)                                                                          aggg.g.g.gCTCTGCCCCCTCT.G.C.TGCGCCTTCACCCCCCCCGCCATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGc                 >  2:2043596/3‑101 (MQ=255)                                                                          aggg.g.g.gCTCTG.CCCCCCTCT.GCTGCGCATTCGCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTg                  <  2:1218018/99‑1 (MQ=21)                                                                            gg.g.g.gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGcc                >  1:3437789/1‑101 (MQ=21)                                                                             g.g.g.gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCca               >  1:4433538/1‑101 (MQ=25)                                                                                 g.gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCcaca             <  2:2706470/101‑1 (MQ=25)                                                                                   gCTCTG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaa            >  1:170574/1‑101 (MQ=25)                                                                                    ctctG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaaa           <  2:4272398/101‑1 (MQ=25)                                                                                    ctctG.CCCCCCTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaaa           >  1:54755/1‑101 (MQ=25)                                                                                      tctgccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCCTCTTGCCACAAat          >  2:540449/5‑101 (MQ=25)                                                                                      tctgccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAat          >  1:715393/5‑101 (MQ=25)                                                                                      tctgccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAat          <  1:4551485/97‑1 (MQ=25)                                                                                      tctgccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAat          >  2:1630873/5‑101 (MQ=25)                                                                                       ctgccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAata         <  1:4664263/98‑1 (MQ=25)                                                                                        tgccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAatat        >  2:2957751/3‑101 (MQ=25)                                                                                         gccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATc       >  1:2297257/2‑101 (MQ=37)                                                                                          ccccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATcc      <  2:3944307/101‑1 (MQ=37)                                                                                           cccccTCT.GCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATccc     <  1:4233699/101‑1 (MQ=37)                                                                                           cccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATccc     <  2:3979555/101‑1 (MQ=37)                                                                                           cccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATccc     >  2:120312/1‑101 (MQ=37)                                                                                            ccccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATCCCg    <  1:221905/101‑1 (MQ=37)                                                                                             cccTCT.GCTGCGCATTCACCCCCCGGGATATTTACGCCAAGGTGAAATGGAAAAGGGGGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATCCCgg   <  1:3987956/101‑1 (MQ=37)                                                                                              ccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATCCCggg  <  2:4587322/101‑1 (MQ=37)                                                                                              ccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATCCCggg  <  1:3739770/101‑1 (MQ=37)                                                                                              ccTCT.GCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAAGGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACAAATATCCCggg  >  2:4030669/1‑101 (MQ=37)                                                                                                     |                                                                                               GCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACAC.GGGCGCA.GG.G.G.GCTCTGCCCCCCCTCT.G.TGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCG.TCCGCT.GGGGCGCCCTTTCATCTTGCCACAAATATCCCGGG  >  NC_014034/1843276‑1843461 Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG