Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,732,868	+G	4.9%	intergenic (‑292/‑111)	RCAP_rcc01601 ← / → tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,732,868	2	.→G	4.9%	6.3 0.0 ‑1.7	81	intergenic (‑292/‑111)	RCAP_rcc01601/tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (12/65):  total (12/69)
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.14e-02

TCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GT...C...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCGGCG  >  NC_014034/1732788‑1732949                                                                                               |                                                                                                tCAACAGGCAGGGGCCCCGCTTTTCCTAGGCAAAATCCGCCTTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GC...GGGTCG.CCTC.GGGCTTG.ccc                                                                        <  2:2033502/101‑1 (MQ=21) tCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GCGGGT...CG.CCTC.GGGCTTG.ccc                                                                        <  2:2757771/101‑1 (MQ=25)  cAACCGGCAAGGGCCCCCCTTTTCCTAGGCAAAATACCCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCg                                                                       <  1:4151551/101‑1 (MQ=18)  cAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTT.GGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GCGGGT...CG.CCTC.GGGCTTG.CCCg                                                                       <  1:7069156/101‑1 (MQ=21)   aaCAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGa                                                                      <  1:5640210/101‑1 (MQ=21)      aGGCAAGGGCCCCGCTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GCGGGT...CG.CCTC.GGGCTTG.CCCGAGGc                                                                   <  1:5649749/101‑1 (MQ=21)      aGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGc                                                                   <  1:4768111/101‑1 (MQ=21)       ggCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCg                                                                  <  2:2976965/101‑1 (MQ=21)          acgggCCCCGCTTTTCCTGGAGAAAAAACCGGCTCTTTTCTTTG.GGGGGGGGCGGGGGTGGAGC...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGaaa                                                               <  1:7117749/99‑1 (MQ=2)             ggCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GCGGGT...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.t                                                         <  2:549009/101‑1 (MQ=18)               cccgGCTTCCCCAAGGCAAACGACGCTCTTCTTTCTTTC.GGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.c                                                      <  2:6875711/101‑1 (MQ=2)                cccGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GC...GGGTCG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.cc                                                     <  1:4999542/101‑1 (MQ=21)                 ccGCTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GC...GGGTCG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.ccc                                                    <  2:4944780/101‑1 (MQ=18)                   gCTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.ccccc                                                  <  1:2940803/101‑1 (MQ=18)                     ttttCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCg                                                <  1:2165169/101‑1 (MQ=21)                     tgttCCTAGGAAAACCACGCCATGCTTTCTTTT.CGGGGGGGCGGGGGGTAACG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCg                                                <  2:5631636/99‑1 (MQ=9)                      tttCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGt                                               <  1:1193580/101‑1 (MQ=21)                          cTAGGCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGaaa                                           <  2:4413121/101‑1 (MQ=18)                          cTAGGCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGaaa                                           <  1:2763998/101‑1 (MQ=18)                          cTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGAAGAC.CCCC...C..CG...GG...G...CG.CCCC.GGGCGTG.CCCCCGGCGGAA.G.CGTT.cccccccgcaa                                            >  2:4421561/1‑96 (MQ=255)                            aGGCAAAATCCGCCCTTCTTTCTTTT.GGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAAc                                         <  1:6231759/101‑1 (MQ=18)                            aGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAAc                                         <  2:4964422/101‑1 (MQ=21)                             ggCAAAAGACCCCTTTTTTTCTTTG.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAAcc                                        <  2:3042503/101‑1 (MQ=17)                              gCAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCt                                       <  1:5345149/101‑1 (MQ=18)                              gCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCt                                       <  1:4575363/101‑1 (MQ=21)                               cAAGATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGAAG...CCCC...C..CGCG.GG...G...CG.CCCC.GGGGTTG.CCCGCGGCGGAA.G.C.TT.CCCCCCCCGACAAccct                                      >  1:6849523/1‑99 (MQ=12)                               cAAAATCCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCtt                                      <  2:3048690/101‑1 (MQ=18)                               cAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCtt                                      <  2:4442504/101‑1 (MQ=21)                                aaactgcGCCCTTCTTTGTTTG.GGGGGGGGGGGGGGGGGACG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTg                                     <  1:6238158/95‑1 (MQ=9)                                aaaaTACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTg                                     <  2:5825812/101‑1 (MQ=21)                                 aagtgacgccTTCTTTCTTTG.CGGGGGGGCGGGGGTTGAAG...CCCCTCGC..G....GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGc                                    <  2:1852828/93‑1 (MQ=17)                                 aaaaccccccTTTTTTTTTTT.GGGGGGGGGGGGGGGGGAAC...CCCC...C..CCGCGGG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTg                                     <  1:3359867/94‑1 (MQ=9)                                 aaaTACGCCCTTCTTTCTTTG.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGc                                    <  2:2936750/101‑1 (MQ=18)                                  attccGCCTTTTTTTCTTTG.GGGGGGGGGGGGGGGTGAGC...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGcc                                   <  2:2410244/97‑1 (MQ=9)                                  aatccccccTTCTTTCTTTG.GGGGGGGGCGCGGGGTGAAG...CCCC...CGGCG...GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGcc                                   <  2:5001854/95‑1 (MQ=17)                                   atccGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCt                                  <  2:2620196/98‑1 (MQ=18)                                   aTACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCt                                  <  2:4639020/101‑1 (MQ=21)                                   aTACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCt                                  <  2:6472310/101‑1 (MQ=21)                                    tACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGGAG...CCCC...C..CCCG.GG...G...GG.CCCCGGGGCTGC.CCCCGGGGGAAA.G.C.TT.CCCCCCGCGAAAAccctcccc                                  >  1:4338992/1‑95 (MQ=255)                                    tACGCCCTTCTTTCTTGT.CGGGGGGGGGGGGGGTGAAG...CCCC...CGGCG...GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCtc                                 <  2:6095208/101‑2 (MQ=18)                                     tctcccTTTTTTGTTTTGGGGGGGGGGGGGGGGTGACC...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTg                                 <  1:6857/98‑1 (MQ=9)                                     ccGCCCTTCTTTCTTTT.GGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTgc                                <  1:908665/100‑1 (MQ=18)                                     ccGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTgc                                <  1:6496571/100‑1 (MQ=18)                                     aCGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...G...CG.CCTC.GGGCTTGCCCCGGGGGGAAA.G.G.TC.CCCCCCG.GAAAACCCTGccccg                                >  1:5695201/1‑98 (MQ=9)                                      cGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTgcg                               <  2:5148133/101‑1 (MQ=21)                                        cccTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTgcgcg                             <  2:766313/101‑1 (MQ=21)                                           ttctttctTTT.CGGGGGGGCGGGGGGGGGAG...CCCC...C..CCCG.GG...G...GG.CCCC.GGGCCTG.CCCCAGGGGAAACG.T.TT.CCCCCCCGGAAACCCTTCccgtggccgg                           >  2:240673/1‑93 (MQ=255)                                            tctttcttTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTc                         <  1:7091998/101‑1 (MQ=21)                                            tctttcttTC.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTc                         <  1:2762994/101‑1 (MQ=18)                                            tctttctgTC.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTc                         <  2:6029221/101‑1 (MQ=18)                                             ctttctttT.CGGGGGGGCGGCGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCg                        <  1:5600068/101‑1 (MQ=21)                                              tttctttT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.c                      <  2:3491174/101‑1 (MQ=21)                                              tttctttT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.c                      <  2:1909029/101‑1 (MQ=21)                                              tttctgtG.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.c                      <  2:5634741/101‑1 (MQ=18)                                               ttctttT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.t                    <  1:1041612/101‑1 (MQ=21)                                               ttctttT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.t                    <  1:7137035/101‑1 (MQ=21)                                               ttctttT.CGGGGGGGCGGGGGGGGAAG...CCCC...C..CCCG.GG...G...CG.CCCC.GGGCGTG.CCCCCGGGGGAA.GCC.TT.CCCCCCCGGAAAA.CTTTCCTGCGGGCGTCG.cc                     >  2:1740682/1‑100 (MQ=255)                                                tctttG.CGGGGGGGCGGGGGGTGACG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.tt                   <  2:6075187/101‑1 (MQ=18)                                                 ctttT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.ttt                  <  1:6700531/101‑1 (MQ=21)                                                 ctttG.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.ttt                  <  2:5117901/101‑1 (MQ=18)                                                 ctttC.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.ttt                  <  2:5428037/101‑1 (MQ=18)                                                  ttttgggggggggCGGGGGGTGAAC...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.ttt                  <  2:193796/95‑1 (MQ=17)                                                  tttt.CGGGGGGGCGGGGGGGGAAGCCCCCCC...C..C....GG...G...CG.CCCC.CGGGCTG.CCCCACGCGGAA.G.A.TT.CCCCCCCCGACAAACCTGCCTGCCCCCGTCG.C.CTTc                 >  2:5804603/1‑101 (MQ=255)                                                   ttt.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTcc                <  2:3633527/101‑1 (MQ=21)                                                   ttt.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTcc                <  1:3091435/101‑1 (MQ=21)                                                     t.CGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCt              <  2:164089/101‑1 (MQ=21)                                                       cGGGGGGGCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCtt             <  1:1511676/101‑1 (MQ=21)                                                       cGGGGGGGCGGGGGGGGAAG...CCCC...C..CCCG.GG...G...CGCCCCC.GGGCTTG.CCCCGGGGGGAG.G.C.TT.CCCCCCCGGAAAACC.TGCCCCCGCGCGGCGCC.TTTCCCt              >  2:4218379/1‑101 (MQ=255)                                                         ggggggCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTc           <  1:2005453/101‑1 (MQ=21)                                                         ggggggCGGGGGGTGAAG...CCCC...C..CGCG.GG...G...CG.CCCC.GGGGTTG.CCCGGGGCGGAA.G.C.TTCCCCCCCGGGAAAACC.TGCCTGCGCGCCGCG.C.CTTCCCCTTc           >  2:2763853/1‑101 (MQ=255)                                                         ggggggCGGGGGGGGAAG...CCCC...C..CGCG.GG...G...CG.CCTC.GGGCTTG.CCCCCGGGGGAA.G.C.TT.TCCCCCCCGAAAACCCTGCCTGCGCGCGGCG.CGCTTCCCttc            >  2:68827/1‑100 (MQ=2)                                                          gggggGGGGGGGTGAAC...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCg          <  2:5895170/101‑1 (MQ=21)                                                          gggggCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCg          <  2:5068994/101‑1 (MQ=21)                                                          gggggCGGGGGGTGAAC...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCg          <  2:6458773/101‑1 (MQ=21)                                                           ggggCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCgg         <  1:4952976/101‑1 (MQ=21)                                                             ggCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGcc       <  2:6380384/101‑1 (MQ=21)                                                              gCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGccc      <  1:4235008/101‑1 (MQ=21)                                                              gCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGccc      <  2:1622281/101‑1 (MQ=21)                                                              gCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCACCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGccc      <  2:6075936/101‑1 (MQ=21)                                                              gCGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCCC.GGGCCTG.CCCGAGGGGGAA.A.C.CT.CCCCCCGTGAAAACCCTGCCCGCCCGCGGCG.C.TTTCCCTTTCCGccc      >  2:2229498/1‑101 (MQ=9)                                                              gCGGGGGGTGAAG...CCCC...C..CGCG.GG...G...CG.CCTC.GGGCTTG.CCCCAGGGGAAA.G.C.TT.CCCCCCGCGAAAACCTTGCCTGCGCGCGTCG.C.GTTCCCTTTTGGtcc      >  2:7046876/1‑101 (MQ=16)                                                               cGGGGGGTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.CTTCCCTTTCGGcccc     >  2:6227832/1‑100 (MQ=21)                                                               cGGGGGGGGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCg     <  2:2172497/101‑1 (MQ=21)                                                                ggggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGccccg    >  2:15301/1‑99 (MQ=21)                                                                ggggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGccccg    >  2:2801358/1‑99 (MQ=21)                                                                ggggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCgg    >  1:6850078/1‑101 (MQ=21)                                                                ggggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGGGAAAACCTTCCCT.CGCGCGGCG.C.CTTTCCTTTCGGGCCCgc   >  2:3049960/1‑101 (MQ=9)                                                                ggggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGCGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCgg    >  1:6706643/1‑101 (MQ=21)                                                                ggggggTGAAG...CCCC...C..CGCG.GG...G...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCgg    >  2:3986917/1‑101 (MQ=25)                                                                 gggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCGgc   >  1:592384/1‑101 (MQ=21)                                                                 gggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCGgc   >  1:768010/1‑101 (MQ=21)                                                                 gggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCGgc   >  2:4695729/1‑101 (MQ=21)                                                                  ggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTACCTTTCGGcccctcg  >  2:7048065/1‑97 (MQ=21)                                                                  ggggTGAAG...CCCC...C..CGCG.GG...T...CG.CCCC.GGGCTTG.CCCGCGGCGCAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.A.CTTCCCCTTCGGCCCgggc  >  2:2371003/1‑99 (MQ=16)                                                                                               |                                                                                                TCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCC...C..C....GT...C...CG.CCTC.GGGCTTG.CCCGAGGCGAAA.G.C.TT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCG.C.TTTCCCTTTCGGCCCGGCG  >  NC_014034/1732788‑1732949 Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 20 ≤  ATCG  < 32 ≤  ATCG  < 39 ≤  ATCG