Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,026	C→T	G64S (GGT→AGT)	paeA ←	putative cadaverine/putrescine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,026	0	C	T	100.0%	30.0 / NA	12	G64S (GGT→AGT)	paeA	putative cadaverine/putrescine exporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (6/6);  total (6/6)

GTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCCACTGCGTTCAGACCGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTTGATTTTGCGTGAGGCGGCAAGCGAGATCTCGGACATCGAGAAGAACGCACTTACAGCGATCAAGCAGAGTATGACTAAAATACTGTTTAACATATCTTATCCGACCTTTCAGGCCAGATCCTCGGAAGGGGAAGTTGATTAATTTG  >  NC_000913/4440781‑4441268                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   gTACAGGCAGAAGCGCATCGGTTTGATGATACGCAAAGCCACAGTTTCTGGCGCAATCATACCGATGCTTTTCGTGGTTAAATCCCCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTTCTCAGAGAGCTCTGCCTACATGTATCGGGAGAACATGCTTGGAAAAGCTCGAGAAAATGCCCCATCACCGACGCTACCGCCGAGAATCGCCACTGCGTTCAGACTGAttt                                                                                                                                                                                                                                               <  1:165180/251‑1 (MQ=255)                                              aggacgtgtgccctcccgatctcGTTTCGGGTTTAAATCCGCAAACAGGAGAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCTCATCACCGACGATACCGCCGAGAATCGCCACTTCGTTCAGTCTGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGaa                                                                                                                                                                                                  <  1:170100/229‑1 (MQ=255)                                                                               tAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGTGCTCTGCCGACATATAGCGGGAGAACAGGCTCTGAAAAGCTGGAGAAAATGCCGCATCACCGGCGATACCGCCGAGAATCGCCACTGCGTTCAGACTGATTTGGACCACGGTAAAGAACATGGCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATcagcc                                                                                                                                                                  <  1:188019/250‑1 (MQ=255)                                                                                             ggATAAACATGCCAGTCNCNAACAAGAAAGAGAGAAAAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACNTATAGCTTGANAACNNTCTTNANANAGCTGGAGAAAATGCCGCATCTCCGACGANACCGCCGAGAATCGCCACTGCGTNNAGNNNGANNTGGAc                                                                                                                                                                                                                                           <  2:240852/163‑1 (MQ=255)                                                                                             ggATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCCACTGCGTTCAGACTGATTTGGAc                                                                                                                                                                                                                                           >  1:240852/1‑163 (MQ=255)                                                                                                                                                                                  gCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCCACTGCGTTCAGACTGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTTGATTTTg                                                                                                                                             >  1:190669/1‑172 (MQ=255)                                                                                                                                                                                               ggAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCCACTGCGTTCAGACTGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTTGATTTTGCGTGAGGCGGCAAGCGAGATCTCGGACATCGAGAAGAACGCACTTACAGCGATCAAGCAGAGTATGACTAAAATACTGTTTAACATATCtt                                                  >  1:149581/1‑250 (MQ=255)                                                                                                                                                                                                         ccGCATCACCGACGATACCGCCGAGAATCGCCACTGCGTTCAGACTGATTTGGACCACGGTAAAGACCATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTTGATTTTGCGTGAGGCGGCAAGCGAGATCTCGGACATCGAGAAGAACGCACTTACAGCGATCAAGCAGAGTATGAc                                                                         >  1:163653/1‑217 (MQ=255)                                                                                                                                                                                                                          ccgccgAGAATCGCCTCTGCGTTCAGACTGATTTGGACCACGGTAAAGAACATGCCGGGATTTGCCTGCATATTCAGAACGCGTTGTGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTGGATTTTGCGTGAGGCTGCAAGCGGGATCTCTGAGATCGAGAAGAACGCACTGACATCGATCAAGCATAGTATGACTAAAATACTGTTTAACATATCTTATccca                                            <  1:308369/229‑3 (MQ=255)                                                                                                                                                                                                                            gccgAGAATCGCCACTGCGTTCAGACTGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATcagcn                                                                                                                                                                  <  1:203918/109‑2 (MQ=255)                                                                                                                                                                                                                            gccgAGAATCGCCACTGCGTTCAGACTGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATcagcc                                                                                                                                                                  >  2:203918/1‑109 (MQ=255)                                                                                                                                                                                                                                      gCCACTGCGTTCAGACTGATTTGGACCACTGTAAAGAACATGCCTGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATAGTGCCGTCATCAGCCAGCAGTTTAAGGTTGATTTTGCGTGAGTCTGCAAGTGAGATCTCGGACAGCTAGAAGAACGCCCTTACAGCGATCAAGGAGATTATGACGAAAATACTGTTTAACATTTCTTATCCGACCTTTCAGGCCAGATCCTCGGACGGGGAAGTTGa          <  1:159938/251‑1 (MQ=255)                                                                                                                                                                                                                                              gTTCAGACTGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTTGATTTTGCGTGAGGCGGCAAGCGAGATCTCGGACATCGAGAAGAACGCACTTACAGCGATCAAGCAGAGTATGACTAAAATACTGTTTAACATATCTTATCCGACCTTTCAGGCCAGATCCTCGGAAGGGGAAGTTGATTAATTtg  >  1:158892/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   GTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCCACTGCGTTCAGACCGATTTGGACCACGGTAAAGAACATGCCGGGATTTTCCTGCATATTCAGAACGCGTTGGGCATTTATATTGCCTTCATCAGCCAGCAGTTTAAGTTTGATTTTGCGTGAGGCGGCAAGCGAGATCTCGGACATCGAGAAGAACGCACTTACAGCGATCAAGCAGAGTATGACTAAAATACTGTTTAACATATCTTATCCGACCTTTCAGGCCAGATCCTCGGAAGGGGAAGTTGATTAATTTG  >  NC_000913/4440781‑4441268

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