Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,130,798	T→C	F322L (TTT→CTT)	metL →	fused aspartate kinase/homoserine dehydrogenase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,130,798	0	T	C	100.0%	156.2 / NA	48	F322L (TTT→CTT)	metL	fused aspartate kinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/23);  total (25/23)

CGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACA  >  NC_000913/4130671‑4130927                                                                                                                                |                                                                                                                                  caaaaTCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCg                                                                                                                       <  1:81923/138‑1 (MQ=255) cGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCg                                                                                                                       <  1:392408/140‑1 (MQ=255)  aaaaTCGCCCTGCAACTGCGCTGTAGCTACCCGCCGGATAAAGGTCCCACGCGCTTTGACCGCGTGCTGCCCCCCGGTATTGGTGCGCTTTTTGTCCCCAGCCACGATGATGTCTGTTTGATTGAGCTTAAGGTGCCCGc                                                                                                                      <  2:966258/139‑1 (MQ=255)      ccgcccTGCAACTGCGCTGTAGCTCCACGCCGAATAAGGGTCCCCCGCGCATTGACCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCAAAATGATGTCTTTTTGATTGAGCTTCAGGTGCCCGCCAGt                                                                                                                  <  2:399524/136‑1 (MQ=255)       cGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTG‑CCGCCAGTCa                                                                                                                <  2:422139/140‑1 (MQ=255)                    cgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACt                                                                                                    <  1:429226/140‑1 (MQ=255)                       tGTAGCTACACGCGGGATCAAGGTTCCACGCGCTTTGAACGGGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCCCCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                                                                                 <  2:693848/140‑1 (MQ=255)                        gTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCCCCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgcg                                                                                                <  2:98290/140‑1 (MQ=255)                        gTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgcg                                                                                                <  1:204249/140‑1 (MQ=255)                         tAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgcgc                                                                                               >  2:935408/1‑140 (MQ=255)                                 cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAagag                                                                                       >  1:487744/1‑140 (MQ=255)                                   ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                                                                                     >  1:549483/1‑140 (MQ=255)                                   ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                                                                                     >  2:342386/1‑140 (MQ=255)                                   ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                                                                                     >  1:803458/1‑140 (MQ=255)                                   ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                                                                                     >  1:803317/1‑140 (MQ=255)                                               tCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCt                                                                         <  2:549483/140‑1 (MQ=255)                                                ccACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTg                                                                        >  2:375835/1‑140 (MQ=255)                                                    gcgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc                                                                    <  2:803458/140‑1 (MQ=255)                                                    gcgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc                                                                    >  1:390086/1‑140 (MQ=255)                                                    gcgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc                                                                    <  2:616136/140‑1 (MQ=255)                                                      gcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc                                                                  >  1:78577/1‑140 (MQ=255)                                                       cATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg                                                                 >  2:623725/1‑140 (MQ=255)                                                                cgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTAc                                                        >  1:550214/1‑140 (MQ=255)                                                                  tGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGc                                                      >  2:345746/1‑140 (MQ=255)                                                                  tGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGc                                                      >  2:334996/1‑140 (MQ=255)                                                                           ccGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcgcg                                                               <  2:797169/122‑1 (MQ=255)                                                                           ccGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcgcg                                                               >  1:797169/1‑122 (MQ=255)                                                                            cGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            >  1:894656/1‑140 (MQ=255)                                                                            cGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            <  2:487744/140‑1 (MQ=255)                                                                            cGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            <  2:894656/140‑1 (MQ=255)                                                                            cGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            <  2:78577/140‑1 (MQ=255)                                                                            cGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            <  2:803317/140‑1 (MQ=255)                                                                                    gTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            >  1:851783/1‑132 (MQ=255)                                                                                    gTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCgg                                            <  2:851783/132‑1 (MQ=255)                                                                                      gcgcGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACa                                  >  1:258076/1‑140 (MQ=255)                                                                                       cgcgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGGGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGt                                           <  2:710051/130‑1 (MQ=255)                                                                                       cgcgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGGGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGt                                           >  1:710051/1‑130 (MQ=255)                                                                                       cgcgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACAt                                 >  1:950367/1‑140 (MQ=255)                                                                                         cgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATaa                               <  1:375835/140‑1 (MQ=255)                                                                                           tATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTg                                         <  2:614143/128‑1 (MQ=255)                                                                                           tATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTg                                         >  1:614143/1‑128 (MQ=255)                                                                                                 cACCAGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGcc                       <  1:935408/140‑1 (MQ=255)                                                                                                     aGCCACGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGtt                   >  1:307750/1‑140 (MQ=255)                                                                                                          cGATGATGTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTtgctgc              <  1:345746/140‑1 (MQ=255)                                                                                                                tgTCTGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAAtttt        >  2:29566/1‑140 (MQ=255)                                                                                                                    tGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTa    <  2:307750/140‑1 (MQ=255)                                                                                                                    tGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTa    <  2:950367/140‑1 (MQ=255)                                                                                                                    tGTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTa    <  2:258076/140‑1 (MQ=255)                                                                                                                     gTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTtgctgc              <  1:319543/129‑1 (MQ=255)                                                                                                                     gTTTGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTtgctgc              >  2:319543/1‑129 (MQ=255)                                                                                                                      tttGATTGAGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTaca  >  1:776961/1‑140 (MQ=255)                                                                                                                       ttgattgaGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTa    >  1:584031/1‑137 (MQ=255)                                                                                                                       ttgattgaGCTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTa    <  2:584031/137‑1 (MQ=255)                                                                                                                                |                                                                                                                                  CGAAATCGACCTGCAACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACA  >  NC_000913/4130671‑4130927

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