Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,891,231	G→A	L332L (CTG→CTA)	tnaB →	tryptophan:H(+) symporter TnaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,891,231	0	G	A	100.0%	83.8 / NA	29	L332L (CTG→CTA)	tnaB	tryptophan:H(+) symporter TnaB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/21);  total (8/21)

AATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTGCCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTA  >  NC_000913/3890986‑3891472                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aaTATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTaccacc                                                                                                                                                                                                                                                 <  2:360649/251‑1 (MQ=255)                    gtAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATGATCGAGTTTTGCCTGCTGGTGTTCTCTAACGTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGa                                                                                                                                                                                                                                <  2:153035/248‑1 (MQ=255)                             ttATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGAACTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTCCTTTTAACCTTCCTACCACCCGCGTTGTTGTATCTGATCTTcccacac                                                                                                                                                                                                                     >  2:335000/1‑247 (MQ=255)                                               aCGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATc                                                                                                                                                                                                   <  2:258971/251‑1 (MQ=255)                                                             gTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGGCC‑gggg                                                                                                                                                                                    >  2:365903/1‑250 (MQ=255)                                                                    cgtgccTCGGCTCCAACCAGCACGTCATTATCTAGTTTTGCCTGCTGGTGTTCTCTATCTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGTGTTGTTCGATTATCTGGCGGACCTGTTTCAGTTTGATAACTCCCCCGGCGGTCATTTCAAAGCCTTGCTGTTGACCCTCCTACCTCCTGCGTTTTTGTATCTTATCTTCCCGAACGGCCTTATTTACGGGATCGGCGGTG‑CC‑GGGCTCTGCGcc                                                                                                                                                                            <  2:133712/247‑1 (MQ=255)                                                                     gTTCCTCGGCACCAAACAGCACTGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGcca                                                                                                                                                                           <  2:393942/251‑1 (MQ=255)                                                                          tCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATcggcgggcgcgggccggtcccccccacc                                                                                                                                                                        >  1:121908/1‑229 (MQ=255)                                                                                           ggCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTcc                                                                                                                                                     <  2:249920/251‑1 (MQ=255)                                                                                                 aTCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGt                                                                                                                                               <  1:80110/251‑1 (MQ=255)                                                                                                         ttGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCCCCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCa                                                                                                                                        <  1:143375/250‑1 (MQ=255)                                                                                                                  ggTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCt                                                                                                                               <  1:414772/250‑1 (MQ=255)                                                                                                                            aaCTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGt                                                                                                                     <  1:228685/250‑1 (MQ=255)                                                                                                                              cTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGGCGTGCGCCCCCCTCTTGGGGGGCATTTTTTCCCCCGTGCTTGCCATCACAGCTCGCaaaaaggttt                                                                                                                  >  1:120608/1‑247 (MQ=255)                                                                                                                                tAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTc                                                                                                                  <  1:365903/249‑1 (MQ=255)                                                                                                                                     gTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAAt                                                                                                             <  2:291386/249‑1 (MQ=255)                                                                                                                                              tcgttcTTTGGTGTATTGCTGGGGTGGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCTGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCATCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGGTCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑Ca                                                                                                  <  2:211433/250‑1 (MQ=255)                                                                                                                                               cgttcTTTGGTGTCTCGCTGTGGTTGTTCGATTATCTTGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGTGTTTCAAAACCGTGCTTTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACg                                                                                                <  2:302672/251‑1 (MQ=255)                                                                                                                                                              cGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAAt                                                                                 <  2:302638/251‑1 (MQ=255)                                                                                                                                                                          tCGATTAGCTGGCGGACCTTATTAAGATTGATAACTCCCCCGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GTGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGTCATTCATATTCCGGc                                                                      <  2:295472/250‑1 (MQ=255)                                                                                                                                                                             gttATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTCGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTT‑CC‑GGGCTGTGCTCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAATAAGTTTCCCAATCAGATGTT‑CACGGTCTGTGGCGGCAATCTTATTCCGTCGAtt                                                                  <  2:148192/250‑1 (MQ=255)                                                                                                                                                                                                           aCTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGGG‑CCCGGGTTG‑GCGCCCCCCTCTTGGGGGTGATTTTTCCCGCCGTGCTTGGAATCAAAGGTCGCAAAAAGTTTTCCCATCAGATGTT‑TACGGGCTGGGGGGGCCAACTTTTTTCGGGGATTTTCCTTTTCTTTGGGATAACCGTGGttttt                                    >  1:32512/1‑250 (MQ=255)                                                                                                                                                                                                                   ggcgggcgTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGtt                            <  2:120608/251‑1 (MQ=255)                                                                                                                                                                                                                    gcgggcgTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTc                           <  2:121908/251‑1 (MQ=255)                                                                                                                                                                                                                          gTTTCAAAACCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAAc                     <  2:109281/251‑1 (MQ=255)                                                                                                                                                                                                                               aaaaCCGTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGGTGTGCGCCACCATTTTGGCGGGCATTATTTCCGCCGGGCCTGCAATCAAAGCTCGCAAGAAGTTTCCCCATCAGATGTTTCA‑GGGCTGGGGGGGCCATCTTATTTCGGCGGTTTTTTTTCTCTTTTGTATTAACGGGGTTTTGTGCTGGTTCGGCAACGTCtt                >  2:60152/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTGCTGTTAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCCGTGCGCCCCCATCTGGGCGGTCATTATTTCCGCAGGGCTTGCAATCCAAGCTCGCAAGAAGTTTCCCCATCAGATGTT‑CCCGGGCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCCTTTGTGTAAACGTGATTTTGTGCTGGTTCGTGAACGTGTTTTACgt          >  1:393627/1‑251 (MQ=255)                                                                                                                                                                                                                                           ttAACCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTT‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCCGTGCTTGCAATCAAAGCTCGCAAGAAGTTTTCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGGTTGTCATTCTCTTTGGTATAACCGTGATTTTTTGCTGGGTCGGCCACGTCTTTAACGTGtttcc    >  1:96408/1‑248 (MQ=255)                                                                                                                                                                                                                                             aaCCTTCCTACCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTa  <  1:60152/251‑1 (MQ=255)                                                                                                                                                                                                                                                cTTCCTACCACCTGCGTTGTTGTATCTGATCTGCCCGAACGGGTTTATTTACGGGATCGGCGGTG‑CC‑GTGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAAt                                                                                                                                      <  1:245582/117‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     AATATTCCGCGAGAAAGCTTTAAGGCGATTATCTCCTCAGGCGGCAACGTTGATTCGCTGGTGAAATCGTTCCTCGGCACCAAACAGCACGGCATTATCGAGTTTTGCCTGCTGGTGTTCTCTAACTTAGCTGTTGCCAGTTCGTTCTTTGGTGTCACGCTGGGGTTGTTCGATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGTTTCAAAACCGTGCTGTTAACCTTCCTGCCACCTGCGTTGTTGTATCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTG‑CC‑GGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTT‑CACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTA  >  NC_000913/3890986‑3891472

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