Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,599,441	T→C	D73G (GAC→GGC)	livJ ←	branched chain amino acid/phenylalanine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,599,441	0	T	C	100.0%	79.1 / NA	24	D73G (GAC→GGC)	livJ	branched chain amino acid/phenylalanine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/6);  total (18/6)

TCGTCTTCGTAGATGTCAGACGCAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGGTCACCGTACTGCGCAACCGGACCGGACATTGCGC  >  NC_000913/3599315‑3599564                                                                                                                                |                                                                                                                            tcgtcttcgtAGATGTCAGACGCAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGGCACCGCAATCGCATGTTTCGGGTCAGAGGCGTCGCCATAtt                                                                                                                       >  2:143195/1‑133 (MQ=255) tcgtcttcgtAGATGTCAGACGCAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                                                       >  1:296750/1‑133 (MQ=255) tcgtcttcgtAGATGTCAGACGCAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                                                       >  1:475141/1‑133 (MQ=255)   gtcttcgtAGATGTCAGCCGCAGGCTGCGTTGATGAAGAACAGAGGTGCCCAATCACATATTTAAGGCCGTCGTTAA‑CGCCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtttt                                                                                                                     <  2:434295/133‑1 (MQ=255)            gATGTCAGACGCAGGCTGCTTTGATAAAGAACAGAGGCGACCAATCACATATTTAACGCCGTCTTTAA‑CGATTTTGTTCCCCCCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGc                                                                                                            <  2:474435/133‑1 (MQ=255)                  aGACGCAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                                                      >  2:29734/1‑133 (MQ=255)                       cAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGcc                                                                                                 >  1:18641/1‑133 (MQ=255)                          gCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTGTACGATTTTCAGTATGTTGCCttt                                                                                              >  2:142149/1‑133 (MQ=255)                                     gaagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                                                   <  2:475141/133‑1 (MQ=255)                                     gaagaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCtt                                                                                   <  1:29734/133‑1 (MQ=255)                                                                 tAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                                                       >  1:380194/1‑133 (MQ=255)                                                                      ccGTCGTTAACCGTTTTTTTCCCCCCCCGAACCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgcg                                                   <  2:5021/133‑1 (MQ=255)                                                                      ccGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgcgc                                                  >  2:353069/1‑133 (MQ=255)                                                                        gtcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCg                                                >  1:5021/1‑133 (MQ=255)                                                                         tcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAGCCGCCTGCTCTGCGCCgg                                               >  1:249853/1‑133 (MQ=255)                                                                         tcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCgg                                               >  1:256776/1‑133 (MQ=255)                                                                         tcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCAGCAACCGCCTGCTTTGCGCCgg                                               >  1:211204/1‑133 (MQ=255)                                                                         tcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTCATACACGCAACCACCTGAGCTGCGccag                                               >  2:245413/1‑131 (MQ=255)                                                                         tcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGGTGATATCCGCAACCTCCTGCTCTGCGCCgg                                               >  1:236031/1‑133 (MQ=255)                                                                         tcgtTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGCTGATATCAGCAACCTCCTGCTTTGCGCCgg                                               >  1:184014/1‑133 (MQ=255)                                                                                 cGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTc                                        >  1:466553/1‑133 (MQ=255)                                                                                                      ccGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGGTCACCGTACTGCGCAac                   <  2:26148/133‑1 (MQ=255)                                                                                                          cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGGTCACCGTACTGCGCAaccgga               >  1:386569/1‑133 (MQ=255)                                                                                                                       caGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGGTCACCGTACTGCGCAACCGGACCGGACATTgcgc  >  2:67594/1‑133 (MQ=255)                                                                                                                                |                                                                                                                            TCGTCTTCGTAGATGTCAGACGCAGGCTGCGTTGATGAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAA‑CGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGGTCACCGTACTGCGCAACCGGACCGGACATTGCGC  >  NC_000913/3599315‑3599564

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