Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,116	T→C	Y594H (TAT→CAT)	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,116	0	T	C	100.0%	174.5 / NA	55	Y594H (TAT→CAT)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/31);  total (24/31)

CAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTT  >  NC_000913/1869‑2359                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     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)                                                                                                                                                                                                                                                  ggtGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCttt  <  2:157664/249‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     CAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAT‑GACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTT  >  NC_000913/1869‑2359

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: ‑