Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	909,056	T→C	Q80R (CAG→CGG)	ltaE ←	L‑allo‑threonine aldolase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	909,056	0	T	C	100.0%	193.0 / NA	54	Q80R (CAG→CGG)	ltaE	L‑allo‑threonine aldolase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/30);  total (24/30)

CGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGT  >  NC_000913/908909‑909200                                                                                                                                                    |                                                                                                                                                 cGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGAc                                                                                                                                               <  2:379896/151‑1 (MQ=255)       aaaTGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACaa                                                                                                                                         <  2:47177/151‑1 (MQ=255)       aaaTGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACa                                                                                                                                          >  2:44899/1‑150 (MQ=255)       aaaTGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACa                                                                                                                                          <  1:285419/150‑1 (MQ=255)             atatCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATAc                                                                                                                                    <  1:519003/150‑1 (MQ=255)                   tcgGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTcgc                                                                                                                             >  2:153697/1‑151 (MQ=255)                     gggTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTcgccg                                                                                                                           >  2:204121/1‑151 (MQ=255)                     gggTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTcgcc                                                                                                                            >  1:312947/1‑150 (MQ=255)                       gTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCcgcg                                                                                                                         <  2:91011/151‑1 (MQ=255)                       gTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCcgcg                                                                                                                         <  2:271661/151‑1 (MQ=255)                          tGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTc                                                                                                                      >  1:317082/1‑151 (MQ=255)                            aTTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGc                                                                                                                    <  2:617532/151‑1 (MQ=255)                               ttCATCGCCCCTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAg                                                                                                                  <  1:166531/150‑1 (MQ=255)                                tCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTg                                                                                                                >  1:602689/1‑151 (MQ=255)                                    cGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTg                                                                                                            >  2:159710/1‑151 (MQ=255)                                    cGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTg                                                                                                            <  2:395052/151‑1 (MQ=255)                                     gCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTg                                                                                                            >  1:374378/1‑150 (MQ=255)                                          tttATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCAGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGc                                                                                                         >  2:15964/1‑148 (MQ=255)                                              tCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCg                                                                                                   >  2:327263/1‑150 (MQ=255)                                                   cGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAgg                                                                                             <  1:631183/151‑1 (MQ=255)                                                         cGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGcc                                                                                       >  1:422468/1‑151 (MQ=255)                                                               gTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGt                                                                                  <  1:562121/150‑1 (MQ=255)                                                                  ggCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCg                                                                              <  2:679142/151‑1 (MQ=255)                                                                  ggCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCg                                                                              <  1:94245/151‑1 (MQ=255)                                                                        cgcgTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGc                                                                        >  2:275408/1‑151 (MQ=255)                                                                         gcgTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCa                                                                       <  2:668141/151‑1 (MQ=255)                                                                               aTGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGaaaa                                                                  >  2:673750/1‑150 (MQ=255)                                                                                tGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGaa                                                                    <  1:102107/147‑1 (MQ=255)                                                                                 gggTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATg                                                               <  1:681608/151‑1 (MQ=255)                                                                                 gggTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATg                                                               <  1:384949/151‑1 (MQ=255)                                                                                   gTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATggc                                                             <  1:210054/151‑1 (MQ=255)                                                                                     tGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATggcg                                                            <  2:192259/150‑1 (MQ=255)                                                                                       cGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAGAATGGCGGct                                                         >  2:605051/1‑151 (MQ=255)                                                                                           tGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGcttctt                                                     >  2:442034/1‑151 (MQ=255)                                                                                            gAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCttt                                                    >  1:458620/1‑151 (MQ=255)                                                                                                aCTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCg                                                <  2:370950/151‑1 (MQ=255)                                                                                                  tGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGGCTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCgg                                               >  2:682906/1‑150 (MQ=255)                                                                                                   gCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGaa                                             <  2:88315/151‑1 (MQ=255)                                                                                                    cccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGaa                                             <  1:634329/150‑1 (MQ=255)                                                                                                     ccAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGaa                                             >  1:667487/1‑149 (MQ=255)                                                                                                      cAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGc                                          >  1:492557/1‑151 (MQ=255)                                                                                                         cACCGCCGCGCCGCCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGctct                                       >  1:347617/1‑151 (MQ=255)                                                                                                                 cgcCACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGt                                >  1:538252/1‑150 (MQ=255)                                                                                                                   ccACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGt                             <  2:267532/151‑1 (MQ=255)                                                                                                                   ccACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGt                             <  2:585553/151‑1 (MQ=255)                                                                                                                    cACCGGCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTc                            >  2:411732/1‑151 (MQ=255)                                                                                                                         gCTTCAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCa                       <  1:480703/151‑1 (MQ=255)                                                                                                                             cAAACAGATAGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCagag                    <  1:411719/150‑1 (MQ=255)                                                                                                                                      aGTTATGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACgg           <  1:649630/150‑1 (MQ=255)                                                                                                                                           tGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAggg      >  2:600772/1‑150 (MQ=255)                                                                                                                                           tGCGCGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAgg       <  1:422214/149‑1 (MQ=255)                                                                                                                                            gcgcGGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGgt     <  1:175715/150‑1 (MQ=255)                                                                                                                                              gcgGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGgtcgt  <  2:502235/151‑1 (MQ=255)                                                                                                                                              gcgGCCCGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGgtcg   >  1:397341/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CGGGCGAAATGGATATCGTCGGGTTTGATTTTCATCGCCACTTTATCCAGCGGTAGCGTGCCGTCGGCAGCCGCGTCTATGGGTTGCGGTTGAATACTGCCCAGCACCGCCGCGCCACCGGCTTCAAACAGATAGTTATGCGCGGCCTGACCGACAATATACTCTTCGCCGCGTTCGCAGTGACTGAGCAGAGCGACCAGGTTGGCCTGAGTGCCGGTAGGCAGAAAAATGGCGGCTTCTTTACCGGAAAGCTCTGCTGCGTAGTCCTGCAGAGCATTAACGGTAGGGTCGT  >  NC_000913/908909‑909200

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