Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,211	A→T	A625A (GCA→GCT)	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,211	0	A	T	100.0%	172.0 / NA	56	A625A (GCA→GCT)	thrA	fused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (27/29);  total (27/29)

ATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGC  >  NC_000913/1967‑2456                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      atcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAACTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgat                                                                                                                                                                                                                                                 <  1:391179/251‑1 (MQ=255)  tcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCTGGCAGTGGCGGATCAATTTGCCGACTTCTTGCGC‑GAAGGGTTCCACGTTGTTACGCCGAACAAAAAGGCCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATTGCGGCGTTAATTATTCCATGCCACCAACGTTGGGGATTGCTTACCGGTTATTATGAACCTGCAAAAACTTCTCAATGCTGGtgt                                                                                                                                                                                                                                                  >  2:37647/1‑248 (MQ=255)                cGGTCACTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATAAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCATCGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGtgc                                                                                                                                                                                                                                  <  2:58181/251‑3 (MQ=255)                                cACTTCCAGCCAGGCAGTGGCTGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCCCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                                                                                                                                                                                                                   <  1:540709/250‑2 (MQ=255)                                 aCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGAGGCGTTAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTTAGGACCTGCAAAATCTGCCCAATGCTGGTGATGAATTGATGATGTTCTCCCGCATTCTTTTTgg                                                                                                                                                                                                                 >  2:264255/1‑251 (MQ=255)                                    tCCAGCCAGTCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCTAGGTGGTCACGCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAGCCCGGCGTAAATTCCTCTATGACTCCAACGTTGGGGCTGGATTTCCTGTTATTGAGACCCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATTAATTTCGCCGGcgtggtttcaggtt                                                                                                                                                                                                               <  2:112668/250‑14 (MQ=255)                                         ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                                                                                                                                                                         <  1:673986/251‑1 (MQ=255)                                         ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                                                                                                                                                                         <  1:563405/251‑1 (MQ=255)                                             gtaGTGGCGGATCAATCTGCCGTCCTCCTGCGC‑GAATGTTTCCATGTTGTCTCGTCGAACTAAAAGGCCAACACCTCCTCGATGGATTACTACCAACTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATCCCTCTGTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTta                                                                                                                                                                                                     <  2:240711/249‑1 (MQ=255)                                                   gCGGATCAATAAGCCGACTTCTTG‑GCAGAAGGTTTCCACGTTGTCACTCCGAACAAAATGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCtt                                                                                                                                                                                               <  2:344015/251‑1 (MQ=255)                                                     ggATCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCg                                                                                                                                                                                             <  1:163745/251‑1 (MQ=255)                                                       aTCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc                                                                                                                                                                                           >  1:459564/1‑251 (MQ=255)                                                       aTCAATATGCCGACTTCTTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGAGGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGc                                                                                                                                                                                           <  1:291906/251‑1 (MQ=255)                                                                      tcttGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCACTGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGACGGc                                                                                                                                                                            >  2:237226/4‑251 (MQ=255)                                                                            gc‑gAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATAGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGt                                                                                                                                                                      >  2:612517/1‑251 (MQ=255)                                                                               gAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGttt                                                                                                                                                                    >  1:313831/1‑251 (MQ=255)                                                                                     ttCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATtgatga                                                                                                                                                                                                                                       <  1:10103/178‑1 (MQ=255)                                                                                        cACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                                                                                                                                      >  1:182474/1‑192 (MQ=255)                                                                                        cACGGTGTCACGCCGAACAAAAAGGCCAACACCTCGGCGATGGATTACTACCATCTGGTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTGATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                                                                                                                                      <  2:182474/192‑1 (MQ=255)                                                                                                  cGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGActactctg                                                                                                                                                 <  1:278493/251‑8 (MQ=255)                                                                                                        aCAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATAAGCTTCTCCGAGGCTACTACTCTGGCGCgg                                                                                                                                           >  2:460618/1‑251 (MQ=255)                                                                                                             aaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAAt                                                                                                                                      >  1:549721/1‑251 (MQ=255)                                                                                                                 ccAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGAGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGCTGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGt                                                                                                                                  >  1:614326/1‑251 (MQ=255)                                                                                                                        cTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCg                                                                                                                           >  1:327009/1‑251 (MQ=255)                                                                                                                        cTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCg                                                                                                                           >  1:373143/1‑251 (MQ=255)                                                                                                                             cgATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGGTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGGTATACCGAACCg                                                                                                                      >  2:147732/1‑251 (MQ=255)                                                                                                                                   aTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCcg                                                                                                                >  1:573551/1‑251 (MQ=255)                                                                                                                                                ggtGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCttt                                                                                                   <  2:157664/249‑1 (MQ=255)                                                                                                                                                               ggAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTGGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGAtgt                                                                                      >  1:448640/1‑249 (MQ=255)                                                                                                                                                               ggAAAAAACGCGGCGTAAATTCCTCTATGACTCCATCGTTGGGGCTGGATTCCCGGGTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTGCGGCAAGTTAGACGAAGGCATGAGGTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGAtgtag                                                                                    <  2:439176/251‑3 (MQ=255)                                                                                                                                                                               aaaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAACTCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATAAGGTTATCCGAGACGGCTACTCTAGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGGAAGCTATTGAt                                                                    >  1:53080/1‑251 (MQ=255)                                                                                                                                                                               aaaTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGAt                                                                    <  1:444113/251‑1 (MQ=255)                                                                                                                                                                                 aTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATtc                                                                  <  1:460618/251‑1 (MQ=255)                                                                                                                                                                                 aTTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATtc                                                                  <  1:440010/251‑1 (MQ=255)                                                                                                                                                                                    ccgctATGACACCATCGTTGTGGCTGGATTAACGGTTATTGAGAACCTGCAAAATCTGCTCAATGTTGGGGATGAATGGATGAAGTTCTCCGGCATTCTTTCATGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGGTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATtctc                                                                <  2:222564/247‑1 (MQ=255)                                                                                                                                                                                       ctATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGCGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTctcgctc                                                            >  1:230280/1‑251 (MQ=255)                                                                                                                                                                                         aggacacgtaCGTTGGTGCTGGTTTACAGGTTATTGATAACCTTCAAAATCTGCACAATGCTGGTGATGAAATGATGAAGTTCTCCGGCATGCTTTCATGTTCGCTTGCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAGCCGGATTCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGt                                                          <  2:363332/242‑1 (MQ=255)                                                                                                                                                                                         aTGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTctcgct                                                             <  1:261741/248‑1 (MQ=255)                                                                                                                                                                                          tGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAGCCTGCAAAATCTGCTCAATGCTGGTTATGAATTGATGAATTTCTCCGGCATTCTTTCAGGTGCGCTTTCTTATCTCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGAGCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGt                                                          <  2:363351/250‑1 (MQ=255)                                                                                                                                                                                               ccAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGAATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACg                                                    >  1:107464/1‑251 (MQ=255)                                                                                                                                                                                                    gTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGCTGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGCAACGGGACg                                               >  1:33220/1‑251 (MQ=255)                                                                                                                                                                                                     ttGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGt                                              >  1:291882/1‑251 (MQ=255)                                                                                                                                                                                                              gATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTgg                                       >  1:286622/1‑249 (MQ=255)                                                                                                                                                                                                              gATTACCGGTTATTGAGAACCTGCAAAAGCTGCTCAATGCTGGTGATGAATTGATGAAGTTCGCCGGCATTCTTTCAGGTTCGCCTTCTTATATCTTCGGCAAGGTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAg                                     <  2:566851/251‑1 (MQ=255)                                                                                                                                                                                                                atACCGGGTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATGCTTGCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTGTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGAACTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCt                                   <  2:113085/250‑1 (MQ=255)                                                                                                                                                                                                                        tatggagaACCTGCCAAATCTGCTCCAGGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGATTCGCTTTCTTATATCTTCGGCACGTTAGACGATGGCATGCGATTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTAGACCGAACCGGCGCCGCGCGCTGCTCTTTCTGGTATGGATGAAGAGCGTAAGCTATTGTTGCTCGATCGTGAAACGGGACGTGAACTGGAGCTGACGGat                            >  2:151121/5‑250 (MQ=255)                                                                                                                                                                                                                          ttGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATatt                         >  1:224017/1‑251 (MQ=255)                                                                                                                                                                                                                                   tGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGCTattgaaattga                 >  1:247719/1‑250 (MQ=255)                                                                                                                                                                                                                                    gCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATattgaaattgaa                <  1:420680/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATattgaaattga                 <  1:147732/249‑1 (MQ=255)                                                                                                                                                                                                                                    gCAAAATCTGCGCAATGCTGGTGATGTATTGATTAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGGTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGAGCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCGGGCGGATatttaaattgaa                <  2:459564/250‑1 (MQ=255)                                                                                                                                                                                                                                     cAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAAcc              <  1:196533/251‑1 (MQ=255)                                                                                                                                                                                                                                          tCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAAc                                                                                                                        >  2:362454/1‑140 (MQ=255)                                                                                                                                                                                                                                          tCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAAc                                                                                                                        <  1:362454/140‑1 (MQ=255)                                                                                                                                                                                                                                              cTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGcc     >  1:177021/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaTGCTGGGGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCTGGTATGGATGTAGCGCGTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGATGGCGGATATTGAAATTGAACCTGTGCTGCCCGc  >  2:479909/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGC‑GAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGC  >  NC_000913/1967‑2456

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