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

TTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCA  >  NC_000913/4359072‑4359559                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      tttCGGAATACCGTTGCTGCCGCCTTCACCAAAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACAGCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAGGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCa                                                                                                                                                                                                                                                <  2:28037/250‑1 (MQ=255)            cGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTCAACCTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGAACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCCGCCCAGTTACCGAGGATAGTTGAGCCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAAc                                                                                                                                                                                                                                     >  2:244278/1‑250 (MQ=255)                           cACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAAcc                                                                                                                                                                                                                      <  1:127632/250‑1 (MQ=255)                           cACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAAcc                                                                                                                                                                                                                      <  1:4266/250‑1 (MQ=255)                                  aaCTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCAACGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTAccc                                                                                                                                                                                                               >  1:183916/1‑250 (MQ=255)                                              gttaaagtcgtTAGCGGCACGTACACCGGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGcc                                                                                                                                                                                                   >  2:64344/7‑250 (MQ=255)                                                                  gTACACCTGCCTGTCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCTGAACGGTACGTTTCggg                                                                                                                                                                              <  1:236560/251‑1 (MQ=255)                                                                       cctgcctgGCCTGCCAACATCAGCCAGGAGCCCAGAGAAGTCAGGAACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACGGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGttttt                                                                                                                                                                         <  2:213309/251‑1 (MQ=255)                                                                                  tACCAACATCATCCAGGAGACCAGAGAAGTCAGGCAAGCAAAGGCGGAGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGAGAGATGACGCACAGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCAAAGGCGCAGCGATGAAAACAATACATGCTAAACAAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCAGGTTTTTAACCATAcc                                                                                                                                                                >  2:230705/1‑249 (MQ=255)                                                                                          tcctcccggcaccgagagAAGTCAGGCCCGCAACCGCGGTGAATGCAGAAACCAGCGGCGCAGCCCCGTTACCGAGGATAGTTGAAGCAATGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGCTACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTAc                                                                                                                                                       >  2:47065/15‑250 (MQ=255)                                                                                                aGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCg                                                                                                                                                >  1:53449/1‑251 (MQ=255)                                                                                                    gCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGAtt                                                                                                                                            >  1:202338/1‑251 (MQ=255)                                                                                                       cAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCaa                                                                                                                                         <  1:84827/251‑1 (MQ=255)                                                                                                         gagaAGTCAGGCACGCAAAGGCGGTGAATGCAGAACCCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAaca                                                                                                                                       <  2:25183/251‑1 (MQ=255)                                                                                                                                                    cgcAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGt                                                                                                                                                                                        <  1:40334/159‑1 (MQ=255)                                                                                                                                                    cgcAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGt                                                                                                                                                                                        >  2:40334/1‑159 (MQ=255)                                                                                                                                                          ccAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGTGCACCGGAAGCCGCCATTACAGAAGACGGCTACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTACACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGATGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGAc                                                                                      >  1:140130/1‑251 (MQ=255)                                                                                                                                                                aCCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGAACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                                                                                  >  1:225830/1‑249 (MQ=255)                                                                                                                                                                         aGTTAAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCg                                                                       <  1:219159/251‑1 (MQ=255)                                                                                                                                                                                gCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACAAACCGGAAAGCACCTAAGTCGAAGCGATGTAAACAATACCTGCTAAACCAGAACCCAGCATGGTTGCCAGCGGAACGGTAAGTTTCGGGTTTTTAACCATACCAGTACTTACAGCAGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGccgaatt                                                                >  2:271736/1‑245 (MQ=255)                                                                                                                                                                                 cACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTc                                                               >  1:198415/1‑251 (MQ=255)                                                                                                                                                                                          cAAACGGAGCACCGGCAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCt                                                      >  1:78520/1‑251 (MQ=255)                                                                                                                                                                                          aaaaCGGAGCCCCGGAGGCCGCCATTACATAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCt                                                      <  2:99637/250‑1 (MQ=255)                                                                                                                                                                                                                   tACAGAAGACGGCTACATACCGGAAAGCACCTACGTCGCAGCGATGTAAACAATACCTGCTAAACCCGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGACGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGAGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAACCCAAt                             >  2:90499/1‑251 (MQ=255)                                                                                                                                                                                                                     cAGAAGACGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATgc                           >  1:226746/1‑251 (MQ=255)                                                                                                                                                                                                                        aagaCGGATACATACCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATgcc                          >  1:5828/1‑249 (MQ=255)                                                                                                                                                                                                                                    tacCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGc             >  2:105460/1‑250 (MQ=255)                                                                                                                                                                                                                                    tacCGGAAAGCACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGc             <  1:105460/250‑1 (MQ=255)                                                                                                                                                                                                                                              cACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATcacca  <  1:165334/251‑1 (MQ=255)                                                                                                                                                                                                                                              cACCTAAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGGTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATcacca  <  2:146168/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      TTTCGGAATACCGTTGCTGTCGACTTCACCATAAACTTTCGGGAAGTTACCGTCGTTAGCGGCACGTACACCTGCCTGGCCTACCAACATCATCCAGGAGCCCAGAGAAGTCAGGCACGCAAAGGCGGTGAATGCAGAAACCAGCGGCGCAGCCCAGTTACCGAGGATAGTTGAAGCACTGATTGCAAACGGAGCACCGGAAGCCGCCATTACAGAAGACGGATACATACCGGAAAGCACCTGAGTCGCAGCGATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAGCTGCATAAGTTGCCGCATCAAACCAATGCCAGCCAACAATAGCAGTCATCACCA  >  NC_000913/4359072‑4359559

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