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%	90.6 / NA	27	Q240* (CAG→TAG)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/17);  total (10/17)

TCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGC  >  NC_000913/4359074‑4359548                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           tCGGAATACCGTTGCTGTCGACTTCACCATAAACTTGCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCACCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCATTTACCGAGGATAGTTGAAGCACGGATTTCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGc                                                                                                                                                                                                                                   <  1:246684/250‑1 (MQ=255)                tGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGACAGCACCTAAGTCGCAGCGATGTAAACAATACCt                                                                                                                                                                                                                   >  2:147178/1‑251 (MQ=255)                               aaaCTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACATGCCTCACCTACCAACATCATCCAGGAGCCCAGAGAAATCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAGGCCGCCATTACAGAAAGCGGCTACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGtcacc                                                                                                                                                                                                    >  1:40184/1‑247 (MQ=255)                                           aGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGcc                                                                                                                                                                                        >  1:187200/1‑251 (MQ=255)                                               cccGTCGTTAGCGGCACGTACACCTTCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGTAGCCCCGTAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGGGATGTAAACAATACCGGCTAGACCAGTAGCCAGCAGGGTTGGCAGCg                                                                                                                                                                                    <  2:35260/250‑1 (MQ=255)                                                                                     aCATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGGGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTGGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACtt                                                                                                                                              <  1:14033/251‑1 (MQ=255)                                                                                              aGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCg                                                                                                                                     <  2:207960/251‑1 (MQ=255)                                                                                                   cccAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTc                                                                                                                                <  1:96462/251‑1 (MQ=255)                                                                                                     cAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCACAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACAAGTACTTACAGCTGCGGATTCa                                                                                                                               >  2:34158/1‑250 (MQ=255)                                                                                                       gagaAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAac                                                                                                                             >  2:10136/1‑250 (MQ=255)                                                                                                                       aaaGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGAATAGTGGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCcaca                                                                                                            <  2:106400/251‑1 (MQ=255)                                                                                                                        aaGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAg                                                                                                           <  1:215753/251‑1 (MQ=255)                                                                                                                           gCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGgca                                                                                                        <  1:24125/251‑1 (MQ=255)                                                                                                                                            cAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTa                                                                                        >  1:48857/1‑250 (MQ=255)                                                                                                                                                     aGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTTCCCAGCAGGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACTAATGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCAt                                                                              <  1:139363/251‑1 (MQ=255)                                                                                                                                                                 gAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTgg                                                                  >  1:118604/1‑251 (MQ=255)                                                                                                                                                                    gATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTAGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACATGCAGCGCAGAATACTTTTAATGATCGCATGACCATCAGTggaat                                                               >  2:185233/1‑248 (MQ=255)                                                                                                                                                                               cACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTc                                                    <  1:233042/251‑1 (MQ=255)                                                                                                                                                                                ccTGATTGCAAACGGCGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCAGGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTTGTATCCGCAGTATTcc                                                   <  2:217718/250‑1 (MQ=255)                                                                                                                                                                                   gATTGAAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCa                                                  <  2:46795/249‑1 (MQ=255)                                                                                                                                                                                      tGCAAAAGGAGCCCCGGAAGCAGCTATTACAGAAGACGGATACATACCGAAAAGCACATAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGtt                                               >  2:213614/1‑249 (MQ=255)                                                                                                                                                                                                            gccATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGc                                                           <  1:34507/215‑1 (MQ=255)                                                                                                                                                                                                            gccATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGc                                                           >  2:34507/1‑215 (MQ=255)                                                                                                                                                                                                             ccATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGCACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAAc                      <  1:163207/251‑1 (MQ=255)                                                                                                                                                                                                               aTTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCAGACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCa                    <  1:147388/251‑1 (MQ=255)                                                                                                                                                                                                                  aCAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAAt                  <  2:63087/250‑1 (MQ=255)                                                                                                                                                                                                                                 atacCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGc  <  1:99221/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           TCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGC  >  NC_000913/4359074‑4359548

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: ‑