Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,492	T→C	A307A (GCA→GCG)	yjeH ←	L‑methionine/branched chain amino acid exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,492	0	T	C	100.0%	168.5 / NA	51	A307A (GCA→GCG)	yjeH	L‑methionine/branched chain amino acid exporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/29);  total (22/29)

CCATTGGCATAAATAATAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGTGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATAACGCTC  >  NC_000913/4369367‑4369616                                                                                                                              |                                                                                                                             ccATTGGCATAAATAATAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGt                                                                                                                       >  2:258692/1‑133 (MQ=255)   aTTGGCATAAATAATAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAg                                                                                                                     >  1:275906/1‑133 (MQ=255)       gtataaataaTAAGAGCATCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGTATTATGCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGt                                                                                                                 <  2:354237/131‑1 (MQ=255)       gCATAAATAATAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGt                                                                                                                 <  2:448686/133‑1 (MQ=255)       gCATAAATAATAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGt                                                                                                                 <  2:337957/133‑1 (MQ=255)         ataaataaTATGAGCCTCCAGATTCATCGCTAAAGCATGAACCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCa                                                                                                               <  2:152844/133‑1 (MQ=255)         ataaataaTAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCa                                                                                                               <  1:463897/133‑1 (MQ=255)             ataataAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGAt                                                                                                           <  1:258692/133‑1 (MQ=255)                     agCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTgc                                                                                                   >  1:430989/1‑133 (MQ=255)                       cGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTgcgc                                                                                                 >  1:40788/1‑133 (MQ=255)                       cGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTgcgc                                                                                                 >  2:20950/1‑133 (MQ=255)                       cGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGCAGTAGTCAGGATTCTGTTgcgc                                                                                                 >  1:217170/1‑133 (MQ=255)                       cGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAAGTAGTGGTCAGGAATATGTTgcgc                                                                                                 >  2:156868/1‑133 (MQ=255)                        gTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGcc                                                                                                >  1:422215/1‑133 (MQ=255)                         tCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGCGAGGCGCGCCAGGTAGTGGCCAGGCTTATGTTgccgct                                                                                               >  2:353947/1‑129 (MQ=255)                         tCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCt                                                                                               >  1:354508/1‑133 (MQ=255)                         tCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCt                                                                                               >  2:506336/1‑133 (MQ=255)                         tCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCt                                                                                               >  1:82652/1‑133 (MQ=255)                         tCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCACCAGGTAGTGGTCAGGATTATGTTGCGCCt                                                                                               >  1:359118/1‑133 (MQ=255)                                  gtctcTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGACGCACGCCAGGTAGTGGTCAGGATTATGTAGCGCCTGCgaccaga                                                                                      >  2:232995/2‑133 (MQ=255)                                            cATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATGATGTTGCGCCTACGACCAGACCAggcgggc                                                                            >  1:202890/1‑133 (MQ=255)                                             aTGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAggcgggcg                                                                           <  2:505667/133‑1 (MQ=255)                                              tGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGa                                                                          >  1:368923/1‑133 (MQ=255)                                                aaTCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAg                                                                        <  2:39787/133‑1 (MQ=255)                                                aaTCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAg                                                                        <  2:40788/133‑1 (MQ=255)                                                 aTCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGc                                                                       <  2:433250/133‑1 (MQ=255)                                                  tCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGct                                                                      <  2:354508/133‑1 (MQ=255)                                                  tCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGct                                                                      <  2:350213/133‑1 (MQ=255)                                                  tCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGct                                                                      <  1:506336/133‑1 (MQ=255)                                                  tCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGct                                                                      <  2:40695/133‑1 (MQ=255)                                                   cACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGctc                                                                     <  2:275906/133‑1 (MQ=255)                                                    aCCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGctct                                                                    <  2:305618/133‑1 (MQ=255)                                                     ccAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTg                                                                   <  2:49122/133‑1 (MQ=255)                                                      cAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGGGAAACCCTGt                                                                  >  2:89737/1‑133 (MQ=255)                                                      cAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGt                                                                  <  2:368923/133‑1 (MQ=255)                                                              tCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGGAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGGACTGTATATAAAt                                                          >  2:281369/1‑133 (MQ=255)                                                                accacACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGTGCGGAGATCTGTaaataaatat                                                        >  1:169775/1‑131 (MQ=255)                                                                accacACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGta                                                                 <  2:21939/124‑1 (MQ=255)                                                                accacACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGta                                                                 >  1:21939/1‑124 (MQ=255)                                                                         cagcCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTgg                                               <  2:430989/133‑1 (MQ=255)                                                                              gAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAg                                          <  1:115111/133‑1 (MQ=255)                                                                               aGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAAAGTTGAGACTGGCAAAGc                                         >  1:369972/1‑133 (MQ=255)                                                                                   ccGCATTGAGGGCATAATTCGGGATATGGCGAAAAGAGAGGTGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAATCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGc                                     <  2:169775/133‑1 (MQ=255)                                                                                       aTTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGa                                 <  2:82652/133‑1 (MQ=255)                                                                                              gCATTATTCGGGATATGGCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCaa                          >  1:41029/1‑133 (MQ=255)                                                                                                   attCGGGAGATGGCGCGACGAGAGGCGCGCTAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACg                     <  1:156868/133‑1 (MQ=255)                                                                                                              ggCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGGAGGTCTGTATATAAATGTTGAGACTGGCAAAACAGGCCAAATAGCCAAACACGGCGGCAATccc          >  2:205232/1‑132 (MQ=255)                                                                                                              ggCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCCATCACGCAGGGAAACCa          >  2:66586/1‑133 (MQ=255)                                                                                                               gCGAGAAGAGAGGCGCGCCAGGTAGTGGTCAGGAATATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCAt         <  1:326815/133‑1 (MQ=255)                                                                                                                 gagaAGAGGGGCGCGGCAGGTATTGGTCAGAATTATGTTGCGCCTGCGACCAGACCAGGCGCGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATaa       <  1:232995/133‑1 (MQ=255)                                                                                                                 gagaAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATaa       <  2:481408/133‑1 (MQ=255)                                                                                                                 gagaAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATaa       <  1:89737/133‑1 (MQ=255)                                                                                                                 gagaAGAGAGGCGCGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATaa       <  1:353947/133‑1 (MQ=255)                                                                                                                 gagaAGAGAGGCGCACCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATaa       <  2:359118/133‑1 (MQ=255)                                                                                                                      gagagGCGCGCCAGGCAGTGGTCAGGATTATGCTGCGCCTGAGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATAACGCTc  <  2:173202/133‑1 (MQ=255)                                                                                                                              |                                                                                                                             CCATTGGCATAAATAATAAGAGCGTCCAGATTGATCTCTAAAGCATGAATCACCAAAGTGCTCACCACACAGCAGCCGAGCACCGCATTGAGGGCATTATTCGGGATATGGCGAGAAGAGAGGCGTGCCAGGTAGTGGTCAGGATTATGTTGCGCCTGCGACCAGACCAGGCGGGCGAAGCTCTGTATATAAATGTTGAGACTGGCAAAGCAGGCCAGATAGCCAATCACGCAGGCAATCCATAACGCTC  >  NC_000913/4369367‑4369616

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