Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,359,909	C→T	W41* (TGG→TGA)	cadB ←	lysine:cadaverine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,909	0	C	T	100.0%	52.9 / NA	19	W41* (TGG→TGA)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/9);  total (10/9)

CGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAA  >  NC_000913/4359669‑4360144                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            cgccgGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCaa                                                                                                                                                                                                                                   <  1:248523/251‑1 (MQ=255)             aTCATTTAATACTGCGAAGAAGGCGGAACGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCACTCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAACTGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACATCAGATAGCAATACCATCGATAc                                                                                                                                                                                                                       >  2:254105/1‑251 (MQ=255)                        ctgggaagatggaggaccgctaagataCCGCGGTCATACCAATCGGCAGGTTCCCACTCCAGATACAATAGTCATAAAGAACCCCTGTCTGAAAACCAAATGCAGGGGACATTTCTCCGGCATAAACACTTGGGCCACCTTGTTGCGGGTTTCTTGTTGCCAGTCGGGCATATACATACGCCAGCGCCATTGCACCAATACTAGAGATAATCCAACCTAAGCTACCAATACCACCGCTACTTGCTAGGTTc                                                                                                                                                                                                            >  2:372147/21‑251 (MQ=255)                          gggAAGAAGGTTGAAAGATAAGATACAGCGGGAATACCAATCGACAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCGGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCGCCTTGTTGCGTGTTTTTATTTGCCAGTCGGGCATATACATACGCCATCTACATTGCACCAATAATAGAGATAATCCAACCTCAGCTAGCAATACCCCCGATTCttccgatgttcg                                                                                                                                                                                                           <  2:382232/250‑11 (MQ=255)                                           ataagataCAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGTATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCCCCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCCCCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTATGTTCGCAGGTAATAATGCAATAc                                                                                                                                                                                         <  1:102549/251‑1 (MQ=255)                                                  aCAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGAAATAATGCAATACCGCTccc                                                                                                                                                                                  >  1:101547/1‑251 (MQ=255)                                                     gCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcac                                                                                                                                                                               >  1:7505/1‑250 (MQ=255)                                                        gTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatcat                                                                                                                                                                            >  1:339973/1‑251 (MQ=255)                                                                    gCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGcaac                                                                                                                                                                >  1:231217/1‑251 (MQ=255)                                                                                                                                     aaTTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTAGCCAGTCGGGCATATACATACGCCAGCGCCATTGAACCAATAATCGAGATAATCCAACCTAAGATAGAAATACCACCGATACTTGCTAGGTTCGCAGGCACTAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTCCAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCAtt                                                                                               >  1:398038/1‑251 (MQ=255)                                                                                                                                                       aagtGGGCCACCTTGGTGCGGTTTTTTTGTTGCCAGTCGGGCATATACAGACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATGCCACCGATACTTGCTAGGTTCGCATGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAacac                                                                                                                                                           <  1:262027/170‑1 (MQ=255)                                                                                                                                                       aaTTGAGCCACTTTGTTGCGGGTTTTTTGGTGCCAGTCTGGCATATACACACGCCAGCGACCTTGAACCCCTAATTGAGATAATCCAACCTAAGAAAGCAATACCACCGATACTTGCTAGGATCGCAGGCAATAATGAAATACCGCTCCCCATCATATTACCGGCAACAacac                                                                                                                                                           >  2:262027/1‑173 (MQ=255)                                                                                                                                                          tGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACGTACGCCAGCGACATTGCACCAATAATAGAGATAACCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGt                                                                           <  2:272217/250‑1 (MQ=255)                                                                                                                                                          tGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGt                                                                           <  2:98126/250‑1 (MQ=255)                                                                                                                                                              cctccTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCatga                                                                      <  2:99405/248‑1 (MQ=255)                                                                                                                                                                                   gttgCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGGCCGGGTCATGATGTATAACTATGTCCTgacca                                                 <  1:337054/251‑1 (MQ=255)                                                                                                                                                                                                   atacatacGCCAGCGACATTGCAACAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGAt                                  >  1:115033/1‑250 (MQ=255)                                                                                                                                                                                                               gCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACg                     <  2:84030/251‑1 (MQ=255)                                                                                                                                                                                                                                    agaTAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTAGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTAATGATGTATAACTATTTCCTGACCAGACAAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTaa  >  1:53385/1‑249 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            CGCCGGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAA  >  NC_000913/4359669‑4360144

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