Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,108	A→T	Q591L (CAG→CTG)	thrA →	fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,108	0	A	T	100.0%	183.6 / NA	58	Q591L (CAG→CTG)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (26/32);  total (26/32)

GCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCG  >  NC_000913/1863‑2346                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                                                                                             <  1:217660/250‑1 (MQ=255) gCTCAACAATGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                                                                                             <  1:17525/250‑1 (MQ=255)     aCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGtc                                                                                                                                                                                                                                         >  1:274721/1‑249 (MQ=255)     aCCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTAGTCGATGGATTACTACCATCTGTTGCGtc                                                                                                                                                                                                                                         >  1:125735/1‑249 (MQ=255)       cAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATg                                                                                                                                                                                                                                      <  1:599376/251‑1 (MQ=255)                              aaaaCTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaaa                                                                                                                                                                                                               <  1:497827/251‑1 (MQ=255)                              aaaaCTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGATGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaa                                                                                                                                                                                                                >  1:157664/1‑250 (MQ=255)                                    ggCAGTAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCct                                                                                                                                                                                                          <  1:489880/250‑1 (MQ=255)                                     gCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTCACATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCctct                                                                                                                                                                                                        >  1:222564/1‑251 (MQ=255)                                       aGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAt                                                                                                                                                                                                      <  1:477947/251‑1 (MQ=255)                                              aCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGACAGCATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGAGGCTCAATATGCCGGCTTCTTGCGC‑GAAGGCTTCCACGTTGTAACGCCGACCAAAAAGGCCCACACCTCGTCGCTGGATTACTACCATATGTTGCGTCTTGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACAccc                                                                                                                                                                                              >  1:616665/1‑250 (MQ=255)                                              aCTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCa                                                                                                                                                                                              <  1:667552/251‑1 (MQ=255)                                               ctagcgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTCACGCCGAACCAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGGAAAATCGCGGCGTAAATTTCTCTAT‑GACACCa                                                                                                                                                                                              >  1:363332/4‑250 (MQ=255)                                               cTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTCAATTCCTCTAT‑GACACCa                                                                                                                                                                                              >  1:363351/1‑250 (MQ=255)                                                   cgcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGtt                                                                                                                                                                                         <  1:386768/251‑1 (MQ=255)                                                     cAAGCCAAAGAGCCGTTTAAGCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTgg                                                                                                                                                                                       <  2:274721/251‑1 (MQ=255)                                                       aGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCTGTCAATGTTTACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTgggg                                                                                                                                                                                     <  2:329781/251‑1 (MQ=255)                                                            aatagctgTTTAATCTCGGGTGCTTAATTCACCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAGTGTTGACTGCACTTCCAGCCAGGCAGTTGCGGATCAATATGCCGACTTCTTGCGC‑GAACGTTTCCACGTTGTCACGCCGAACAAAACGGCCAACACCTCGTCGATGTATTACTACCATCTGTTGCGTCATGCGGCGGATAAATCGCGGCGTAATTTCCTCTAT‑GACTACAACGTTGGGGCTGGa                                                                                                                                                                                <  2:634879/244‑1 (MQ=255)                                                             agtgCCGTTTAATCGCGGGCGCTTAACTCGCCTCGTGCAAGAATATCATCTGCTGAACCCTGTCTTTGTTGACTGCACGTCCAGCCAGGCAGTGGTGGATCAATATGCCGACTTCTTGAGC‑GAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACTTCGGCGATCGATGACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACAAACGTTGGGGCTGGAt                                                                                                                                                                               <  2:378867/248‑1 (MQ=255)                                                              gagCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCCACGTTGGGGCTGGAtt                                                                                                                                                                              >  1:566851/1‑251 (MQ=255)                                                                 ccGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTAcc                                                                                                                                                                           <  1:384125/251‑1 (MQ=255)                                                                   gTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCgg                                                                                                                                                                         <  1:280962/251‑1 (MQ=255)                                                                                        tCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAAt                                                                                                                                                     <  1:576554/250‑1 (MQ=255)                                                                                           ccTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGc                                                                                                                                                 <  1:378470/251‑1 (MQ=255)                                                                                                 gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTaaa                                                                                                                                                                                                               <  1:460972/184‑1 (MQ=255)                                                                                                 gAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAAAATGCCGACTTCTTGCGC‑GAAGGTTTCCCCGTTGTCCCGCCGAAAAAAAAGGCCAACACCTCGTCGCTGGATTACTACCAACTGGTGCGTCATGCTGCGGAAAAATCGCGGCGTaaa                                                                                                                                                                                                               >  2:460972/1‑184 (MQ=255)                                                                                                    aGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGctg                                                                                                                                        <  1:17787/251‑3 (MQ=255)                                                                                                         atcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAACTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgat                                                                                                                                   <  1:391179/251‑1 (MQ=255)                                                                                                          tcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCTGGCAGTGGCGGATCAATTTGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTTACGCCGAACAAAAAGGCCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATTGCGGCGTTAATTATTCCAT‑GCCACCAACGTTGGGGATTGCTTACCGGTTATTATGAACCTGCAAAAACTTCTCAATGCTGGtgt                                                                                                                                    >  2:37647/1‑248 (MQ=255)                                                                                                                gctgAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTCACGCCGAACAAACAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATAGCGGCGAGAATTCCTCTATTCCCACCA‑CGTTGGGGCTGGATTACCGGTTATTTAGAACCTGCAAAATCTGCGCaaagatgggggagacttgg                                                                                                                            >  2:196533/1‑234 (MQ=255)                                                                                                                        cGGTCACTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATAAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCATCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGtgc                                                                                                                    <  2:58181/251‑3 (MQ=255)                                                                                                                                        cACTTCCAGCCAGGCAGTGGCTGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCCCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                                                                                                     <  1:540709/250‑2 (MQ=255)                                                                                                                                         aCTTCCAGCCAGGCAGTTGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGGCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTGGCGTCATGCGGCGTAAAAATCGCGGCGTAAATTCCTCTAT‑GTCACCAACGTTGGGGCTGGAGTACCGGTTATTGAGaat                                                                                                                                                              <  2:628942/192‑2 (MQ=255)                                                                                                                                         aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGCACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTCT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAAc                                                                                                                                                              >  1:628942/1‑192 (MQ=255)                                                                                                                                         aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGAGGCGTTAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTTAGGACCTGCAAAATCTGCCCAATGCTGGTGATGAATTGATGATGTTCTCCCGCATTCTTTTTgg                                                                                                   >  2:264255/1‑251 (MQ=255)                                                                                                                                            tCCAGCCAGTCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCTAGGTGGTCACGCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCCCGGCGTAAATTCCTCTAT‑GACTCCAACGTTGGGGCTGGATTTCCTGTTATTGAGACCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATTAATTTCGCCGGcgtggtttcaggtt                                                                                                 <  2:112668/250‑14 (MQ=255)                                                                                                                                                 ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                                                           <  1:673986/251‑1 (MQ=255)                                                                                                                                                 ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                                                           <  1:563405/251‑1 (MQ=255)                                                                                                                                                     gtaGTGGCGGATCAATCTGCCGTCCTCCTGCGC‑GAATGTTTCCATGTTGTCTCGTCGAACTAAAAGGCCAACACCTCCTCGATGGATTACTACCAACTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATCCCTCTGT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTta                                                                                       <  2:240711/249‑1 (MQ=255)                                                                                                                                                           gCGGATCAATAAGCCGACTTCTTG‑GCAGAAGGTTTCCACGTTGTCACTCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCtt                                                                                 <  2:344015/251‑1 (MQ=255)                                                                                                                                                             ggATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCg                                                                               <  1:163745/251‑1 (MQ=255)                                                                                                                                                               aTCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc                                                                             >  1:459564/1‑251 (MQ=255)                                                                                                                                                               aTCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGAGGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc                                                                             <  1:291906/251‑1 (MQ=255)                                                                                                                                                                 cAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATACGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAAt                                                                                                                                            >  1:96955/1‑186 (MQ=255)                                                                                                                                                                              tcttGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCACTGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGACGGc                                                              >  2:237226/4‑251 (MQ=255)                                                                                                                                                                                    gc‑gAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATAGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGt                                                        >  2:612517/1‑251 (MQ=255)                                                                                                                                                                                       gAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt                                                      >  1:313831/1‑251 (MQ=255)                                                                                                                                                                                             ttCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATtgatga                                                                                                                         <  1:10103/178‑1 (MQ=255)                                                                                                                                                                                                cACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                        >  1:182474/1‑192 (MQ=255)                                                                                                                                                                                                cACGGTGTCACGCCGAACAAAAAGGCCAACACCTCGGCGATGGATTACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTGATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                        <  2:182474/192‑1 (MQ=255)                                                                                                                                                                                                          cGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGActactctg                                   <  1:278493/251‑8 (MQ=255)                                                                                                                                                                                                                aCAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATAAGCTTCTCCGAGGCTACTACTCTGGCGCgg                             >  2:460618/1‑251 (MQ=255)                                                                                                                                                                                                                     aaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAAt                        >  1:549721/1‑251 (MQ=255)                                                                                                                                                                                                                         ccAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGAGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGCTGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGt                    >  1:614326/1‑251 (MQ=255)                                                                                                                                                                                                                                cTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCg             >  1:327009/1‑251 (MQ=255)                                                                                                                                                                                                                                cTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCg             >  1:373143/1‑251 (MQ=255)                                                                                                                                                                                                                                     cgATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGGTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGGTATACCGAACCg        >  2:147732/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCcg  >  1:573551/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCG  >  NC_000913/1863‑2346

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