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	98.4%	192.3 / NA	63	D305A (GAC→GCC)	aes	acetylesterase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (32/30);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTAC  >  NC_000913/498815‑499305                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ccgtagaagaAGGCGTTGGGCATATTCAGTGGATGGGCCCGGGCTTTGCTGCGTACTGTATGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCGTAATGCCACGCC‑GGAACATATCGAAAGGATGGCTAATCGTGTTGCCGCGTAGCGCTAAAGCTGAGCGGTAAAGAACGGAGC‑GCCGTCGCCAATAGCCTCGTcggcg                                                                                                                                                                                                                                                  <  2:494116/247‑1 (MQ=255)  cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCggcggc                                                                                                                                                                                                                                                   <  1:56117/249‑1 (MQ=255)     ggaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGttt                                                                                                                                                                                                                                              <  2:470502/250‑1 (MQ=255)        aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTACGGCGCTTAATACCAAGCC‑GGAACATAACGCAATGATGGCTAATCTTGTTGTCGCCTATCGCTAAAACTGAGCGGGAAAGAACTGAGCCACCG‑CGCGAAGAGCCTCGGCGGCGGTTTtc                                                                                                                                                                                                                                            >  1:413621/1‑250 (MQ=255)          aaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTAATACCACGCC‑GGGAAATATAGAAATGATGGCTTATCTTGTTTCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCTGCGGCGGGTTttatc                                                                                                                                                                                                                                         >  1:565868/1‑251 (MQ=255)               gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTGCATCATCCg                                                                                                                                                                                                                                    <  1:444636/251‑1 (MQ=255)                      atatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCTTAATGCCACGCC‑GGAAAATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAACTGAGGGGTAAAGAAATGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCCTCAACCGTGtagtc                                                                                                                                                                                                                             >  1:154256‑M1/1‑246 (MQ=255)                      atatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGAAAAAAATATGGATATTTTCCGGCGCGTTATGCCACGCC‑GGAACCAATCGAAAAGATGGCTAATCTTTGTGCCGCGTATCGCCAAAACTGAGCGGGAAAGAACAGAGC‑GCCGTCGCGGAGAGGCCCGGCGGCGGgttttttcatcagtgtggtc                                                                                                                                                                                                                             >  1:319103/1‑235 (MQ=255)                        atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgt                                                                                                                                                                                                                           <  2:592173‑M1/251‑8 (MQ=255)                           cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtag                                                                                                                                                                                                                         <  2:455277‑M1/250‑10 (MQ=255)                           cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAATATTTCGGCGCTTAATGCCCCGCCCGGAACA‑AACGAAATGATGGATTCACTTTTTTCCGCGTAACGCTAAAAACGAGCGGTGAAGCACACAGC‑GCCCTCGCGAGGAGCCAcggcggcgggcttcacctaccaggtcaacctatg                                                                                                                                                                                                                        >  2:462400/1‑226 (MQ=255)                                              gtgctgcgGATTGTCGGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcc                                                                                                                                                                                                       <  2:467579‑M1/248‑28 (MQ=255)                                                    cGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGTTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAGGCCACGCC‑GGAACATATCGAAATGATGTCTAATCTTGTTGCCGCGTAGCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCTGTTTTTATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcga                                                                                                                                                                                               <  1:310686‑M1/251‑36 (MQ=255)                                                    cGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGGTTTTATCATCCGTGtagtcggaggcattagacatttgtgcctggtgggg                                                                                                                                                                                               >  1:166265‑M1/1‑216 (MQ=255)                                                     ggATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATTATGGCTAATCTTGTTGCCGCGTAGCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtatgcattagacatttgtgcctgatgcgac                                                                                                                                                                                              <  2:12422‑M1/251‑37 (MQ=255)                                                              ggAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgc                                                                                                                                                                                     <  1:232706‑M1/251‑46 (MQ=255)                                                              ggAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgac                                                                                                                                                                                       <  1:411057‑M1/249‑44 (MQ=255)                                                                             gaTTGCCGAGCAAAACCGTGAGGTCTGCCTCGGTGCCGGCGGGAAAAAATATGAATATAGTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggccta                                                                                                                                                                      <  2:181612‑M1/251‑61 (MQ=255)                                                                                 gCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAAGCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaac                                                                                                                                                                  <  2:455654‑M1/251‑65 (MQ=255)                                                                                     aGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggt                                                                                                                                                              <  1:7517‑M1/251‑69 (MQ=255)                                                                                        aaaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTCCATCAGCCGTGtagtcgtaggcattagtcatttgttcctgatgcgacgcttgacgcgtcgtatcaggcctacaaccggtgcc                                                                                                                                                           <  2:545411‑M1/251‑72 (MQ=255)                                                                                        aaaaCCGTGAGGTCTGCCTCGGTGCCGGCGGGAAAAAATATGAATATAGGCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAGCTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttaccaggcctacaaccggtgcc                                                                                                                                                           <  1:362463‑M1/251‑72 (MQ=255)                                                                                              gTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAATATTCCGGCGCTTAATGCCACGCC‑GGAACAAATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GGCGTCGCGCAGAGGCTCGGCGGCGGTTTTCATCATCCGTGgagtcgtaggccattgacatttgtgcctgatgcgacgcttgacgcgtcttattaggcctacaacccgtgccgcctcc                                                                                                                                                     >  1:429020‑M1/1‑174 (MQ=255)                                                                                                gAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCGCCTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGGTTTTTTCATCCGTGtagtcgtaggcattagaaatttttgcctgatgcgacgcttgacgcgtcttttcaggcctacaaccggtgccgcatccgg                                                                                                                                                   >  1:85949‑M1/1‑172 (MQ=255)                                                                                                     cTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGGTGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgctttacgcgtcttatcaggcctaccaccggtgccgcatccggccatt                                                                                                                                              >  1:270062‑M1/1‑167 (MQ=255)                                                                                                       tCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAAGGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaat                                                                                                                                               <  2:192837‑M1/248‑84 (MQ=255)                                                                                                                   gcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcct                                                                                                                                >  1:502714‑M1/1‑153 (MQ=255)                                                                                                                          aaaaTATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCCAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaacctgtgccgcatccggcaattggtgcacaatgcctgatgcga                                                                                                                         <  1:340258‑M1/251‑106 (MQ=255)                                                                                                                           aaaTATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacccttgacgcgtcttatccggcctacaaccggggccgcatcctgcaattggggcacaatgcctgatgggtg                                                                                                                        >  2:207454‑M1/1‑145 (MQ=255)                                                                                                                             atatGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgc                                                                                                                      <  1:550812‑M1/251‑109 (MQ=255)                                                                                                                             atatGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgc                                                                                                                      >  1:147545‑M1/1‑143 (MQ=255)                                                                                                                                 gAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GACGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacacttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgacgcttgt                                                                                                                  >  1:379157‑M1/1‑139 (MQ=255)                                                                                                                                  aaTATATGCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgac                                                                                                                 <  1:69902‑M1/251‑114 (MQ=255)                                                                                                                                         tCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatctt                                                                                                          >  2:117702‑M1/1‑131 (MQ=255)                                                                                                                                           cGGCGCTTAATGCCACGCC‑GGAACCTATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggccttagacctttgtgcctgatgcgacgcttgacgtgtcttatcaggcctccaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttat                                                                                                        >  2:347984‑M1/1‑129 (MQ=255)                                                                                                                                              cgcTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatc                                                                                                       >  1:431249‑M1/1‑126 (MQ=255)                                                                                                                                                   aaTGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcct                                                                                                 >  1:119647‑M1/1‑121 (MQ=255)                                                                                                                                                    aTGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctac                                                                                               >  1:562677‑M1/1‑120 (MQ=255)                                                                                                                                                    aTGCCACGCC‑GGAACATATCGAAATGATGGCGAATCTGGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcggcgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctac                                                                                               <  2:501066‑M1/251‑132 (MQ=255)                                                                                                                                                           gCC‑GGAACATATCGAAATGATGGCTCATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatccggcctacaatgggt                                                                                        >  2:592262‑M1/1‑113 (MQ=255)                                                                                                                                                                            aggatggcTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGCGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTGCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttggcgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgctcaatgcctgatgcgatggtttacgcatcttatcaggcttacaatggggaccggatcggtaggcc                                                                        <  1:228947‑M1/249‑155 (MQ=255)                                                                                                                                                                              gagggctaATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccg                                                                       <  2:502714‑M1/247‑156 (MQ=255)                                                                                                                                                                                    tAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTAGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatccggcctacaaccggtgccgcctccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataagg                                                                >  1:398168‑M1/1‑89 (MQ=255)                                                                                                                                                                                     aaTCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataag                                                                 <  2:166265‑M1/249‑162 (MQ=255)                                                                                                                                                                                      agctgggtgccgcggaTCGGTAAAGCGGTGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtaggcgtaggcattagacatttgtgcctgatgcgacgcttgtggcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctggtgcgatgcttgacgcatcttatcaggcctacaatgggtaccggagcggtaggccggataaggcg                                                              <  2:14671‑M1/244‑165 (MQ=255)                                                                                                                                                                                              gCCGCGTATCGCTAAAGCTGATCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgcc                                                      <  1:561831‑M1/251‑173 (MQ=255)                                                                                                                                                                                              gCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgcc                                                      <  2:147545‑M1/251‑173 (MQ=255)                                                                                                                                                                                               ccGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcgggaggccggataaggcgtttacgccg                                                     <  2:270062‑M1/251‑174 (MQ=255)                                                                                                                                                                                                   gTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggcccgataaggcgtttacgccgcatc                                                 >  1:299992‑M1/1‑74 (MQ=255)                                                                                                                                                                                                       cGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCAGCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccgg                                              <  1:235105‑M1/250‑181 (MQ=255)                                                                                                                                                                                                            aaGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcctccggcaaga                                         >  2:467956‑M1/1‑65 (MQ=255)                                                                                                                                                                                                             aGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGACTAGGCGGCGGTTTTAATCATCCGTGtagtcgtaggcattagacatttgggcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcagacctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaat                                       >  1:192125‑M1/1‑64 (MQ=255)                                                                                                                                                                                                                   gCGGTAAAGAACTGAGC‑GCCGTCGCGGAGAGCCGCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagca                                 <  2:582916‑M1/251‑194 (MQ=255)                                                                                                                                                                                                                      gTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctaca                                                                                              >  2:243934‑M1/1‑55 (MQ=255)                                                                                                                                                                                                                      gTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctaca                                                                                              <  1:243934‑M1/187‑133 (MQ=255)                                                                                                                                                                                                                         aaGAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcatacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagtt                           >  1:78266‑M1/1‑52 (MQ=255)                                                                                                                                                                                                                           gAACTGAGC‑GCCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatacgacgcttgccgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcctccggcaagaatagagcaccagtta                          >  1:590815‑M1/1‑50 (MQ=255)                                                                                                                                                                                                                                gAGC‑GCCGTCGAGAAGAGCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtacgcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggactacaaacggcgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccgggtcggtaggccggataaggcgtttccgccccatccggcaagaatagagcaccagttaaccga                     >  2:46209‑M1/1‑45 (MQ=255)                                                                                                                                                                                                                                     gcCGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaactt                 >  1:66671‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                                       cGTCGCGAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcaaaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttact              >  2:588219‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                                             gAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgc        >  1:282318‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                             gAAGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgc        >  1:406849‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                              aaGAGCCTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaagagagcaccagttaaccgaacttactctgcgcc       <  2:59430‑M1/251‑220 (MQ=255)                                                                                                                                                                                                                                                  gCCTCGGCGGCGGGTTTCATCATCCGTGtagtcgtaggctttagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggagccacatccggcaattggtgcaccatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtagggcgggaaaggcgtttacgcagcatcccgcaagaatagagcaacagttaaccgaacttactattcgcccaac   >  2:517517‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                                   ccTCGGCGGCGGTTTTCATCATCCGTGtagtcgtaggcattagacatttgtgcctgatgcgacgcttgacgcgtcttatcaggcctacaaccggtgccgcatccggcaattggtgcacaatgcctgatgcgatgcttgacgcatcttatcaggcctacaatgggtaccggatcggtaggccggataaggcgtttacgccgcatccggcaagaatagagcaccagttaaccgaacttactctgcgcccaaat  >  1:144709‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCC‑GGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGC‑GCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAGGGCGGAACTTCGCGAGTGAGATCATTATTAAACAGACAGTAATACGGCGACTCGCGGTCCGCGTCGTTGCTTAAATATGCCTCTTCGTAC  >  NC_000913/498815‑499305

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