Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,843,359	Δ1 bp	100%	intergenic (+21/+127)	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,359	0	C→.	100.0%	228.8	69	intergenic (+21/+127)	ppiB/tyrS	peptidyl‑prolyl cis‑trans isomerase B/tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  new base (7/61):  ref base (1/0):  total (8/61)

AGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.G.GGGGCTCT.GCCCCCCCTCTG.TGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCG.TCCGCT.GGGGCGCCCTTTCATCTTGCCACAAA  >  NC_014034/1843262‑1843452                                                                                                      |                                                                                                 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:1649648/101‑8 (MQ=255)           gggCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctg                                                                                      <  1:2331650/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:4007934/101‑10 (MQ=255)             gcgcgcCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.Gcccccctctgctgcg                                                                                    <  1:4061688/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:472628/101‑1 (MQ=255)                   ccGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCa                                                                             <  2:636661/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                                                           <  2:3472161/101‑1 (MQ=37)                                     tCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCg                                                           <  1:1436576/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:2664365/101‑1 (MQ=25)                                              ccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGa                                                  <  1:1490545/101‑1 (MQ=25)                                              ccGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGa                                                  <  1:1271243/101‑1 (MQ=25)                                               cGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaa                                                 <  1:3001069/101‑1 (MQ=25)                                               cGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaa                                                 <  2:3213368/101‑1 (MQ=25)                                                gATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaaa                                                <  2:255021/101‑1 (MQ=25)                                                gATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGaaa                                                <  1:113669/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:47970/101‑1 (MQ=25)                                                      aTCGTCAAGGCCACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGaaa                                          <  1:3409060/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:2001518/101‑5 (MQ=25)                                                                ccACGATCGAGTGACACGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCgctcc                                <  1:1434091/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:3050912/101‑1 (MQ=21)                                                                                cGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCttt                <  2:186366/101‑1 (MQ=21)                                                                                cGGGCG.CAGG.GGGGCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCttt                <  2:3729234/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)                                                                                      gcagg.ggggCTCT.G.CCCCCCTCTGCTGCGCATTCACCTCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGGGCTCCGCTGGGGGGGCCCTTTCATCt           >  1:2254495/5‑101 (MQ=21)                                                                                       cagg.ggggCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCtt          <  2:3943579/98‑1 (MQ=21)                                                                                        agg.ggggCTCT.G.CCCCCCGCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTg         >  2:207364/3‑101 (MQ=21)                                                                                            ggggCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCca      >  2:3391235/1‑101 (MQ=25)                                                                                            ggggCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCca      <  2:1135816/101‑1 (MQ=25)                                                                                             gggCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCcac     <  1:857396/101‑1 (MQ=25)                                                                                             gggCTCT.G.CCCCACTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCcac     <  1:952929/101‑1 (MQ=25)                                                                                              ggCTCT.G.CCCCCCTCTGCGGCGCCATCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCcaca    <  1:1651313/101‑1 (MQ=21)                                                                                               gCTCT.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaa   >  1:2107761/1‑101 (MQ=25)                                                                                                ctct.G.CCCCCCTCTGCTGCGCATTCCCCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaaa  <  1:1541924/101‑1 (MQ=25)                                                                                                ctct.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaaa  <  1:2583253/101‑1 (MQ=25)                                                                                                ctct.G.CCCCCCTCTGCTGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCGCTCCGCTGGGGGCGCCCTTTCATCTTGCCACaaa  <  2:2165508/101‑1 (MQ=25)                                                                                                      |                                                                                                 AGCTGAAAAAGGGCGCGCCCGAGCTGAACGGCCTTGTCGAAGACCCCGATTACATCGTCAAGGCCACGATCGAGTGACACGGGCG.CA.G.GGGGCTCT.GCCCCCCCTCTG.TGCGCATTCACCCCCCGGGATATTTTCGCCAAGGTGAAATGGAAAAGGGCG.TCCGCT.GGGGCGCCCTTTCATCTTGCCACAAA  >  NC_014034/1843262‑1843452 Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG