Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,374,924	C→A	R245L (CGG→CTG)	epmB ←	lysine 2,3‑aminomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,374,924	0	C	A	100.0%	85.2 / NA	25	R245L (CGG→CTG)	epmB	lysine 2,3‑aminomutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (14/11);  total (14/11)

TATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGCCTCGGTGATACGTGCCGGGA  >  NC_000913/4374793‑4375056                                                                                                                                    |                                                                                                                                     tATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAAc                                                                                                                              >  1:86788/1‑140 (MQ=255)           aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCg                                                                                                                    >  2:341098/1‑140 (MQ=255)                     tACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGaa                                                                                                          >  2:83684/1‑140 (MQ=255)                       cGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATg                                                                                                        <  1:269081/140‑1 (MQ=255)                             tACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCAt                                                                                                  <  1:306176/140‑1 (MQ=255)                                 ccGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAc                                                                                              <  2:185966/140‑1 (MQ=255)                                     cATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCa                                                                                          <  1:239322/140‑1 (MQ=255)                                      aTCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCAt                                                                                         >  1:228404/1‑140 (MQ=255)                                           acaacGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCa                                                                                    <  1:58285/140‑1 (MQ=255)                                                         aGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGAtt                                                                      >  2:77606/1‑140 (MQ=255)                                                              tGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCa                                                                 >  1:289896/1‑140 (MQ=255)                                                                                aTCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAg                                               <  2:357134/140‑1 (MQ=255)                                                                                 tCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGa                                              >  2:259464/1‑140 (MQ=255)                                                                                              gTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGCAGAACGGGCAAAGCAt                                 >  1:95693/1‑140 (MQ=255)                                                                                                     aaCGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCa                          >  1:43371/1‑140 (MQ=255)                                                                                                      aCGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAg                         <  1:89688/140‑1 (MQ=255)                                                                                                      aCGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAg                         <  2:228404/140‑1 (MQ=255)                                                                                                      aCGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAg                         >  1:70675/1‑140 (MQ=255)                                                                                                      aCGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAg                         >  1:333183/1‑140 (MQ=255)                                                                                                       cGCTCTGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAgc                        >  2:295767/1‑140 (MQ=255)                                                                                                          tctGGTTCAGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGcc                     >  2:278374/1‑140 (MQ=255)                                                                                                                  aGCAAAGTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGCCTCGGTGAt             <  2:108828/140‑1 (MQ=255)                                                                                                                        gTAACGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGCCTCGGTGATACGTGc       <  2:96482/140‑1 (MQ=255)                                                                                                                            cGCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGCCTCGGTGATACGTGCCggg   >  1:121755/1‑140 (MQ=255)                                                                                                                             gCCTACCAGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGCCTCGGTGATACGTGCCGGGa  <  1:341098/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     TATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCATTGGCATGGTTGATGTGATTCACCAGCAAGATTTGCAGCGTAGAACGGGCAAAGCATTCAACCAGCGCCTCGGTGATACGTGCCGGGA  >  NC_000913/4374793‑4375056

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