Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	81,498	C→T	G154D (GGT→GAT)	leuB ←	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	81,498	0	C	T	100.0%	33.1 / NA	14	G154D (GGT→GAT)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/10);  total (4/10)

CCGTGGCGATCTCGTTAACGATCTCCCGCCATAAAATAGAGGATTGCAGCACGTTGGCTTTATCGATCGACGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTGAATAATTTGAAGTGCTTACGCAGAGGCAGCAGCGCGCCGCGTTCTGGTTGCTGGTCTGGTGGTAAATGTTCCCACTTCGGGCCGCC  >  NC_000913/81316‑81730                                                                                                                                                                                       |                                                                                                                                                                                                                                         cgggggcaagctcGTTAACGTTCTCCCGCCATAAAATAGAGGATTGCAGCACGTTGGCTTTATCGATCGACGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCTGTCAGTTCGCGCACACACAGGATGTCGAAg                                                                                                                                                                      <  1:396962/241‑1 (MQ=255)                                                 gcaCGTTGGCTTTATCGATCGACGTAACTTTGTGGCGACGCTTGCGAGCAGATTCCAACGCGATGCGGGCGATCCGTTCGCTCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTc                                                                                                                             >  1:338401/1‑244 (MQ=255)                                                                    cgaCGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTGTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCGGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTAGCCggg                                                                                                    <  1:222153/250‑1 (MQ=255)                                                                               ngngGNNACGTNTGCAAGCAGAGTCAAACCACATGCNGGCGATACGTTCGATCTCAAAACGGTGATACNCNTCGGTATCAAAGGCTTTTTCATATTGCCCGCTATCTTCGCGTCCTTTTNGCGGACTGAANTAGNNGCCGNCNGNCAGTTCGCGCACACACAGGATGTCGAAGCNGTTTGCGGCAATGTCTGCACGCANNGGNNNGANNGCTTc                                                                                                                             <  2:338401/213‑1 (MQ=255)                                                                                                                           cGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAACTAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAgg                                                                                              >  1:102987/1‑201 (MQ=255)                                                                                                                                   nccNAAANGNNGNNACANCTCGGTATCAAAGGCTGTGTCNNATTGTCCGCNATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCANACACAGGATNNCGNAGCCGNTNNNGNNNATNTCTGCACGCAGCGGACAGANTGNNTCCAGCNCCTGATACAGTTTTGCCGGGCGCAgg                                                                                              <  2:102987/191‑1 (MQ=255)                                                                                                                                                  cacCTNGGTATNANNGGCTTTTTCATATTGTCCGNTATCTTCGCGGCCt                                                                                                                                                                                                                               >  2:178852/1‑49 (MQ=255)                                                                                                                                                  cacCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCt                                                                                                                                                                                                                               <  1:178852/49‑1 (MQ=255)                                                                                                                                                  cacCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTg                                                                                        <  1:243451/184‑1 (MQ=255)                                                                                                                                                    ccTCGGTATCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACNNACAGGann                                                                                                                                                                            >  2:21427/1‑96 (MQ=255)                                                                                                                                                           aTCAAAGGCTTTTTCATATTGTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGc                                                                                                                                                  <  1:374501/117‑1 (MQ=255)                                                                                                                                                                      tttCATATTGTCCTCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCTCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGTGCAGGTTGCTGAATAATTTGAAGGGCTTACGCAGATGCAGCATCGCTCCGCGTTCTGGTTGCTGGTCTGGTGGTAAATGTTCCCACTTCGGgccgac  <  1:420224/250‑3 (MQ=255)                                                                                                                                                                               gTCCGCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGACTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTGAATAAGTTGAAGGGCTTACGCAGAGGCAGCAGCGCGCCGCGTTCTGGTTGCTGGTCTGGTGGTAAATGTTCCCACTTCGGgccgcc  <  1:449729/241‑1 (MQ=255)                                                                                                                                                                                   gCTATCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCTGCAATGTCTGCACGCAGCGGan                                                                                                                                       <  1:358083/104‑2 (MQ=255)                                                                                                                                                                                       |                                                                                                                                                                                                                                         CCGTGGCGATCTCGTTAACGATCTCCCGCCATAAAATAGAGGATTGCAGCACGTTGGCTTTATCGATCGACGTCACTTTGTGGCGACGCTTGCGAGCAGATTCAAACGCGATGCGGGCGATACGTTCGATCTCAAAACGGTGATACACCTCGGTATCAAAGGCTTTTTCATATTGTCCGCTACCTTCGCGGCCTTTTGGCTGACCGAAATAGATGCCGCCGGTCAGTTCGCGCACACACAGGATGTCGAAGCCGTTTGCGGCAATGTCTGCACGCAGCGGACAGAATGCTTCCAGCCCCTGATACAGTTTTGCCGGGCGCAGGTTGCTGAATAATTTGAAGTGCTTACGCAGAGGCAGCAGCGCGCCGCGTTCTGGTTGCTGGTCTGGTGGTAAATGTTCCCACTTCGGGCCGCC  >  NC_000913/81316‑81730

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