Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	369,202	A→G	I198V (ATT→GTT)	mhpA →	3‑(3‑hydroxyphenyl)propanoate hydroxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	369,202	0	A	G	100.0%	240.3 / NA	73	I198V (ATT→GTT)	mhpA	3‑(3‑hydroxyphenyl)propanoate hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/42);  total (31/42)

CAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGATTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGTTGAATTGATTCGCCAGCGTGTCTACACCCACAACGCGC  >  NC_000913/368962‑369448                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       cAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGGAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGACTGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGtcgg                                                                                                                                                                                                                                              >  2:19052/1‑248 (MQ=255)  aGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGt                                                                                                                                                                                                                                                 >  1:109583/1‑247 (MQ=255)      cGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCg                                                                                                                                                                                                                                         <  1:437975/251‑1 (MQ=255)          gCGGTGCTGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGGAAAACTGCGCCAAATCAGTGGGTTGTGGGAGATATCGCCa                                                                                                                                                                                                                                      >  2:63681/1‑250 (MQ=255)           cGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGGTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCCGAATGTGCCGTTGGAGGGTAAAACTGCGCTAAATCAGTGGGTTGTGGTAGATATCGCCAAc                                                                                                                                                                                                                                    <  1:26336/251‑1 (MQ=255)            ggTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTTCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACg                                                                                                                                                                                                                                   <  2:511078/251‑1 (MQ=255)                     ggAAGGGTTGTCGCTTGTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAATATGACGAAGTGACCTTGCGCCTGAAAACGGCATAAGGGCAGCGGGAAATAGGCAAAGCCCAGTGGCTGGTAGCCTTTGACGGTGGAGCATGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAATCTGCGCCAAATCAGTGGGTTGTGTTAGATATCGCCAACGATCCGtt                                                                                                                                                                                                                            <  2:470894/249‑1 (MQ=255)                              gtCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTGTTGTCCGTCGCACTCTGAATGTGCCGTATGAAGGTAAAACTGCGCCAAATCAGTGGGGTGTGGTAGATATCGCCAACGATCCGTTAAGTAcgccg                                                                                                                                                                                                                  <  2:280381/250‑1 (MQ=255)                               tCGCGTTTTCCGAACGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCACAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAATAcgccgcg                                                                                                                                                                                                                >  2:88214/1‑250 (MQ=255)                                cgcgTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCat                                                                                                                                                                                                               >  1:134168/1‑251 (MQ=255)                                cgcgTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCat                                                                                                                                                                                                               <  2:38934/251‑1 (MQ=255)                                  cgTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGGAGATATCGCCAACGATCCGTTAAGTACGCCGCatat                                                                                                                                                                                                             >  2:521358/1‑251 (MQ=255)                                   gTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATc                                                                                                                                                                                                            <  1:421349/251‑1 (MQ=255)                                        ccGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCCGTGGGTTGTGGTAGATATCGCCAACGAACCGTTTAGTGCGCCGCATTTCTAttt                                                                                                                                                                                                       >  2:22328/1‑251 (MQ=255)                                        ccGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTAttt                                                                                                                                                                                                       <  1:483815/251‑1 (MQ=255)                                        ccGAATGTGCGCTGCTTCTTTTCCCGCGAGCGGGAGGCCTACATGCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGTCAGCGGGAAATAGTCAAAGCCCAGTCGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAGCTGCGCCAAATCAGTGGGTTGTGATAGATATCGCCAACGATCCGGGAAGTACGCCGCATATCTAttt                                                                                                                                                                                                       <  2:77950/251‑1 (MQ=255)                                           aaTGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTtgtg                                                                                                                                                                                                    <  1:365311/251‑1 (MQ=255)                                              gtgCGCTGCTTGTTTTCCCCAGAGTTGGCGGACTTCAGTCCGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAACTCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGCACGCCGCATATCTATTTGTGTTg                                                                                                                                                                                                 >  2:335039/1‑251 (MQ=255)                                                gcgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTAAGCCAGATGAAGCAGTGACCTTGCACCTGAAAACGAAAGAAGGACAGCGGGAAATAGTAAAAGCCCAGTGACTGGTAGCCTGTGACGGCGGAGCACGTTTTGTCCGTCGCACTCTGAACGTGCCGTTTAACGCTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATACGTTAAGTACGCCGCATATCTATTTGTGTTGCg                                                                                                                                                                                               >  1:505944/1‑251 (MQ=255)                                                 cgcTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGCGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGa                                                                                                                                                                                              <  2:379778/251‑1 (MQ=255)                                                       ttGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGt                                                                                                                                                                                        >  2:117059/1‑251 (MQ=255)                                                         gTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTgc                                                                                                                                                                                      <  1:22328/251‑1 (MQ=255)                                                                   gAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGGGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGt                                                                                                                                                                            <  1:347081/251‑1 (MQ=255)                                                                                tcagtcagCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTGTCTGCCGCGCt                                                                                                                                                                <  2:70423/250‑1 (MQ=255)                                                                                 cagtcacCAAGATGACGAAGTGACCTTGCACCTGAACACGGCAGACGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGAAGCCCGTGACGGCGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCTTATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGc                                                                                                                                                              >  1:240634/1‑251 (MQ=255)                                                                                  agtcagCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTTTTTGTGTTGCGATCCGGTGCGCCCGGATGTTTCTGCCGCGGTGcc                                                                                                                                                             >  2:296295/1‑251 (MQ=255)                                                                                     cagCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTc                                                                                                                                                           <  2:188902/250‑1 (MQ=255)                                                                                      atcAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTTTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCAt                                                                                                                                                         <  1:534354/249‑1 (MQ=255)                                                                                      agCAAGAGGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGGCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAAGCAGTGGGTTGTGGTAGATAGCGCCAACGATCCGGTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCAt                                                                                                                                                         <  2:458218/251‑1 (MQ=255)                                                                                                 gAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCg                                                                                                                                              <  2:387499/251‑1 (MQ=255)                                                                                                 gAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGCCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGATTTGTCCGTCGCACTCTGAATGTGACGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATAGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCg                                                                                                                                              >  1:76806/1‑251 (MQ=255)                                                                                                 gAAGTGACCTTGCACCTGAAAAAGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTAAGCCGCCTATCTATTTGTGTTGCGATCCGGGGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCg                                                                                                                                              >  2:250380/1‑251 (MQ=255)                                                                                                  aaGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGc                                                                                                                                             <  1:176850/251‑1 (MQ=255)                                                                                                               ccTGAAAACGGCATAAGGGCAGCGGGAAATAGGCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCTTTTGAAGGTAAAACTGCGCCAAAGCAGTGGGTTGTGGTAGATATCGCCAACTATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCTCTGCCTCATGCGGTACGTCTCTTTTAATTTATgg                                                                                                                                <  2:478087/251‑1 (MQ=255)                                                                                                                   aaaaCGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTGCGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCGGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGAt                                                                                                                            <  1:341592/251‑1 (MQ=255)                                                                                                                           aGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCggg                                                                                                                      >  1:117748/1‑249 (MQ=255)                                                                                                                            gAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCGCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGgaga                                                                                                                   <  1:244965/251‑1 (MQ=255)                                                                                                                               gggCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAAc                                                                                                                <  2:379878/251‑1 (MQ=255)                                                                                                                                   aGCGGGAAATCGTAAAGGACCAGTGGCTGGTAGCCTGTGACGATGGAGGAAGTTATATACGTAGCCATCTAACTATCCAGATTGCAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCACCGATACGTTAAGTACGCCGCATATATATTTGTGTTGCGATCCGGTGCGCCCGCATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGTTGCCGGgagacaacgc                                                                                                             >  2:506977/1‑244 (MQ=255)                                                                                                                                           aTAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCt                                                                                                    <  1:519857/251‑1 (MQ=255)                                                                                                                                                cAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTg                                                                                               <  2:413689/251‑1 (MQ=255)                                                                                                                                                  aaGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCCAGTTTTGTCCGTCGCACTGTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAACTCAGTGGGTTGTGGTAGATATCGCCAACGAACCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTTTGGTGATGCggggcga                                                                                                                   >  2:519777/1‑226 (MQ=255)                                                                                                                                                  aaGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGgaga                                                                                                                   <  1:519777/229‑1 (MQ=255)                                                                                                                                                  aaGCCCAGTGGCTGGTAGCATGTGACGGTGGAGCACGATCTGTCCGTCGCACTCTGCATGTGCCGTTTGAAGGAAAAACTGAGACAACTCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTg                                                                                               >  2:323742/1‑249 (MQ=255)                                                                                                                                                      ccAGTGGCTGGTAGCCTGTGACGGTGGAGCCCGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGc                                                                                         >  1:103377/1‑251 (MQ=255)                                                                                                                                                      ccAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCGGCCTCAGGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGc                                                                                         <  1:323742/251‑1 (MQ=255)                                                                                                                                                         gTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCaaa                                                                                      <  1:326839/251‑1 (MQ=255)                                                                                                                                                            gCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCAGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGAGGGTTGTGGTAGATCTCACCAACGATCCGTTAAGTACGACGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGCGATGCCGGGAGAAACCGACGAGCAGCTGCGTGAGCCGCAAAat                                                                                    >  1:257615/1‑250 (MQ=255)                                                                                                                                                              tGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATg                                                                                 >  1:382206/1‑251 (MQ=255)                                                                                                                                                                         gACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGtt                                                                      <  1:296295/251‑1 (MQ=255)                                                                                                                                                                         gACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTAGGGTGAGGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGtt                                                                      <  2:257615/251‑1 (MQ=255)                                                                                                                                                                           cGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTaa                                                                    <  1:63681/251‑1 (MQ=255)                                                                                                                                                                                gAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGCAGGTACAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGAGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGCGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCaa                                                                >  1:267369/1‑250 (MQ=255)                                                                                                                                                                                    aaGTTTTGTCCGTCGCACTCTGAATGTGACGTTTGACGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGCAACCGAAGAGCAGCTGCGTGCGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGtgc                                                           >  2:210401/1‑251 (MQ=255)                                                                                                                                                                                         ttGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGCGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGcc                                                       <  1:27919/250‑1 (MQ=255)                                                                                                                                                                                             ccGTCGCACTCTGAATGTGCCGTTTGAAGGTATAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGACGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGATTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATc                                                  >  2:72078/1‑251 (MQ=255)                                                                                                                                                                                               gTCGCACTCTGAATGGGCCGTTGGAAGGGAAAACTGCGGCAAATCAGTGGGTTGTGGGAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGAGCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGAGGGCGGGAGAAACCGAAGAGCAGCTGCGGGAGCCGCAAAATATGCGCAATCTGTTAATCAAAGTGCTGCCTAATCCg                                                <  2:324915/251‑1 (MQ=255)                                                                                                                                                                                                tCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCgg                                               >  2:319229/1‑251 (MQ=255)                                                                                                                                                                                                  gCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGAc                                             <  1:117059/251‑1 (MQ=255)                                                                                                                                                                                                   cACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGGAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCTCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAAGATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACa                                            <  1:335039/251‑1 (MQ=255)                                                                                                                                                                                                        ggAATGTGCCGTGTGAAGGTAAAACTGCGCCAAATCAGGGGGTTGTGGTAGATAGCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTGCGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTGATCCGGACAATGtt                                       <  2:65796/250‑1 (MQ=255)                                                                                                                                                                                                         gAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCTGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATAGGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGt                                        <  1:250380/249‑1 (MQ=255)                                                                                                                                                                                                               gCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGGTTGTGGTAGATAGCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTAGGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTTTTAAGCAAAGTGCGGCCTAAACCTGACAATGTTGAATTg                                 <  2:375506/250‑1 (MQ=255)                                                                                                                                                                                                                              aaCTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGTTGAATTGATTCGCCAGCGTGTCt                 >  1:49939/1‑251 (MQ=255)                                                                                                                                                                                                                               aCTGCGCCAAATCAGTGGGTTGTGGTAGATTTCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTGGCGAGCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCGGTTAAGCAAAGTGCTGCCTAAGCTGGACAATGGTGAGTTGATTCGCCAGCGTGTCTa                <  2:267369/251‑1 (MQ=255)                                                                                                                                                                                                                               aCTGCGCCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGTTGAATTGATTCGCCAGCGTGTCTa                <  1:88214/251‑1 (MQ=255)                                                                                                                                                                                                                                 tGCGCCAAATCAGTGGGTTGTGGTCGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGTTGAATTGATTCGcc                          >  1:438660/1‑239 (MQ=255)                                                                                                                                                                                                                                 tGCGCCAAATCAGTGGGTGGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCGTTGAATTTAGGGTGAGGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAAGCCGGACAATGTTGAATGGATTCGcc                          <  2:438660/239‑1 (MQ=255)                                                                                                                                                                                                                                   cgcCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATACGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGCGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGCCAATGTTGAATTGATTCGCCAGCGTGTCTACAcc            >  1:201254/1‑251 (MQ=255)                                                                                                                                                                                                                                    gcCAAATCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGAATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCACGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGTTGAATTGATTCGCCAGCGTGTCTACAccc           >  2:194994/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCg                                                                                                                        <  1:127165/137‑1 (MQ=255)                                                                                                                                                                                                                                         aTCAGTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCg                                                                                                                        >  2:127165/1‑137 (MQ=255)                                                                                                                                                                                                                                             gTGGGTTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGAAAATGTTGAATTGATTCGCCAGCGTGTCTACACCCCCAAcgcgc  >  2:231283/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CAGGTCGATGCGGTGATGCTGGAAGGGGTGTCGCGTTTTCCGAATGTGCGCTGCTTGTTTTCCCGCGAGCTGGAGGCCTTCAGTCAGCAAGATGACGAAGTGACCTTGCACCTGAAAACGGCAGAAGGGCAGCGGGAAATAGTCAAAGCCCAGTGGCTGGTAGCCTGTGACGGTGGAGCAAGTTTTGTCCGTCGCACTCTGAATGTGCCGTTTGAAGGTAAAACTGCGCCAAATCAGTGGATTGTGGTAGATATCGCCAACGATCCGTTAAGTACGCCGCATATCTATTTGTGTTGCGATCCGGTGCGCCCGTATGTTTCTGCCGCGCTGCCTCATGCGGTACGTCGCTTTGAATTTATGGTGATGCCGGGAGAAACCGAAGAGCAGCTGCGTGAGCCGCAAAATATGCGCAAGCTGTTAAGCAAAGTGCTGCCTAATCCGGACAATGTTGAATTGATTCGCCAGCGTGTCTACACCCACAACGCGC  >  NC_000913/368962‑369448

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