Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,627,625	T→C	L37L (TTG→CTG)	ydfG →	NADP‑dependent 3‑hydroxy acid dehydrogenase, malonic semialdehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,627,625	0	T	C	100.0%	287.8 / NA	78	L37L (TTG→CTG)	ydfG	NADP‑dependent 3‑hydroxy acid dehydrogenase, malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/47);  total (31/47)

ACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGG  >  NC_000913/1627491‑1627751                                                                                                                                       |                                                                                                                               acacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCa                                                                                                                            <  2:1256603/139‑1 (MQ=255) acacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCa                                                                                                                            >  1:917309/1‑139 (MQ=255) acacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACCGGCCGTCGCCAGGAACGGCTGCa                                                                                                                            >  1:467706/1‑139 (MQ=255)  cacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAg                                                                                                                           >  1:105231/1‑139 (MQ=255)  cacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAg                                                                                                                           >  2:949963/1‑139 (MQ=255)  cacaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAg                                                                                                                           <  2:410082/139‑1 (MQ=255)   acaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAgg                                                                                                                          >  1:324303/1‑139 (MQ=255)   acaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAgg                                                                                                                          >  1:934143/1‑139 (MQ=255)   acaAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAgg                                                                                                                          >  2:189380/1‑139 (MQ=255)               aGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACg                                                                                                              <  2:324303/139‑1 (MQ=255)               aGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACg                                                                                                              <  1:949963/139‑1 (MQ=255)                 gAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGaa                                                                                                            >  2:330865/1‑139 (MQ=255)                   aCAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                          <  1:801276/139‑1 (MQ=255)                   aCAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                          <  2:590112/139‑1 (MQ=255)                     aGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTgg                                                                                                        <  1:32444/139‑1 (MQ=255)                      gTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTggg                                                                                                       >  2:305880/1‑139 (MQ=255)                        gagaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                     <  2:934143/139‑1 (MQ=255)                        gagaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                     <  1:816290/139‑1 (MQ=255)                        gagaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                     <  2:615330/139‑1 (MQ=255)                        gagaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgag                                                                                                     <  1:666248/139‑1 (MQ=255)                         agaTGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGgaga                                                                                                    >  2:748068/1‑139 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGTAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:1388989/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:1383987/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:1431776/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:1380247/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:429365/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  1:639210/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  1:305880/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACGGGGAGAt                                                                                                   <  2:8206/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACCGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:467706/139‑1 (MQ=255)                          gatgatCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGAGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  2:917309/139‑1 (MQ=255)                          gatgatCGTTTAAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  1:1222397/139‑1 (MQ=255)                              atCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATc                                                                                               >  2:3894/1‑139 (MQ=255)                              atCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATc                                                                                               <  2:105231/139‑1 (MQ=255)                                   tttAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGtat                                                                                          >  2:332009/1‑139 (MQ=255)                                                gCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTgg                                                                             <  2:1365648/139‑1 (MQ=255)                                                              tGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                          <  2:1172857/96‑1 (MQ=255)                                                              tGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                          >  1:1172857/1‑96 (MQ=255)                                                               ggTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACcgcg                                                              >  1:1361524/1‑139 (MQ=255)                                                                gTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCgcgc                                                             <  2:404074/139‑1 (MQ=255)                                                                gTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCgcgc                                                             <  2:1361524/139‑1 (MQ=255)                                                                gTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCgcgc                                                             <  1:748068/139‑1 (MQ=255)                                                                gTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCgcgc                                                             <  2:341134/139‑1 (MQ=255)                                                                  gAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                           <  1:332009/139‑1 (MQ=255)                                                                               gtcgtTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGAt                                              >  1:465420/1‑139 (MQ=255)                                                                                           aacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAg                                                                                                                 <  1:1018771/60‑1 (MQ=255)                                                                                           aacaaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAg                                                                                                                 >  2:1018771/1‑60 (MQ=255)                                                                                             caaGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTc                                >  2:1192079/1‑139 (MQ=255)                                                                                                gggCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTgg                                                                                                        <  1:370209/64‑1 (MQ=255)                                                                                                gggCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTgg                                                                                                        >  2:370209/1‑64 (MQ=255)                                                                                                gggCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   >  2:936074/1‑69 (MQ=255)                                                                                                gggCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  1:936074/69‑1 (MQ=255)                                                                                                gggCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTg                             >  1:1133941/1‑139 (MQ=255)                                                                                                gggCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTg                             <  1:189380/139‑1 (MQ=255)                                                                                                     tAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                        >  1:878201/1‑139 (MQ=255)                                                                                                     tAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                        <  2:878201/139‑1 (MQ=255)                                                                                                       aaGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                      >  2:406684/1‑139 (MQ=255)                                                                                                       aaGTTATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtg                      >  2:936359/1‑139 (MQ=255)                                                                                                          ttATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggtg                   <  1:3894/139‑1 (MQ=255)                                                                                                          ttATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggtg                   <  1:748704/139‑1 (MQ=255)                                                                                                          ttATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAgtggtg                   <  2:1133941/139‑1 (MQ=255)                                                                                                           tATCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGc                  <  1:330865/139‑1 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   <  1:391144/57‑1 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGAt                                                                                                   >  2:391144/1‑57 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGta                                                                                           <  2:798850/65‑1 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGta                                                                                           >  1:798850/1‑65 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                           <  1:718888/97‑1 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGcgccg                                                           >  2:718888/1‑97 (MQ=255)                                                                                                            aTCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCa                 >  1:735305/1‑139 (MQ=255)                                                                                                             tCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                          >  2:806214/1‑49 (MQ=255)                                                                                                             tCGCCACTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACt                                                                                                          <  1:806214/49‑1 (MQ=255)                                                                                                                  aCTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                        <  1:1293811/126‑1 (MQ=255)                                                                                                                  aCTGGCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAg                        >  2:1293811/1‑126 (MQ=255)                                                                                                                      gCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCgcgc                                                             <  1:187567/85‑1 (MQ=255)                                                                                                                      gCCGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCgcgc                                                             >  2:187567/1‑85 (MQ=255)                                                                                                                       ccGTCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATc      <  2:678366/139‑1 (MQ=255)                                                                                                                          tCGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTg   <  2:465420/139‑1 (MQ=255)                                                                                                                           cGCCAGGAACGGCTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTgg  >  2:110430/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                               ACACAAAAGCGTTGAGGAACAGTGAGATGATCGTTTTAGTAACTGGAGCAACGGCAGGTTTTGGTGAATGCATTACTCGTCGTTTTATTCAACAAGGGCATAAAGTTATCGCCACTGGCCGTCGCCAGGAACGGTTGCAGGAGTTAAAAGACGAACTGGGAGATAATCTGTATATCGCCCAACTGGACGTTCGCAACCGCGCCGCTATTGAAGAGATGCTGGCATCGCTTCCTGCCGAGTGGTGCAATATTGATATCCTGG  >  NC_000913/1627491‑1627751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑