Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	81,102	G→A	S286L (TCG→TTG)	leuB ←	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	81,102	0	G	A	100.0%	59.3 / NA	22	S286L (TCG→TTG)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/15);  total (7/15)

AGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATT‑TCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCG  >  NC_000913/80859‑81318                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         aGCCATGATTACACCCCGTCTGCTCCATAGCGGGCAAGTAAACTGCCCATTGTC‑GCGGTACGAACGGCGGTAGCGCCACGGGCTAAATCCCCGGTGCGAATGCTTTCTTCTAATGCGCGGTGAAT‑GGCGCTTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCTGCAGTGCCAGCGAAAGGATGTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAgccgcc                                                                                                                                                                                                                   <  2:378303/251‑1 (MQ=255)  gCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATT‑TCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGccgccc                                                                                                                                                                                                                  <  1:369437/251‑1 (MQ=255)  cccatcttttccccccTTCTGCGCCATAGCGGTCAATGATATCCCCCATT‑TCACTGTTACTAACGGCGGCAGCTCCGCTGCCTAAATTCCCGGTGCGAATTCCTTCTCCTACTGTGCGGTTAAG‑GGCGGGTTCAAGTGCGCAAGCCGCATCCTCGGCATGCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTCGAATCGGGTTGGCGAGGTTTTTGCCTGCGATTTCTGGTGCCAAGcctccc                                                                                                                                                                                                                  <  2:93150/239‑1 (MQ=255)               cccTTCTGCTACATAGCGGGCAATGATATAGCCCATT‑TACTCGGTACTAACGGCGGCAACGCCCCGGGCTAAATCCCCGGTGCGAATGCCTTCTTATAATGCGCGGTTACT‑GGCACGTTCAATAGCGCAAGCCGCATCATCGGCCTCCACACTGTAACGCAGCAGAAGAGCCAGCGAAAGGATAGTTGCACACGGGTTGTCGAAGTTTTTGCCTGCAACATCTGGGTCCAAGCCGCCGCCCGGGTCATAc                                                                                                                                                                                                      >  2:283618/1‑250 (MQ=255)                        tACATAGCGGGCAATGATATCGCCCATT‑TCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAAc                                                                                                                                                                                            >  1:251814/1‑251 (MQ=255)                                 ggCAATGATACCGCCCATT‑TCATCGGTACTAACGGCGGCAGCGCCACGGGCTAATTCCCCGGGGCGAATGCCTTCTTCTAAGGCGCGGTTAAT‑GGCGGGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGt                                                                                                                                                                                   <  1:104/251‑1 (MQ=255)                                               ccATT‑TCCTCGTCACTATCGGCGGCAGCGCCTCGGGCTAAATCCCCGGTGCGAATGCCTCCTTCTAATGCGCGGTTAAG‑GGCGCGTACAATGGTGAAAGCCGCATCATCGGCATCCAGGTTGTAACGCAGCGGCAGTGCCAGCGAAAGGATTTGTGAAACCGGGTTGGAGATGTTTTTGCCTGCGATATCGGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGaa                                                                                                                                                                     <  1:371588/251‑1 (MQ=255)                                                  tt‑tCACCGGTACAACCGGCGGAAGCGCCACGGGCTAAATCCCCGGGGCGAATGCTTTCTTCAAATGCGCGGTTAAG‑GGCGCGTTAAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAgg                                                                                                                                                                   <  1:203020/250‑1 (MQ=255)                                                                aaCGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCGGTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGGGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGGTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATcccccacg                                                                                                                                                     >  1:187564/1‑247 (MQ=255)                                                                                              ggTGCGAATGCTTTCTTCTAAGCCGCG‑TTAATTGGCGCGTCCAATGGCGCAAGCCGCATCAGCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATTTTTTTGCCTGCGATATCGGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTcagaca                                                                                                                       <  1:95375/251‑1 (MQ=255)                                                                                              ggTGCAAATGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCTTACAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTACCGCAGCAGCAGTGCCGGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCAAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTcagaca                                                                                                                       <  1:231844/251‑1 (MQ=255)                                                                                                      tGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCg                                                                                                               >  1:91401/1‑251 (MQ=255)                                                                                                              ttgaagtcggtgTTAAT‑GGCGCGTTCAATGGCGCAAGCCGCATCATCGGCAGCCGGTCTGTAACGCAGCAGCGGGGCCAGCGAAAGGATTTGGGCAATGGGGTGGGCGATGTTTTTGCCTGGGATATCTGGTGCTAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCaaa                                                                                                          <  1:45855/237‑1 (MQ=255)                                                                                                                    gcgcGGTTAAT‑GGCGCGTTAAAGGGCGCAAGCCGCATCAGCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGa                                                                                                 <  1:142484/251‑1 (MQ=255)                                                                                                                                  ggcgtcaagtGGCCCAAGCCGCACCACGGGCATCCAGGCTGTAACGCAGCGCCAGTGCCGGCGAAAGGATTTGGGCAACCGTGTTGGCGATTTTTTTGCCGGCGATATCTGGGGCCAGGCCTCCCGCCGGTTCATACAGGCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACg                                                                                    <  1:163496/242‑1 (MQ=255)                                                                                                                                     gTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTCCAGGCTGGCGGCAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTcc                                                                                 >  2:163830/1‑250 (MQ=255)                                                                                                                                                   gCCGCATCACCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTCCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATc                                                                   <  1:163830/251‑1 (MQ=255)                                                                                                                                                          cagcggattCCAGGCTGAACCGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATTTCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATc                                                            <  1:177889/243‑1 (MQ=255)                                                                                                                                                                  cccAGGCTGTAACGCAGCAGCATTGCCAGCGAAAGGTTTTGTGCAATCGGGTTGGCGATTTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCAt                                                    <  1:104795/250‑1 (MQ=255)                                                                                                                                                                                                 aaaGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGt                     >  1:73299/1‑251 (MQ=255)                                                                                                                                                                                                              aaTCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGt        >  1:97938/1‑251 (MQ=255)                                                                                                                                                                                                                    gtgggCGTTGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCg  <  1:220482/248‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                         AGCCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATT‑TCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAAT‑GGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGCAGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCG  >  NC_000913/80859‑81318

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