Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,150,316	G→A	E219K (GAG→AAG)	yghA →	NADP(+)‑dependent aldehyde reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,150,316	0	G	A	100.0%	116.0 / NA	37	E219K (GAG→AAG)	yghA	NADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/17);  total (20/17)

GGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTT  >  NC_000913/3150071‑3150560                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ggAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCgaagaa                                                                                                                                                                                                                                                  >  2:330235/1‑250 (MQ=255)        ggTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTc                                                                                                                                                                                                                                          <  2:99259/251‑1 (MQ=255)        ggTTGCCATTACGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCAGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTc                                                                                                                                                                                                                                          >  1:294035/1‑251 (MQ=255)           tGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCgg                                                                                                                                                                                                                                        <  1:215686/250‑1 (MQ=255)                                         cGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGa                                                                                                                                                                                                         <  1:220445/251‑1 (MQ=255)                                                tttCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACt                                                                                                                                                                                                  >  2:299692/1‑251 (MQ=255)                                                     aaaGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCaa                                                                                                                                                                                              >  2:315811/1‑250 (MQ=255)                                                              tGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcggcg                                                                                                                                                                                    >  1:176503/1‑251 (MQ=255)                                                                            tCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGAt                                                                                                                                                                      >  2:124271/1‑251 (MQ=255)                                                                             cgcgCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCTCTTCGTCAAGCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCATGATa                                                                                                                                                                     <  2:193612/251‑1 (MQ=255)                                                                                gcTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTTCAAGTATCATCACCACTTCGTCAATCCAGGCATTCCAGCCAAGTCCGCATTTACTGTACTATGCGGCTTCGAAGGCGGCGATACTGAAGTACAGCCGTGGCTTGTCAAAGCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGa                                                                                                                                                                  <  1:397012/251‑1 (MQ=255)                                                                                 cTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCCGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATa                                                                                                                                                                     >  2:349214/1‑247 (MQ=255)                                                                                     tCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGGCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATccc                                                                                                                                                              <  2:311751/250‑1 (MQ=255)                                                                                      cTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGCCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAACCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGc                                                                                                                                                            >  1:282183/1‑251 (MQ=255)                                                                                           tAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACg                                                                                                                                                                            >  2:359381/1‑230 (MQ=255)                                                                                           tAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACg                                                                                                                                                                            <  1:359381/230‑1 (MQ=255)                                                                                           tAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGttt                                                                                                                                                       >  2:194323/1‑251 (MQ=255)                                                                                                   gAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGc                                                                                                                                                >  1:304213/1‑250 (MQ=255)                                                                                                                     ccGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGGGAATATTGTCGCGCCAGGCCCGATCGGGACAGCACTGCAAATTTCGGGCGAACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGc                                                                                                                             <  2:176503/251‑1 (MQ=255)                                                                                                                                    atcatcACCTCGTCGTCAATCCAGGCATACCAGCCATGTCCGCATGTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAACACAGGTCGCGAAGAAAGGTAGTCGGGTGAATAGTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCTAAcggcgg                                                                                                               <  2:100950/250‑1 (MQ=255)                                                                                                                                       cgccccccTTCGTCAAGCCAGGCAAACCAGCCAAGTCCGAATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTGCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGAGGAACt                                                                                                           <  2:45970/244‑1 (MQ=255)                                                                                                                                         caccacTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTgg                                                                                                         <  1:368496/251‑1 (MQ=255)                                                                                                                                                tCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCt                                                                                                    <  1:349214/249‑1 (MQ=255)                                                                                                                                                     aaTCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATg                                                                                             >  1:107282/1‑251 (MQ=255)                                                                                                                                                      aTCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGt                                                                                            >  1:80890/1‑251 (MQ=255)                                                                                                                                                         cAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGttt                                                                                          >  1:247560/1‑250 (MQ=255)                                                                                                                                                                     ccAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTc                                                                              >  1:184728/1‑250 (MQ=255)                                                                                                                                                                              cATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAgacg                                                                      >  1:108815/1‑246 (MQ=255)                                                                                                                                                                                 ttACTGGACTATGCGGCTACGAAGGCGGCTATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCt                                                                  <  2:12897/250‑1 (MQ=255)                                                                                                                                                                                           aggCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTGTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCg                                                         <  2:304213/247‑1 (MQ=255)                                                                                                                                                                                                  tACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGa                                                                          <  1:94947/225‑1 (MQ=255)                                                                                                                                                                                                  tACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGa                                                                          >  2:94947/1‑225 (MQ=255)                                                                                                                                                                                                                    gAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAAGTATTCGGGTGAATATTGTCGCGCCAGGCACGATCTGGACAGAACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATACCGAAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCACTGTATATGTTTATCTGGCAAGACAGGAATCGAGCTACGTCACCGAAGATGTTGACGGCGTGTgcggcagcg                              >  1:1277/1‑251 (MQ=255)                                                                                                                                                                                                                                 ggCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGtt                  <  1:299692/250‑1 (MQ=255)                                                                                                                                                                                                                                    ttGGCAAAACAGGTCGCGAAGCAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTaaa               >  2:339495/1‑250 (MQ=255)                                                                                                                                                                                                                                             cAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCg     <  2:247560/251‑1 (MQ=255)                                                                                                                                                                                                                                                gTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGtt  >  1:147037/1‑251 (MQ=255)                                                                                                                                                                                                                                                gTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATATGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAAATGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTAGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGGGAGCATTTAGGTTAAAAAAATGACCGGtt  >  1:148/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     GGAAACAGGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGTT  >  NC_000913/3150071‑3150560

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