Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,338,465	G→A	P686L (CCA→CTA)	adiA ←	arginine decarboxylase, degradative

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,338,465	0	G	A	100.0%	50.2 / NA	19	P686L (CCA→CTA)	adiA	arginine decarboxylase, degradative
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/9);  total (10/9)

CACGCACATAACGTGGTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTGGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGCCAGGCAAACATGGTGTCACCCAGATCGTGAATCCCCATGTTCGCGTAAGTGTCAGGATATTGTTCAACAAGTTCCGGCATCACCTGCGCCAGCGGTG  >  NC_000913/4338263‑4338696                                                                                                                                                                                                            |                                                                                                                                                                                                                                        cACGCACATAACGTGGTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCGCCAGACAGCAGCATCGGGATTCCTGGCGGATACAGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATTCCAGTTCGACAttgttn                                                                                                                                                                                                      <  1:39155/238‑2 (MQ=255)      acaTAACGTGGTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGGTTTCACCAGACAGCAGCATCGGGATTCCTGGTGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAAGGGATTCCAGTTCGACATTGTTGTCGACAATCGCGTTGTAc                                                                                                                                                                                    <  1:274228/251‑1 (MQ=255)      acaTAACGTGGTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACGTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTTAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAAGGCGACAGTGTTGTCGACAATCGCGTTGTAc                                                                                                                                                                                    <  1:117565/251‑1 (MQ=255)               ggTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCCcgc                                                                                                                                                                           >  1:329282/1‑251 (MQ=255)                gTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGCGGCCCCAGGCTTGCAGCGAGCGTCAAACAA‑TTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCCTCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGCTGCGTCCTAGCAGAt                                                                                                                                                                                                                                   >  1:160395/1‑194 (MQ=255)                       ccGTCAATAATTTCAGTCCCTTCAGTTGCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCAGCGGGATGCCTGGCGGATACGGGATAACTGAGTTTCCCGCGATGCGTCCTAGCAGATTTTCAATGGAGACCAGTTCGACATTGTTGTCGACAATCTCGTTGTACGCTTCACGCGGGGTGAc                                                                                                                                                                   <  1:283415/251‑1 (MQ=255)                                                                          tggtCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAACAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATa                                                                                                                                                                                                                       >  1:131300/1‑148 (MQ=255)                                                                          tggtCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCAGAAGTTTTCACCAGAAAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCAAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCgg                                                                                                                                                                         >  1:170881/1‑194 (MQ=255)                                                                              cccAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGCTAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGCGCGCCAGGGTTGttt                                                                                                            >  1:73267/1‑251 (MQ=255)                                                                                       gCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATGTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACTCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTGTCTTTCAGc                                                                                                   <  1:169944/251‑1 (MQ=255)                                                                                           cgagcgT‑AAATNACTTACTTGCGGACTGTTTTTATCGCCGAAGNTNTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGNTTGCCGCGATNCGTNNTAGCNGNTNTTCAATGGATACCAGTTCGACATTGTTGTNGACAATCGCGTTGTACGCTTCACNNGGNNNGANNTCCGc                                                                                                                                                             <  2:236171/189‑1 (MQ=255)                                                                                           cgagcgT‑AAATAACTTACTTTCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGTCATTGTTGTCGACAATCGCGTTGTACGCTTCATGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGCCAgg                                                                                               <  1:284107/251‑1 (MQ=255)                                                                                           cgagcgT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATCCGGGATAACTGAGTTTGCCGCGCTGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGc                                                                                                                                                             >  1:236171/1‑189 (MQ=255)                                                                                            gagcgT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCgg                                                                                                                                               >  1:161961/1‑202 (MQ=255)                                                                                                                                             ccAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGCCAGGCAAACATGGTGTCACCCAGATCGTGAATCCCCATGTTCGCGTAAGTGTc                                              >  1:311671/1‑251 (MQ=255)                                                                                                                                                                       cGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGGTCGACATGGTTGTCGACAATCGCGTTGTACGCTGCACGCGGGGTGACTTCCGCTACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTtttctttc                                                                                                      <  1:161945/169‑1 (MQ=255)                                                                                                                                                                               ggATAACTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGc                                                                                                   >  1:162194/1‑164 (MQ=255)                                                                                                                                                                                     cTGAGTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGc                                                                                                   <  1:7975/158‑1 (MQ=255)                                                                                                                                                                                         gTTTGCCGCGATGCGTCCTAGCAGATTTTCAATGGATACCCGTTCGACATTGTTGTCGACAATCGCGCTGTAAGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGCCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGCCAGGCAAACATGGTGTCACCCAGATCGTGCATCCCCATGTTCGCGTAAGTGTCAGGATATTGTTCACCAAGTTCCGGCATCACCTGCGCCCGCGgtg  >  1:177521/1‑251 (MQ=255)                                                                                                                                                                                                            |                                                                                                                                                                                                                                        CACGCACATAACGTGGTAAATACCGTCAATAATTTCAGTCCCTTCAGTTTCGTGTTCAAATCCAGGGAAATGGTGGTCCCAGGATTGCAGCGAGCGT‑AAATAACTTACTTGCGGACTGTTTTTATCGCCGAAGTTTTCACCAGACAGCAGCATCGGGATTCCTGGCGGATACGGGATAACTGAGTTTGCCGCGATGCGTCCTGGCAGATTTTCAATGGATACCAGTTCGACATTGTTGTCGACAATCGCGTTGTACGCTTCACGCGGGGTGACTTCCGCCACCGGCAGGCCGGAATAGGCTTCGTTCAACCGTGCGCCAGGGTTGTTTTCTTTCAGCCAGGCAAACATGGTGTCACCCAGATCGTGAATCCCCATGTTCGCGTAAGTGTCAGGATATTGTTCAACAAGTTCCGGCATCACCTGCGCCAGCGGTG  >  NC_000913/4338263‑4338696

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