Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,495,621	C→T	A197T (GCC→ACC)	ahr ←	NADPH‑dependent aldehyde reductase Ahr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,495,621	0	C	T	100.0%	188.6 / NA	55	A197T (GCC→ACC)	ahr	NADPH‑dependent aldehyde reductase Ahr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/32);  total (23/32)

GCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGAT  >  NC_000913/4495333‑4495888                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              gCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGc                                                                                                                                                                                                                                                                    >  1:335306/1‑300 (MQ=255) gCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCAACTCGc                                                                                                                                                                                                                                                                    >  2:282904/1‑300 (MQ=255) gCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGAAGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCAACACTTTCTCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTc                                                                                                                                                                                                                                                                      >  2:200851/1‑298 (MQ=255)   tCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCAAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCAt                                                                                                                                                                                                                                                                  >  2:730567/1‑300 (MQ=255)      tAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATccc                                                                                                                                                                                                                                                               <  1:645662/300‑1 (MQ=255)            gTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGc                                                                                                                                                                                                                                                          <  1:330621/299‑1 (MQ=255)                        gAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGtttt                                                                                                                                                                                                                                             >  1:400892/1‑300 (MQ=255)                             atagaCGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAg                                                                                                                                                                                                                                         <  1:626024/297‑1 (MQ=255)                             agagACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAg                                                                                                                                                                                                                                         <  1:625996/299‑1 (MQ=255)                                  cGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTAAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTAGCATCCCATTGCGTGCAGAAGTTTTATAGCGatat                                                                                                                                                                                                                                   >  1:90288/1‑300 (MQ=255)                                        gATCGCCCGCAATTGGCGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCGTAGGTCAGCGCATCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGTAGAAGTTTTATAGCGGTATGcccc                                                                                                                                                                                                                              <  2:483745/299‑1 (MQ=255)                                         aTCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAg                                                                                                                                                                                                                            <  2:465280/300‑1 (MQ=255)                                                gCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTTTCACTTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGcc                                                                                                                                                                                                                          <  2:85394/295‑1 (MQ=255)                                                   aTTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATa                                                                                                                                                                                                                  >  1:675338/1‑300 (MQ=255)                                                             aGGCCGGAACAGACAACGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTTACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCTCCAGCCCGCCAATACCAATTAccc                                                                                                                                                                                                        <  2:220872/300‑1 (MQ=255)                                                                               gcgTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGTTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAg                                                                                                                                                                                      <  2:647084/300‑1 (MQ=255)                                                                                     gaACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGANGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑atat                                                                                                                                                                               <  1:151465/300‑1 (MQ=255)                                                                                     gaACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑atat                                                                                                                                                                               <  1:629926/300‑1 (MQ=255)                                                                                     gaACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTAACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑atat                                                                                                                                                                               <  1:236237/300‑1 (MQ=255)                                                                                     gaACCGCACCGACCGTATGGAAATGACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAANGAGAGCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑atat                                                                                                                                                                               <  2:608435/300‑1 (MQ=255)                                                                                       aCCGCACCTACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTTACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCACCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATgg                                                                                                                                                                             <  2:23307/300‑1 (MQ=255)                                                                                       aCCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAACACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATgg                                                                                                                                                                             <  2:26181/300‑1 (MQ=255)                                                                                                accgTATGGAAATTACCGCCATAGGTCAGCGACTCAAAATAGGGCGGCCAGGCGAGGCTGACGTTGATGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGAGCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCTGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACTCGGCTGGTAGCAGTG‑ATATGGTGCATCAAc                                                                                                                                                                    <  2:193994/300‑1 (MQ=255)                                                                                                   gTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAg                                                                                                                                                                  <  2:690277/299‑1 (MQ=255)                                                                                                            ttatcaccATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCGGACGTTGACGGTGTTGATAATGAGATCAAACNGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTTTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAg                                                                                                                                                        <  1:149424/294‑1 (MQ=255)                                                                                                            ttACCGCCATAGGTCAGCGCCTCAAAATAGTGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAg                                                                                                                                                        <  1:94294/300‑1 (MQ=255)                                                                                                            ttACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACNGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAg                                                                                                                                                        <  1:149566/300‑1 (MQ=255)                                                                                                             tACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAg                                                                                                                                                        <  1:78987/299‑1 (MQ=255)                                                                                                                      aGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATAc                                                                                                                                              <  1:194783/300‑1 (MQ=255)                                                                                                                       ggTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTCGTAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATAcc                                                                                                                                             >  2:784433/1‑300 (MQ=255)                                                                                                                               ccTCAAAATAGGGCTGNCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCNTCAACAGTGGTTTAAAGANCGTGATACCGCCGcaca                                                                                                                                     >  1:667109/1‑300 (MQ=255)                                                                                                                                   aaaaTAGTGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑g                                                                                                                                <  1:247030/300‑1 (MQ=255)                                                                                                                                      aTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATGGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCgg                                                                                                                             <  1:746592/300‑1 (MQ=255)                                                                                                                                      aTAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCgg                                                                                                                             <  1:635433/300‑1 (MQ=255)                                                                                                                                                 aGTCGAGGCTGGCGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCCTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCcggcgg                                                                                                                     <  2:63411/297‑1 (MQ=255)                                                                                                                                                       ggCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCAGGCGGACTCGata                                                                                                             >  2:501710/1‑299 (MQ=255)                                                                                                                                                           gACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAat                                                                                                        >  1:198255/1‑300 (MQ=255)                                                                                                                                                           gACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAat                                                                                                        >  2:628507/1‑300 (MQ=255)                                                                                                                                                           gACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCAAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAat                                                                                                        >  1:198230/1‑300 (MQ=255)                                                                                                                                                                     gtgtTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGc                                                                                              >  1:88645/1‑300 (MQ=255)                                                                                                                                                                                 agaTCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCAcc                                                                                  <  1:330102/300‑1 (MQ=255)                                                                                                                                                                                 agaTCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTGCGCCGGATTAGAACTAAAGGTTTTCACCTCGCATCCCATTGCGTTCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCAcc                                                                                  <  2:353457/300‑1 (MQ=255)                                                                                                                                                                                  gaTCAAACTGCCCCGCCAGTGTTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGATGGTAGCAGTG‑ATCTGGTGCATCAACAGTGGTTTAACGACCGTGATACCGCCGAACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGCATCAccc                                                                                 >  2:322189/1‑300 (MQ=255)                                                                                                                                                                                                    gTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCcgcacgca                                                               >  1:140242/1‑300 (MQ=255)                                                                                                                                                                                                    gTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCA‑TGCCTCAACCCCCCTCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCACTGTGCGCCGAGCGTCTATAGCGATACTTCCCAGCACGCCAATAATGACTCGTCGAACGTGTCTGGTAAC‑GTGGCGATTGTGCATCAACCTTGGTGTAAAGTCCGCCATACCGCCTCACA‑CACGCGGCCCGGCGGCATCGATCTCAATAGTTTCTGGCAGTGGAATCACCCAtggaccgtcagcac                                                                  >  1:355561/1‑284 (MQ=255)                                                                                                                                                                                                        tttCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACtt                                                           >  2:504464/1‑300 (MQ=255)                                                                                                                                                                                                             gTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACNACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTtct                                                         >  1:171744/1‑297 (MQ=255)                                                                                                                                                                                                             gTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACNACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTtct                                                         >  1:577838/1‑297 (MQ=255)                                                                                                                                                                                                                                aTTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCAt                                   >  2:681051/1‑300 (MQ=255)                                                                                                                                                                                                                                aTTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTAGGCAAAGCCACCGCGATTCAt                                   >  2:421697/1‑300 (MQ=255)                                                                                                                                                                                                                                    accacTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCNCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAt                                <  2:246850/299‑1 (MQ=255)                                                                                                                                                                                                                                    accacTTTATCGGCACCCATCACCAGCACTTCCT‑GCTTTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTTCAGAAGTTTTATAGCGAGATGCCCCAGCCCGCCAAGACCAATTACCNCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCTTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAt                                <  2:246791/299‑1 (MQ=255)                                                                                                                                                                                                                                                  aCCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGgcaccgcac                  <  1:285129/299‑1 (MQ=255)                                                                                                                                                                                                                                                             gCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGt      <  2:524064/300‑1 (MQ=255)                                                                                                                                                                                                                                                                 ttCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGTTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGCATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCACACTCGTCGGCACCGCACCTTTCTCTCAgcttat  >  2:360810/1‑295 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                              GCTCATAAGGCGTGCCGGTAGCAGAACCAGAGACGCTGCGATCGCCCGCAATTAACGTAAAGGCCGGAACAGACAGCGGCGTGAGAACCGCACCGACCGTATGGAAATTACCGCCATAGGTCAGCGCCTCAAAATAGGGCTGCCAGTCGAGGCTGACGTTGACGGTGTTGATAATGAGATCAAACTGCCCCGCCAGTGCTTTCAGTGCCTGCGGATCGCGGCTATTCACCACTTTATCGGCACCCATCGCCAGCACTTCCT‑GCTCTTTCGCCGGATTAGAACTAAAGGCTGTCACCTCGCATCCCATTGCGTGCAGAAGTTTTATAGCGATATGCCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTG‑ATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACA‑GCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGAT  >  NC_000913/4495333‑4495888

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