Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,721,519	+GG	100%	intergenic (‑414/+32)	RCAP_rcc01588 ← / ← RCAP_rcc01589	GSCFA family protein/ErfK/YbiS/YcfS/YnhG family protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,721,519	1	.→G	100.0%	212.3	74	intergenic (‑414/+32)	RCAP_rcc01588/RCAP_rcc01589	GSCFA family protein/ErfK/YbiS/YcfS/YnhG family protein
Reads supporting (aligned to +/- strand):  new base (36/36):  ref base (1/0):  total (38/36)
*	NC_014034	1,721,519	2	.→G	100.0%	218.4	74	intergenic (‑414/+32)	RCAP_rcc01588/RCAP_rcc01589	GSCFA family protein/ErfK/YbiS/YcfS/YnhG family protein
Reads supporting (aligned to +/- strand):  new base (37/36):  ref base (1/0):  total (38/36)

TCGGGTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGCAGAGCGGATAAAGTTCCTCGACC  >  NC_014034/1721419‑1721605                                                                                                         ||                                                                                            tcgGGTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACccgcc                                                                                               >  1:4013403/1‑101 (MQ=255)  cgGGTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACccgcc...g                                                                                           >  2:6010889/1‑101 (MQ=255)  cgGGTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACccgcc...g                                                                                           >  1:5367520/1‑101 (MQ=255)   gggTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...gg                                                                                          >  2:3311972/1‑101 (MQ=255)   gggTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...gg                                                                                          >  2:4138225/1‑101 (MQ=255)    ggTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...ggg                                                                                         >  2:3125433/1‑101 (MQ=255)    ggTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...ggg                                                                                         >  1:5992864/1‑101 (MQ=255)    ctttCCAAACAACCGCGCGGGGCCCCCCCCCCCCGCTCCGTTTTTTTTTGGGCGGGTCGCTTGGGGGTTTC.CGCCTC.GGCTGCGCCGAGGCGCCCCGCC...gg                                                                                          <  2:1369518/99‑1 (MQ=9)      ttCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...ggggg                                                                                       >  1:6069440/1‑99 (MQ=255)      ttCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...ggaga                                                                                       >  1:1369569/1‑98 (MQ=255)        ccAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACccgcc...gcggggg                                                                                     >  1:5923400/1‑95 (MQ=255)        ccAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...gggggct                                                                                     >  1:2277713/1‑97 (MQ=255)           aaTAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...ggcggctgcg                                                                                  >  2:5365579/1‑93 (MQ=255)            aTAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...ggggggggcgc                                                                                 >  1:1941321/1‑93 (MQ=255)             tAA.CGCGCTGCGGGCCTCCCCCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...gggggctgcgcc                                                                                >  2:4348565/1‑92 (MQ=37)              aa.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...GGGGG.GCGcgccc                                                                              >  1:4464478/1‑101 (MQ=255)                  gcgcTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCGCCCg                                                                             >  2:4528072/1‑101 (MQ=25)                   cggtGCGGCCCCCCCCCCCGCCTCTTTTTTTTTTGGGCCGGTCGCCTTGGGGTTTC.CGCCTC.GGCGGCGCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCgg                                                                            <  2:46444/98‑1 (MQ=255)                    gatGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCggg                                                                           >  2:4524291/3‑100 (MQ=25)                      tGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..cc                                                                       >  2:6444945/1‑101 (MQ=255)                       gCGGCCCCCCCCCCCGCTGCTTTTTTTTTTGGGCCGGTCGCTTTGGGGTTTT.CCCCGC.GGCCGCCCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..ccc                                                                      <  1:5449330/101‑1 (MQ=255)                       gCGGCCCCCCCCCCCCCGGCTTTTTTTTTTGGGCGGGTCGCCTGGGGTATTC.CGCCGC.GGCGGCGCGGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..ccc                                                                      <  2:6831088/101‑1 (MQ=255)                        cGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..cccc                                                                     >  2:6811243/1‑101 (MQ=255)                          ggCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCCGA..cccccc                                                                   >  1:3866120/1‑101 (MQ=25)                          gcccccccccccGCGGCGTTTTTTTGTTGGCCGGCCCCCTGGGGGTTTC.CCCCTC.GGCGGCGCGGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..cccccc                                                                   <  1:414401/96‑1 (MQ=255)                            ccccccccccGCTGCGCTTTTTTTTTGGCGGGCCCCGTTGGGGTTTC.CGCCTC.GCCGCCCCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..cccccccc                                                                 <  1:3791857/98‑1 (MQ=9)                             ctccctccCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCg                                                                >  2:2739686/1‑101 (MQ=255)                             cccccccccGCTGCTTTTTTTTTGGGGCCGGTCGCTTGGGGGTTTC.CCCCCC.GGCGGCCCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCg                                                                <  1:273028/99‑1 (MQ=255)                             cccccccccGCGGCGTTTTTTTTGGGGCGGGCCGTCTGGGGGTTTC.CGCCGG.GGCGGCGCGGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCg                                                                <  2:4823021/99‑1 (MQ=255)                              tccctcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCGGG..cccccccccc                                                               >  2:3791696/1‑99 (MQ=21)                               ccctcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCGGG..CCCCCCCCCaa                                                              >  1:1875450/1‑101 (MQ=21)                               ccctcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGG..cccccccccga                                                              >  2:5069602/1‑98 (MQ=21)                               ccctcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCCaa                                                              >  1:6821110/1‑101 (MQ=21)                               ccctcccGCAGCGCTTTATTTTTGGCCGGCCGCCTTGGGGATTT.CCCCTC.GGCTGGGCCGAGGCGACCCGCCGG.GGGGG.GCCCCCCCGCA..CCCCCCCCCaa                                                              >  2:4078108/1‑101 (MQ=21)                                cctcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GGGCCCCCCGG..ccccccccccca                                                             >  1:2750312/1‑97 (MQ=21)                                ccccccGCTGCTTTTTTTTTTGGGCGGGTCGCCTGGGGGTTTC.CCCCCC.GGCGCCGCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAg                                                             <  2:5899149/98‑1 (MQ=255)                                ccccccGCTGCGTTTTTTTTGTGGCCTGGCGCCTGGGGGTTTC.CGCCCC.GGCGGCGCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAg                                                             <  2:1567174/98‑1 (MQ=255)                                 ctcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGC..CCCCCCCCCAAGa                                                            >  2:516524/1‑101 (MQ=21)                                  tcccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCCGC..CCCCCCCCGAAgca                                                           >  2:474099/1‑99 (MQ=21)                                   cccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGGGGGGCG.GCGCCCCCCGC..CCCCCCCCCAAGac                                                           >  1:2049515/1‑100 (MQ=21)                                   cccGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCCGG..CCCCCCCCGAAGagc                                                          >  2:3100073/1‑99 (MQ=21)                                    ccgcctcgttttTTTTTGGGCGGGCCGCGTGGGGGTTTC.CCCCCC.GGCGGCCCCGGGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGAtct                                                         <  1:6177499/92‑1 (MQ=255)                                     cGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCCGG..CCCCCCCCCCAGAtctc                                                        >  2:3164264/1‑101 (MQ=21)                                      gCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCGGG..CCCCCCCCCAAGAtctcc                                                       >  1:2925814/1‑100 (MQ=21)                                      gCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCCGG..CCCCCCCCCAAGATCTCa                                                       >  1:1708601/1‑101 (MQ=21)                                      gCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCa                                                       >  2:7104671/1‑101 (MQ=255)                                      cgggcgtTTTTTTTGGGGCCGGCCGCCTGGGGTTTCC.CCCCGC.GGCGGCGCCGGCGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCa                                                       <  1:2222058/98‑1 (MQ=255)                                       cAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCCCAGATCTCAg                                                      >  2:6780166/1‑101 (MQ=21)                                       cAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCCG.GGGCTGGGGCCCCCGGC..CCCCCCCCgcgatctccc                                                       >  1:5244726/1‑92 (MQ=21)                                        aGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCCGG..CCCCCCCCCGCGATCTCCGa                                                     >  2:4077662/1‑101 (MQ=18)                                          gtttttGTTTTGGCCCGTTCGCGTTGGGGTTTC.CCCCCG.GGCGGCGCCGAGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATc                                                   <  2:2944451/98‑1 (MQ=255)                                          cgcTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCCGG..CCCCCCCCCCAGAgctcccagc                                                   >  1:6245541/1‑96 (MQ=18)                                           tttttTTTTTGGGCGGGCCGCTGTGGGGTTTG.CCCCGC.GGCGGCGCCGCGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCa                                                  <  2:5708270/99‑1 (MQ=12)                                           tttttTTTTTGGGCCGGTCGCCTGGGGGTTTC.CCCCCC.GGCGGCGCCGGGCCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCa                                                  <  1:4393628/99‑1 (MQ=12)                                             tttgtttggggCCGGTCGCCTGGGGTTTTC.CGCCGG.GGCCGCGCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.cc                                               <  1:3351083/92‑1 (MQ=255)                                              tttttttgGGCCGGGCGCCGTGGGGTTTC.CCCCCC.GGCGGCGCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCt                                              <  1:1365891/98‑1 (MQ=12)                                              tttttggtGGGCGGTCGCTGGGGGGATTC.CCCCCG.GGCGCCCCCGCGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCt                                              <  2:1637500/94‑1 (MQ=9)                                              ttttgtgtGGCCGGCCGCGTGGGGCTTTC.CGCCCC.GGCGGCGCCGCGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCt                                              <  2:5574641/94‑1 (MQ=12)                                               tttttggggCCGGGCGCGGTGGGGTTTC.CGCCGC.GGCGGCGCCGGGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCtt                                             <  1:2757612/94‑1 (MQ=255)                                               ttttgtgggCCGGTCGCCTGGGGGTTTC.CGCCGC.GGCGGCGCCGAGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCtt                                             <  2:2406478/94‑1 (MQ=255)                                                ttttttGGGCGGGCCGCTTGGGGGTTTCCCCCGC.GGCGGCCCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCtt                                             <  1:6167670/100‑1 (MQ=255)                                                ttttttGGCGGTCCGCGTGGGGTCTTC.CGCCGC.GGCGCCGCGGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTg                                            <  2:2649296/100‑1 (MQ=255)                                                   tttGGCCGGCCGCCGGGGGGTTTC.CGCCCC.GGCCGCGCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGAcc                                         <  1:1863082/101‑1 (MQ=255)                                                   ttgggCCGGGCGGCTTGGGGGTTT.CCCCCCGGGCGGCGCGGGGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGAc                                          <  1:4503109/98‑1 (MQ=255)                                                    ggggCGGGTCGCCGTGGGGTTTC.CCCCGC.GCCGCCGCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCt                                        <  1:3753928/99‑1 (MQ=16)                                                     tGGCCGGTCGCGTGGGGGTTTC.CGCCGC.GCCGGCGCGGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTg                                       <  1:4356023/101‑1 (MQ=255)                                                      ggCGGGTCGCGGGGGGGTTTC.CCCCTG.GGCCGCGCCGGGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGc                                      <  1:6243497/101‑1 (MQ=255)                                                      ggCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGGGGGGGG.GGGCCCCCGGC..CCCCCCCCCAAG.TCTCAGACTA.CCCTTACCtcc                                      >  2:2720943/1‑99 (MQ=9)                                                        ccGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCCAAGATCTCAGATCA.CCTTGACCTGCgg                                    >  1:3133285/1‑100 (MQ=21)                                                         cGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCCGG..CCCCCCCCCCCAATCTCCGATCACCCTTGACCTCCgg                                    >  2:6156591/1‑100 (MQ=9)                                                          ggcggcgtgggggTTTC.CCCCCC.GGCGGCGCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGc                                  <  1:3312088/92‑1 (MQ=16)                                                          ggTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGGCCCCCGG..CCCCCCCCCACGATCCCCGAAGC.CCTTGACCCGCGTGc                                  >  2:5603507/1‑101 (MQ=9)                                                            gggggtgggggCTTC.CGCCGC.GCCCGCCCCGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCg                                <  1:219124/94‑1 (MQ=255)                                                              gCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCCGA..CCCCCCCCCACGAGCTCAGATCC.CCCTGACCCGCGTggccgc                              >  1:354760/1‑96 (MQ=9)                                                              gCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCCCAGATCTCAGATCA.CCTTGACCTGCGTGCCGAc                              >  2:3104035/1‑101 (MQ=21)                                                               gcgtGGGGATTT.CCCCCCGGGCGCCGCGGAGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGAc                              <  2:3626270/98‑1 (MQ=255)                                                               cttgggggGTTC.CCCCTC.GGCCGCGCGGGGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGAcc                             <  2:7012947/97‑1 (MQ=255)                                                                cTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCg                            >  1:2219352/1‑101 (MQ=255)                                                                 gtGGGGGTTC.CCCCCC.GGCCGCGCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCgg                           <  1:6811811/100‑1 (MQ=255)                                                                  tGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGGT.GCGCCCCCCGG..CCCCCCCCCAAGATCTCAGAGTA.CCCTGACCCGCGTGCCCGCCggg                          >  2:2555730/1‑101 (MQ=12)                                                                   ggggCTTT.CCCCTG.GGCGCCGCCGCGCCCCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGc                         <  1:4328247/101‑1 (MQ=21)                                                                   ggggATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.CGGCCCCCCGA..CCCCCCCCCAAGATCTCAGAACA.CCCTGACCCTCGTGcccgccccgc                         >  2:5487667/1‑93 (MQ=9)                                                                        ttt.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCCGATCA.CCTTGACCTGCGTGCCGACCGGGCCGAGc                    >  2:3329249/1‑101 (MQ=21)                                                                         tt.CGCCTC.GGCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGG..CCCCCCCCCCAGAGCCCAGCTCC.CCCTGCACTGCGGGCCGACCGGGCggaggg                   >  2:1764577/1‑99 (MQ=2)                                                                            gcccgc.GGCGGCGCGGCGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGCAGAGCGGa                 <  1:3430281/96‑1 (MQ=255)                                                                              ccgc.GGCCGCGCCGCGGCCCCCCGCCGG.GGGCT.GCGCCCCCGGACCCCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGCAGAGCGGa                 <  1:1977786/98‑1 (MQ=255)                                                                                  gggCGGCGCCGGGGCGCCCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGCAGAGCGGATAAAg            <  1:215054/100‑1 (MQ=255)                                                                                   ggCTGCGCCGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGAACCCAGATCA.CCTTTACCCGCGCGCCCGCCCGGCAGCGGGGGTAAAGt           >  1:2247596/1‑101 (MQ=12)                                                                                           cGAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCCAAGATCTCAGATCA.CCTTGACCTGCGTGCCGCCCGGGCAGAGCGGATAAAGTTCCTCGAc   >  1:1235986/1‑101 (MQ=37)                                                                                            gAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGCAGAGCGGATAAAGTTCCTCGaac  >  2:1024851/1‑99 (MQ=38)                                                                                            gAGGCGACCCGCCGG.GGGCT.GCGCCCCCGGA..CCCCCCCCCAAGATCTCAGATCA.CCTTGACCTGCGGGCCGACCGGGCAGAGCGGATAAAGGTCCCCGAcc  >  2:5427673/1‑101 (MQ=37)                                                                                                         ||                                                                                            TCGGGTTCCAAATAA.CGCGCTGCGGGCCTCCCTCCCGCAGCGCTTTATTTTTGGCCGGTCGCCTTGGGGATTT.CGCCTC.GGCTGCGCCGAGGCGACCCGCC...GGGCT.GCGCCCCCGGA..CCCCCCCCGAAGATCTCAGATCA.CCTTGACCTGCGTGCCGACCGGGCAGAGCGGATAAAGTTCCTCGACC  >  NC_014034/1721419‑1721605 Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 27 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG