Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	499,060	T→G	D305A (GAC→GCC)	aes ←	acetylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	499,060	0	T	G	100.0%	214.8 / NA	61	D305A (GAC→GCC)	aes	acetylesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/31);  total (30/31)

CGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATG  >  NC_000913/498816‑499294                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCggcggc                                                                                                                                                                                                                                        <  2:18831/249‑1 (MQ=255)       aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGCGGTCTTCCTCGGTGCCGGCGGGCAAAAATATGAATATCTTTCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCCCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGGCGCGGAGAGCCTCGGCCGCGGGTTtcc                                                                                                                                                                                                                                >  1:716195/1‑250 (MQ=255)        aaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGGGGGTTtcat                                                                                                                                                                                                                               >  1:481379/1‑251 (MQ=255)                    cATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagt                                                                                                                                                                                                                   <  1:545811‑M1/251‑5 (MQ=255)                          cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtag                                                                                                                                                                                                              <  2:231245‑M1/250‑10 (MQ=255)                          cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtag                                                                                                                                                                                                              <  2:481379‑M1/250‑10 (MQ=255)                             gTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGGGCTTTATGCCCCGCCGGGACATATACAAAAGATGGCTTATCTTTTTTCCCCGTATCGCTACAGCTGAGCGGGAAAGAAATGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGgttttactattcgggtaggccgagggc                                                                                                                                                                                                           >  2:728546/1‑227 (MQ=255)                               gATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAAATGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTTATCATCCGTGtagtcgtaggcattt                                                                                                                                                                                                        >  1:612699‑M1/1‑236 (MQ=255)                                 tgtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAATATATTCCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGGCGCGGATCGATAAAGATGAGCGGGAAAGAAATGAGCGGCGGCGCGAAGAGGCTCGGCGGGGGGTTTCATCATCCGGGtaggcggaggccat                                                                                                                                                                                                         >  1:68629‑M1/1‑234 (MQ=255)                                             ttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCGTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatgtgtgcct                                                                                                                                                                                           <  2:663786‑M1/250‑29 (MQ=255)                                                   cGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcga                                                                                                                                                                                    <  1:31816‑M1/251‑36 (MQ=255)                                                               ggacGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGACGGCGGGAAAAAACAGGAATAGATTCCGGCGCTTAAGGCCACGCCGGAACAGAGCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgg                                                                                                                                                                        <  2:569554‑M1/250‑48 (MQ=255)                                                                  cGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtctt                                                                                                                                                                     <  2:631168‑M1/251‑51 (MQ=255)                                                                         agagATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcggaggcattagacatttgtgcctgatgcgacgctttgcgcgtcttatcaggc                                                                                                                                                              >  1:529246‑M1/1‑194 (MQ=255)                                                                            gaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggccta                                                                                                                                                           <  1:583019‑M1/251‑61 (MQ=255)                                                                                               gAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccg                                                                                                                                         <  1:60670‑M1/250‑79 (MQ=255)                                                                                                       ccTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggt                                                                                                                                <  1:135563‑M1/251‑88 (MQ=255)                                                                                                                cggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGGGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcctccggcaattggggcacaatgc                                                                                                                       >  1:11001‑M1/1‑155 (MQ=255)                                                                                                                 ggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcc                                                                                                                      <  1:707760‑M1/251‑98 (MQ=255)                                                                                                                  gcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcct                                                                                                                     <  1:594521‑M1/251‑99 (MQ=255)                                                                                                                   cggGAAAAAATATGAATATATGCCGGCGCCTAATGCCACGCCGGAACATATCGAAATGATGGCGAATCGTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCAGCATCCGTGtaggcgtaggcattagacagttgtgcctgatgcgaggcttgacgcgtcttatcaggcctacaaccggatccgcatccggcaattggtgttcaaagcctg                                                                                                                    <  2:652676‑M1/251‑100 (MQ=255)                                                                                                                          aaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgat                                                                                                             >  2:83434‑M1/1‑145 (MQ=255)                                                                                                                              atGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgctt                                                                                                         <  2:252159‑M1/251‑111 (MQ=255)                                                                                                                                    atatTCCGGCGCTTAATGCAACGCAGGCACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgc                                                                                                    >  1:535742‑M1/1‑135 (MQ=255)                                                                                                                                          cGGCGATTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttat                                                                                             <  2:488461‑M1/251‑123 (MQ=255)                                                                                                                                           ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTACTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGACGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattcgacatttgtgcctgctgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatc                                                                                            >  1:694763‑M1/1‑128 (MQ=255)                                                                                                                                           ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatccggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgctgctatgcttgacgcatcctctc                                                                                            >  2:10365‑M1/1‑128 (MQ=255)                                                                                                                                           ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGcagtcgtaggcattagacatttgtgcctgatgcga                                                                                                                                                                                    >  1:677164‑M1/1‑128 (MQ=255)                                                                                                                                           ggCGCTTAATGCCACGCCGGAACATAGCGAAATGATGGCTAATCTTGTGGCCGCGTATCGCTAAAGCTTAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacattggtgcctgatgcga                                                                                                                                                                                    <  2:677164‑M1/163‑36 (MQ=255)                                                                                                                                             cgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATATTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatc                                                                                            >  1:696500‑M1/1‑126 (MQ=255)                                                                                                                                                  aaTGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctaccaccggtgacgcatccggcaattggtgcccaatgcctgatgcgatgcttgacgcatcttatcaggcct                                                                                      >  1:130697‑M1/1‑121 (MQ=255)                                                                                                                                                   ctgCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctac                                                                                    >  1:20036‑M1/2‑120 (MQ=255)                                                                                                                                                      ccACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaat                                                                                 <  1:250859‑M1/251‑135 (MQ=255)                                                                                                                                                      ccACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaat                                                                                 >  1:345639‑M1/1‑117 (MQ=255)                                                                                                                                                           cctgaacATAACGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggta                                                                            <  2:612699‑M1/248‑140 (MQ=255)                                                                                                                                                            cGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgca                                                                                                                             >  1:617074‑M1/1‑111 (MQ=255)                                                                                                                                                            cGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgca                                                                                                                             <  2:617074‑M1/201‑91 (MQ=255)                                                                                                                                                             ggAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatccggcctacaatgggtacc                                                                          >  2:23078‑M1/1‑110 (MQ=255)                                                                                                                                                               aaCATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGCAGAGCCTCGGCGGCGGTTTTCATCCTCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatacggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccgg                                                                        >  1:719454‑M1/1‑108 (MQ=255)                                                                                                                                                                 cATATCGAAATGAGGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggat                                                                      <  2:8245‑M1/251‑146 (MQ=255)                                                                                                                                                                  atatCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatc                                                                     <  1:701193‑M1/251‑147 (MQ=255)                                                                                                                                                                      cGAAAGGATGGATAAGCTTGGTGCCGCGTATAGCTAAAGCTGCGCGGTAAAGAACGGAGCGCCGGCGCGAAGCGCCTCGGCGGCGGTTTGCATCAGCCGTGtagtcgtaggcattagacatttgtgcctgatgcggcgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatggggaccggatcggt                                                                  <  2:298028‑M1/250‑150 (MQ=255)                                                                                                                                                                        aaaTGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGAAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtagg                                                               >  1:327100‑M1/1‑99 (MQ=255)                                                                                                                                                                          atgatgGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatcctgcaattggtgcacaatgcctgatgctatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggcc                                                             <  1:470850‑M1/251‑155 (MQ=255)                                                                                                                                                                          atgatgGCTAATCTTGGTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggcc                                                             <  2:20036‑M1/251‑155 (MQ=255)                                                                                                                                                                              tgGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCCTCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggac                                                         >  2:471326‑M1/1‑93 (MQ=255)                                                                                                                                                                                        tgttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtagtcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgaggcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgttt                                                <  1:479245‑M1/250‑168 (MQ=255)                                                                                                                                                                                         gttgCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGGCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcggcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttagcaggcctacaatgggtaccggatcgggaggccggataaggcgtttac                                              <  1:264436‑M1/251‑170 (MQ=255)                                                                                                                                                                                               gcgTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgc                                         >  1:385045‑M1/1‑76 (MQ=255)                                                                                                                                                                                                 gTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatc                                      >  1:40861‑M1/1‑74 (MQ=255)                                                                                                                                                                                                   aTCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttat                                                                                             <  2:188557‑M1/194‑123 (MQ=255)                                                                                                                                                                                                   aTCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttat                                                                                             >  1:188557‑M1/1‑72 (MQ=255)                                                                                                                                                                                                         aaaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGCGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattcgccatttgtgcatgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgctgcttgaagcatcttatcaggcctacaatgggtgccggatcggtaggacggctaaggcgtttacgccgcatccggcacg                               >  1:97756‑M1/1‑66 (MQ=255)                                                                                                                                                                                                         aaaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaa                                <  2:535742‑M1/249‑184 (MQ=255)                                                                                                                                                                                                         aaaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgat                                                                                                                  >  1:95295‑M1/1‑66 (MQ=255)                                                                                                                                                                                                         aaaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgat                                                                                                                  <  2:95295‑M1/167‑102 (MQ=255)                                                                                                                                                                                                          aaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtcggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggt                                                                                                                                                   >  1:354184‑M1/1‑65 (MQ=255)                                                                                                                                                                                                          aaGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCAGCATCCGTGtagtcgtaggcatgagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggt                                                                                                                                                   <  2:354184‑M1/133‑69 (MQ=255)                                                                                                                                                                                                           aGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggagaaggcgtttacgccgcatccggcaagaa                             >  2:20091‑M1/1‑64 (MQ=255)                                                                                                                                                                                                                      aaaGAACTGGGCGCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagt                 <  2:11001‑M1/251‑199 (MQ=255)                                                                                                                                                                                                                                    cGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttacc   >  2:560806‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                                     gTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactc  >  1:46345‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           CGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATG  >  NC_000913/498816‑499294

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