Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,843,348	+G	100%	intergenic (+10/+138)	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,348	1	.→G	100.0%	251.2	77	intergenic (+10/+138)	ppiB/tyrS	peptidyl‑prolyl cis‑trans isomerase B/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  new base (4/72):  ref base (0/0):  total (4/73)

TCTGGTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.G.GGGGCTCT.GCCCCCCCTCTG.TGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCG.TCCGCT.GGGGCGCCCTTTCA  >  NC_014034/1843251‑1843440                                                                                                    |                                                                                                  tCTGGTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gg                                                                                               <  1:2581087/101‑1 (MQ=255) tCTGGTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gg                                                                                               <  1:664602/101‑1 (MQ=255)   tGGTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggg                                                                                             <  2:3953411/101‑1 (MQ=255)    ggTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.ggggg                                                                                            <  2:846027/101‑2 (MQ=255)    ggTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.ggggg                                                                                            >  1:3706958/1‑100 (MQ=255)     gTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggc                                                                                           <  2:4191662/101‑3 (MQ=255)     gTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggc                                                                                           >  2:1649533/1‑99 (MQ=255)      tCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggct                                                                                          <  2:2763421/101‑4 (MQ=255)      tCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggct                                                                                          <  2:2153724/101‑4 (MQ=255)      tCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggct                                                                                          <  2:1387589/101‑4 (MQ=255)       cGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggctc                                                                                         <  2:3536785/101‑5 (MQ=255)       cGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggctc                                                                                         <  2:835696/101‑5 (MQ=255)       cGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggctc                                                                                         <  1:3985357/101‑5 (MQ=255)        gACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.g.gggggctct                                                                                        <  1:3363057/101‑6 (MQ=255)         aCAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.g                                                                                       <  1:647064/101‑1 (MQ=255)            aGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gccc                                                                                    <  1:1904995/101‑1 (MQ=255)             gCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccc                                                                                   <  1:3935128/101‑1 (MQ=255)             gCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccc                                                                                   <  2:3543973/101‑1 (MQ=255)              cTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gccccc                                                                                  <  1:2493233/101‑1 (MQ=255)                gAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccct                                                                                <  2:3043058/101‑2 (MQ=255)                 aagaggggCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctc                                                                               <  1:620350/96‑3 (MQ=255)                 aaaaaGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctc                                                                               <  2:2570047/101‑3 (MQ=255)                  aaaaGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctct                                                                              <  1:197492/101‑4 (MQ=255)                    aaGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgc                                                                            <  2:3458600/101‑6 (MQ=255)                     aGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgct                                                                           <  1:3522693/101‑7 (MQ=255)                      gggCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctg                                                                          <  1:2331650/101‑8 (MQ=255)                      gggCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctg                                                                          <  1:1649648/101‑8 (MQ=255)                      gggCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCGCG.Gccccccccccccc                                                                          >  1:3137940/1‑95 (MQ=255)                       ggCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctgc                                                                         <  2:3830020/101‑9 (MQ=255)                       ggCGCCCCCGAGCTGAACGGCCTTGTCGAAGGCCCCGATTACATCGTCAAGGCCACGATCGAGTGACTCGGGCG.CAGG.GGGGCTCT.Gcccccctctgctgc                                                                         <  1:4168007/101‑9 (MQ=255)                        gcgcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctgcg                                                                        <  1:711157/101‑10 (MQ=255)                        gcgcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctgcg                                                                        <  1:4061688/101‑10 (MQ=255)                        gcgcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctgcg                                                                        <  1:4007934/101‑10 (MQ=255)                          gcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCa                                                                     <  1:277503/101‑1 (MQ=255)                             cccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTc                                                                  <  2:1894959/101‑1 (MQ=255)                              ccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCGCCC.G.GTGGCGCTGGCCCCCCctctccccccctca                                                                   >  2:2166837/1‑91 (MQ=255)                              ccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCa                                                                 <  2:636661/101‑1 (MQ=255)                              ccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCa                                                                 <  2:472628/101‑1 (MQ=255)                               cGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCAc                                                                <  1:1180054/101‑1 (MQ=255)                                gAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCAcc                                                               <  1:2357336/101‑1 (MQ=255)                                 aGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCAccc                                                              <  2:1521522/101‑1 (MQ=255)                                   cTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCAccccc                                                            <  1:2450845/101‑1 (MQ=255)                                   cTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCAccccc                                                            <  2:67092/101‑1 (MQ=255)                                       aCGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCggg                                                        >  2:3560679/1‑101 (MQ=255)                                          gCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGata                                                     <  2:3644915/101‑1 (MQ=38)                                             ttGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATAttt                                                  <  2:1556111/101‑1 (MQ=37)                                              tGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATAtttt                                                 <  2:3758077/101‑1 (MQ=37)                                               gTCGAAGACCCCGATTTCATCGGCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCCCCCCCCGGGATATTTTc                                                <  2:461569/101‑1 (MQ=25)                                                tCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCg                                               <  1:1436576/101‑1 (MQ=37)                                                tCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCg                                               <  2:3472161/101‑1 (MQ=37)                                                    aGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCaa                                           <  1:435492/101‑1 (MQ=37)                                                     gACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAg                                          <  1:2690792/101‑1 (MQ=37)                                                       ccccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGt                                        <  1:729955/101‑1 (MQ=37)                                                        cccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTg                                       <  1:1782448/101‑1 (MQ=37)                                                        cccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGATCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAggag                                       >  2:3116658/1‑99 (MQ=25)                                                         ccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGa                                      <  1:1490545/101‑1 (MQ=25)                                                         ccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGa                                      <  1:2664365/101‑1 (MQ=25)                                                         ccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGa                                      <  1:1271243/101‑1 (MQ=25)                                                          cGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaa                                     <  2:3213368/101‑1 (MQ=25)                                                          cGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaa                                     <  1:3001069/101‑1 (MQ=25)                                                           gATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaaa                                    <  1:113669/101‑1 (MQ=25)                                                           gATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaaa                                    <  2:255021/101‑1 (MQ=25)                                                           gATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.GCGGGGCTCT.G.CCCCCCTCTGCTGCGCCTCCTCCCCCCCGCCCACTTTCGCCAACGTGaaa                                    >  2:2158679/1‑101 (MQ=255)                                                            aTTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCGCCCCCCGGGATATTTTCGCCAAGGTGAAAt                                   <  2:2773259/101‑1 (MQ=25)                                                            aTTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAAt                                   <  2:463467/101‑1 (MQ=25)                                                             ttACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATg                                  <  1:333997/101‑1 (MQ=25)                                                              tACATCGTCAAGGCCCCGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATgg                                 <  2:3339214/101‑1 (MQ=25)                                                                 aTCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaa                              <  1:3409060/101‑1 (MQ=25)                                                                 aTCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaa                              <  1:47970/101‑1 (MQ=25)                                                                  tCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaaa                             <  2:3875275/101‑1 (MQ=25)                                                                  gcgtggcggCCACGATCGAGTGAGACGGGCG.CAGG.GGGGCGCT.G.CCCCCCTCTGCTGCGCCTTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaaa                             <  1:1271752/94‑1 (MQ=18)                                                                   cGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAg                            <  1:1041243/101‑1 (MQ=25)                                                                    ggcAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAgg                           <  2:3072138/99‑1 (MQ=25)                                                                       aaGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCg                        <  2:350191/101‑1 (MQ=25)                                                                        aGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgc                       <  1:2714826/101‑2 (MQ=25)                                                                         gtccACGATCGAGTGACACGTGCG.CATG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgct                      <  1:1018393/99‑3 (MQ=21)                                                                           ccACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctcc                    <  2:252900/101‑5 (MQ=25)                                                                           ccACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctcc                    <  1:1434091/101‑5 (MQ=25)                                                                           ccACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctcc                    <  1:2001518/101‑5 (MQ=25)                                                                            cACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctccg                   <  2:2887677/101‑6 (MQ=25)                                                                             aCGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGc                  <  2:2415826/101‑1 (MQ=25)                                                                                  cGAGTGACACGGGCG.CAGG.GGGGCTCT.GCCCCCCC.CTGCTGCGCATTCAACCCCCGGGATATTTTCGCCAAGGTGAAATCGAAAAGGGGGCTCCGCT.gggg            >  2:3224611/1‑101 (MQ=21)                                                                                  cGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggg            <  2:3267290/101‑1 (MQ=25)                                                                                  cGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggg            <  1:2603355/101‑1 (MQ=25)                                                                                    aGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATCTTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggc          <  2:1628834/101‑3 (MQ=21)                                                                                    aGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggc          <  2:1613083/101‑3 (MQ=21)                                                                                      tGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggcgc        <  1:3696712/101‑5 (MQ=21)                                                                                       gACACGGGCG.CAGG.GGGGCTCT.G.CTCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCT.gggggcgcc       >  2:3796306/1‑96 (MQ=21)                                                                                          ccGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCtt     <  2:862563/100‑1 (MQ=18)                                                                                          acGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCtt     <  1:3532560/101‑1 (MQ=21)                                                                                           ggggCG.CAGG.GGGGCTCT.G.CCCCCGTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCttt    <  2:2395717/100‑1 (MQ=18)                                                                                           cGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCttt    <  2:186366/101‑1 (MQ=21)                                                                                           cGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCttt    <  2:3729234/101‑1 (MQ=21)                                                                                           cGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCttt    <  2:3050912/101‑1 (MQ=21)                                                                                            gggCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTc   <  1:1659703/101‑1 (MQ=21)                                                                                             ggCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCa  <  1:2626057/101‑1 (MQ=21)                                                                                                    |                                                                                                  TCTGGTCGACAAGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.G.GGGGCTCT.GCCCCCCCTCTG.TGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCG.TCCGCT.GGGGCGCCCTTTCA  >  NC_014034/1843251‑1843440 Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG