Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,255,904	+GT	100%	intergenic (+74/+168)	RCAP_rcc01184 → / ← RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,255,904	1	.→G	100.0%	137.3	82	intergenic (+74/+168)	RCAP_rcc01184/RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (36/46):  ref base (0/0):  total (36/46)
*	NC_014034	1,255,904	2	.→T	100.0%	139.5	82	intergenic (+74/+168)	RCAP_rcc01184/RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (36/46):  ref base (0/0):  total (36/46)

TGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGG....CTG...TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCAAATATCCCGGGGGGGGGGAACTTCCCG  >  NC_014034/1255828‑1256004                                                                                     ||                                                                                                     tGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTccccc                                                                                     >  2:2482499/1‑101 (MQ=9)  gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                                                                                    >  2:1214859/1‑101 (MQ=9)  gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                                                                                    >  2:4023558/1‑101 (MQ=9)   aCCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCg                                                                                   >  2:1741696/1‑101 (MQ=9)    ccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgg                                                                                  <  2:2266284/101‑1 (MQ=2)    ccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgg                                                                                  <  2:277196/101‑1 (MQ=2)    ccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgg                                                                                  <  2:4124869/101‑1 (MQ=2)     cTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCggg                                                                                 <  2:1347009/101‑1 (MQ=2)      tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                                                                <  1:4032098/101‑1 (MQ=2)      tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                                                                >  2:4260507/1‑101 (MQ=2)      tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGC..CCGCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGa                                                                                <  2:911090/101‑1 (MQ=2)       aGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGCCTCCCCCCGGGat                                                                               >  2:2089905/1‑101 (MQ=2)       aGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                                                               >  2:486652/1‑101 (MQ=2)       aGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                                                               >  2:633933/1‑101 (MQ=2)       aGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGat                                                                               <  2:374610/101‑1 (MQ=2)         gggTGCCGGACCGGACCCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGatat                                                                             <  1:4227425/101‑1 (MQ=0)         gggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGGCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGatat                                                                             <  2:2966379/101‑1 (MQ=0)         gggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGatat                                                                             <  1:473656/101‑1 (MQ=2)           gTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                                                           <  1:1481920/101‑1 (MQ=0)           gTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                                                           >  1:2077205/1‑101 (MQ=0)           gTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                                                           <  2:2422724/101‑1 (MQ=0)           gTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                                                           <  2:3015245/101‑1 (MQ=0)           gTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt                                                                           >  2:2587775/1‑101 (MQ=0)            tGCCGGCCCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc                                                                          <  2:1934760/101‑1 (MQ=2)            ggCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc                                                                          <  2:3523985/100‑1 (MQ=2)             gCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTcc                                                                         <  2:1711321/101‑1 (MQ=2)               cggaccggaACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc                                                                       <  2:1495320/101‑1 (MQ=2)               cggaccggaACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc                                                                       <  2:508957/101‑1 (MQ=2)               cggaccggaACCGGGGGCCTTCGGGACCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc                                                                       >  1:2699348/1‑101 (MQ=2)                ggaccggaACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc                                                                      <  2:1943684/101‑1 (MQ=2)                 gaccggaACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCa                                                                     <  2:1426592/101‑1 (MQ=2)                 gaccggaACCGGGGGCCTTCGGGCCCCCGC..TTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCa                                                                     >  1:2367910/1‑101 (MQ=2)                                         gcccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGgcgc                                               >  1:4056982/11‑101 (MQ=17)                                         gcccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGgcgc                                               <  1:3667969/91‑1 (MQ=17)                                         gcccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCCGGATATTTCCGCCAAGGGGAAATGAACCCGGGgcgc                                               >  2:818458/11‑101 (MQ=12)                                          cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                              <  1:4132193/92‑1 (MQ=17)                                          cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                              <  1:2335226/92‑1 (MQ=17)                                          cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                              <  2:1570050/92‑1 (MQ=17)                                          cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                              >  2:2987681/10‑101 (MQ=17)                                          cccccgctttcGGTCCGACCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                              >  2:2482111/10‑101 (MQ=17)                                           ccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCCAGGCGCgg                                             >  2:1097102/9‑101 (MQ=9)                                            cccgcttccGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGc                                            <  2:846604/93‑1 (MQ=17)                                               gctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCttt                                         >  2:2300133/5‑101 (MQ=18)                                               gctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCttt                                         >  2:2761367/5‑101 (MQ=18)                                               gctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCttt                                         >  1:994761/5‑101 (MQ=18)                                               gctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCttt                                         >  1:829417/5‑101 (MQ=18)                                                ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc                                        >  2:3056336/4‑101 (MQ=18)                                                ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc                                        <  2:3934153/98‑1 (MQ=18)                                                ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc                                        >  1:2324309/4‑101 (MQ=18)                                                 tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                       >  2:1276039/3‑101 (MQ=18)                                                 tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                       >  2:1145865/3‑101 (MQ=18)                                                 tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                       >  1:377496/3‑101 (MQ=18)                                                 tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                       >  1:3257142/3‑101 (MQ=18)                                                 tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                       >  1:2213172/3‑101 (MQ=18)                                                 ttccGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                       <  2:915557/98‑1 (MQ=18)                                                  ttcGTTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                      <  2:3691824/100‑1 (MQ=18)                                                  ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                      >  2:2606604/2‑101 (MQ=18)                                                  ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                      <  2:2927630/100‑1 (MQ=18)                                                  ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                      <  1:486704/100‑1 (MQ=18)                                                  ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                      <  2:1397427/100‑1 (MQ=18)                                                  tcccgcccGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                      <  1:1196765/95‑1 (MQ=17)                                                   tcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAcc                                     >  2:4229831/1‑101 (MQ=18)                                                    cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                    <  1:1391811/101‑1 (MQ=18)                                                    cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                    <  2:2631301/101‑1 (MQ=18)                                                    cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                    >  2:2407166/1‑101 (MQ=18)                                                     ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt                                   >  1:2771801/1‑101 (MQ=18)                                                     ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt                                   <  2:4080705/101‑1 (MQ=18)                                                      gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg                                  <  1:1981218/101‑1 (MQ=18)                                                      gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg                                  <  1:3169218/101‑1 (MQ=18)                                                       tCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTgg                                 >  2:3207079/1‑101 (MQ=18)                                                       tCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTgg                                 <  2:2949087/101‑1 (MQ=18)                                                        ccGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTggg                                <  2:4068109/101‑1 (MQ=18)                                                         cGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGc                               <  1:3056557/101‑1 (MQ=18)                                                          gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc                              <  2:3485682/101‑1 (MQ=18)                                                          gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc                              <  1:4232886/101‑1 (MQ=18)                                                          gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc                              <  2:2293679/101‑1 (MQ=18)                                                           cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa                             >  2:1647498/1‑101 (MQ=18)                                                           cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa                             <  1:709541/101‑1 (MQ=18)                                                             cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCCCCTTGGGCCaaa                           >  1:3437281/1‑101 (MQ=17)                                                             cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaaa                           >  1:3122483/1‑101 (MQ=17)                                                             cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaaa                           <  2:3881869/101‑1 (MQ=18)                                                             cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaaa                           <  1:2508181/101‑1 (MQ=18)                                                                                cccgc..cccTGCGGGCCCCCCCCCGGGATATTTCCGCAAAGGGGAAATGACCCCGGGGCGCGGCTTTCCCCTTGGGCCAAATATCCCGGGGGGGGGGAACtt      <  1:2300278/97‑1 (MQ=21)                                                                                    c..cccTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCCCCTTGGGCCAAATAGCCCGGGGGGGGGGAACTTCCCg  <  1:2392607/101‑1 (MQ=255)                                                                                     ||                                                                                                     TGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGG....CTG...TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCAAATATCCCGGGGGGGGGGAACTTCCCG  >  NC_014034/1255828‑1256004 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG