Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,359,314	G→A	Q240* (CAG→TAG)	cadB ←	lysine:cadaverine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,314	0	G	A	100.0%	75.6 / NA	27	Q240* (CAG→TAG)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/11);  total (16/11)

TTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGT  >  NC_000913/4359071‑4359551                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              ttttCGGAATGCCGTGGCTGCTGACTTCACCAAAAACTTTCTGGAAGTTACCGTTGTAAGCGGCACTTCCACCTGCCTTGCCTCCCAACACCAGCCAGGAGCCTAGAGAAGTGAGGCACGCAAAGGCGGTGAATACAAAAATCAGCGGTGCAGCCCTGTTACCGATGATAGTTGAAACACTGATTGCAAATGAAACACCGAAAGCCGCCATTACAGAAGACGGATACATACCTGAAAGCACCTAAGTCGc                                                                                                                                                                                                                                         <  1:294167/250‑1 (MQ=255)                                                       cgtTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGCTAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGc                                                                                                                                                                                                                                      >  1:209910/1‑199 (MQ=255)                                                                                gCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGATACAAGCGGCGAAGCTCAGTTACAGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTAGCAGCGATGTAAACa                                                                                                                                                                                                                            >  1:294501/1‑184 (MQ=255)                                                                                     ncAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGAGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGt                                                                                                                                                     >  1:80007/2‑250 (MQ=255)                                                                                       aaCATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCCGCGGCGCAGCCCGGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAACACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                               <  1:432778/206‑1 (MQ=255)                                                                                                      cccAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACa                                                                                                                                              >  1:127496/1‑240 (MQ=255)                                                                                                                         cAAAGGCGGTGAAGGAAGAAACCAGCGGCGCAGCCCGGTTACCGAGGATAGTGGAAGCACTTATGCCAACCGGACCACCGGAAGCCGCCATGACAGAAGACGGATACATACCGGAAATCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACTTTTCGGGTTTTTAACCATACCCGTACTTACAGCTGCGGATTCACCCCCCACGAAGGCCcac                                                                                                                <  1:301921/251‑1 (MQ=255)                                                                                                                                  tGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGACGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAg                                                                                                                                                                                                                                       >  1:444259/1‑123 (MQ=255)                                                                                                                                       gCAGAACCCAGCGGCGCAGACCAGTTACCGAGGATAGTTGAGGCACGGATTGCAAACGGAGCCCCTGAAGCCGCCATTACAGAAGACGGATACATACCTGAAAGCACCTAAGTCGCAGCTATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                               <  1:304620/158‑1 (MQ=255)                                                                                                                                       gCAGAAACNANNNGNNCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGANCACCGGAAGCCGCCANTANAGCANACGNATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGt                                                                                                                                                                                               >  2:304620/1‑158 (MQ=255)                                                                                                                                                aGCGGCGCAGCCCGTTTACCGAGGATAGTTGAAGCACTGATTGCAAACGAAGCACCGGAATCCGCCATGACAGAAGTCGGATTCATACCTGAAAGCACCTAAGTCGCAGCGATGTAAACAATGCCGGCTAAACCAGTACCCAGCATGGTTGCCAGCAGAACGGTACGTTTCTGGTTGTTAACGATACCAGTACGGACGGCTGCGGATTCGACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                                                                         <  1:412680/251‑1 (MQ=255)                                                                                                                                                                               nnccNNNNNTTGNAAACGGCGCCCCGGACGCCGCCATTACAGAAGTCGGATACATACCTGAAAGCACCTAAGTCGCAGCGATGTAAACAATTCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTAGGTCTCGGGtgt                                                                                                                                                                   <  2:265365/143‑3 (MQ=255)                                                                                                                                                                                  cACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCGCCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACTGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTc                                                       <  1:414112/251‑1 (MQ=255)                                                                                                                                                                                          gCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTAc                                                                                                                                                   >  1:174363/1‑151 (MQ=255)                                                                                                                                                                                          gCAAAAGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTAccc                                                                                                                                                                                                       >  1:342903/1‑99 (MQ=255)                                                                                                                                                                                                     aCCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATCCCTGCTAAACAAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACCGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAt                                                         >  1:302193/1‑230 (MQ=255)                                                                                                                                                                                                      ccGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGTAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACGCCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGc                                   <  1:317459/251‑1 (MQ=255)                                                                                                                                                                                                           ggccgccATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCAt                              >  1:124030/2‑251 (MQ=255)                                                                                                                                                                                                             ncgccATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGCAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTAGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATc                             >  1:78660/2‑250 (MQ=255)                                                                                                                                                                                                               gccATTACAGAAGCCGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                                                                      >  1:94767/1‑175 (MQ=255)                                                                                                                                                                                                               gccATTACAGAAGACGGATACNNACCGGAAAGNACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTNGCCAGCGGANNGGNACGTTNCNNNTNNNTANCCATACCAGTACTTACAGCNGCNNATTCAANACCCACGAAGGCCCACATgcagagca                                                                                                      <  2:94767/175‑1 (MQ=255)                                                                                                                                                                                                               gccATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATCCCTGCTAAACCAGTACCAAGCATGGTTGCCAGCGGAACGGAACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCACACAGGCAGAGCAGACTACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGATAGCTGCATAAATTGCCGCATCaa                           >  1:302791/1‑250 (MQ=255)                                                                                                                                                                                                                     aCAGANNACNNNTANNTACCGGAAAGCACCTAAGTCGCANCGATGTAAACAATACCTGCTAAACCAGTANCNANCATGNNTGCNAGCGGAACGGNACGTTTc                                                                                                                                                                         >  2:333978/1‑102 (MQ=255)                                                                                                                                                                                                                     aCAGAAGACGGTTACAGACCCGAAAGCACCTAAGTCGCAGCGATGTACACCATACCTGCTAATCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTc                                                                                                                                                                         <  1:333978/102‑1 (MQ=255)                                                                                                                                                                                                                      cAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGGACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACg                                                                                                                         <  1:240638/149‑1 (MQ=255)                                                                                                                                                                                                                                       ccGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGt  >  1:308515/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGAGGAACGGTAAGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGt  >  1:308506/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              TTTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGT  >  NC_000913/4359071‑4359551

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