Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,255,870	C→T	100%	intergenic (+40/+202)	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,870	0	C→T	100.0%	200.1	82	intergenic (+40/+202)	RCAP_rcc01184/RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (43/39):  ref base (0/0):  total (43/39)

AAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGG....CTG.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAA  >  NC_014034/1255795‑1255938                                                                                 |                                                                       aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggg                                                    >  1:4682701/1‑101 (MQ=255) aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggg                                                    >  1:1675466/1‑101 (MQ=255) aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggg                                                    >  2:2581626/1‑101 (MQ=255) aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggg                                                    >  1:2347690/1‑101 (MQ=255)  agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                   <  1:1204139/101‑1 (MQ=255)  agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                   >  2:3680313/1‑101 (MQ=255)  agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggg                                                   <  1:1897583/101‑1 (MQ=255)   gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggg                                                  >  2:4036338/1‑101 (MQ=255)   gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggg                                                  <  1:665197/101‑1 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggg                                                 <  2:822642/101‑1 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggg                                                 <  1:1894901/101‑1 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAggggggg                                                 >  2:3895144/1‑101 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                <  1:3747550/101‑1 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                <  2:1943650/101‑1 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                <  1:152340/101‑1 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgc                                                >  2:3484863/1‑101 (MQ=255)      ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcg                                               >  2:1274303/1‑101 (MQ=255)      ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCAGCCCGAAGAggggggggg                                               >  1:3593672/1‑99 (MQ=255)       gaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgggc                                              >  1:788473/1‑101 (MQ=255)       gaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcgc                                              >  2:1406797/1‑101 (MQ=255)       gaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcgc                                              >  1:3373694/1‑101 (MQ=255)       gaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGgcgc                                              >  1:4406057/1‑101 (MQ=255)        aagTCGACGACTTCCTTTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCt                                             <  2:2463309/101‑1 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCt                                             <  2:385634/101‑1 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCt                                             >  2:2558179/1‑101 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCt                                             <  2:714756/101‑1 (MQ=255)         agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAgggggggggggg                                            >  2:3675583/1‑96 (MQ=255)         agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtg                                            >  1:2630335/1‑101 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                         >  2:3604684/1‑98 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                         <  1:2234871/101‑4 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                         >  1:1457253/1‑98 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                         >  1:3409880/1‑98 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                         <  2:1696261/101‑4 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                         >  2:1020029/1‑98 (MQ=255)             gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                        <  2:1559349/101‑5 (MQ=255)               cgccTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                      <  1:4036580/98‑7 (MQ=255)               cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                      >  1:205501/1‑95 (MQ=255)               cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                      <  2:4332953/101‑7 (MQ=255)                gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                     <  1:2176029/101‑8 (MQ=255)                gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACAGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                     >  2:2489937/1‑94 (MQ=255)                 acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                    <  1:170926/101‑9 (MQ=255)                   ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                  <  1:743724/101‑11 (MQ=255)                   ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                  >  1:3745582/1‑91 (MQ=255)                                  tGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTccccc                   >  1:3462172/1‑101 (MQ=9)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  <  2:387019/101‑1 (MQ=9)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  >  1:1447001/1‑101 (MQ=9)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  >  2:387091/1‑101 (MQ=9)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  <  2:1498515/101‑1 (MQ=9)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  >  1:3903537/1‑101 (MQ=9)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTGCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  <  1:4123148/101‑1 (MQ=2)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTGCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTcccccc                  <  1:4124795/101‑1 (MQ=2)                                    gccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCg                 <  2:1449457/100‑1 (MQ=2)                                    aCCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCg                 <  1:986809/101‑1 (MQ=9)                                     ccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgg                <  1:2720892/101‑1 (MQ=2)                                     ccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgg                <  2:3806826/101‑1 (MQ=2)                                     ccTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgg                >  2:4419392/1‑101 (MQ=2)                                      cTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCggg               >  1:4723324/1‑101 (MQ=2)                                      cTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCggg               >  1:3181947/1‑101 (MQ=2)                                      cTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCggg               >  2:1248739/1‑101 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCgggc              >  1:3603585/1‑100 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:2827938/1‑101 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:2088336/1‑101 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:4200302/1‑101 (MQ=2)                                        aGGGGTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGat             <  2:2112014/101‑1 (MQ=2)                                         ggggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGata            <  2:911522/101‑1 (MQ=2)                                         ggggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGata            >  2:1912569/1‑101 (MQ=2)                                          gggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGatat           <  2:151761/101‑1 (MQ=2)                                           ggTGCCGGACCGGACCCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt          <  2:4333697/101‑1 (MQ=0)                                           ggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt          >  1:1093040/1‑101 (MQ=2)                                           ggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt          <  1:2087556/101‑1 (MQ=2)                                           ggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt          <  2:3940520/101‑1 (MQ=2)                                           ggTGCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATAtt          <  2:306679/101‑1 (MQ=2)                                              gCCGGACCGGAACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTcc       >  2:3026578/1‑101 (MQ=2)                                                ctgaccggaACCGGGGGCCTTCGGTCCCCTGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     <  2:1361752/99‑1 (MQ=0)                                                cggaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     >  2:157267/1‑101 (MQ=2)                                                cggaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     <  2:743454/101‑1 (MQ=2)                                                 ggaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc    <  2:3720944/101‑1 (MQ=2)                                                 ggaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc    <  2:4588845/101‑1 (MQ=2)                                                 ggaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc    >  1:2050208/1‑101 (MQ=2)                                                  gaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCa   >  1:2334832/1‑101 (MQ=2)                                                  gaccggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGACATTTCCGCCa   >  1:1202385/1‑101 (MQ=2)                                                   accggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  1:4383608/101‑1 (MQ=2)                                                   accggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  >  1:212708/1‑101 (MQ=2)                                                   accggaACCGGGGGCCTTCGGGCCCCCGCTTTCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  2:539925/101‑1 (MQ=2)                                                                                 |                                                                       AAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCGG....CTG.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAA  >  NC_014034/1255795‑1255938 Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG