Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,219,561	A→G	K23R (AAG→AGG)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,219,561	0	A	G	100.0%	134.7 / NA	40	K23R (AAG→AGG)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/25);  total (15/25)

CGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTC  >  NC_000913/3219483‑3219631                                                                               |                                                                       cGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGc                                                                      <  2:599866/81‑1 (MQ=255) cGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGc                                                                      <  2:599865/81‑1 (MQ=255)  ggAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCg                                                                     <  2:599869/81‑1 (MQ=255)  ggAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCg                                                                     <  2:599868/81‑1 (MQ=255)  ggAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCg                                                                     <  2:599867/81‑1 (MQ=255)    aGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGaa                                                                   >  2:599931/1‑81 (MQ=255)     gCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAAc                                                                  >  1:599914/1‑81 (MQ=255)       atatTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGc                                                                >  1:599928/1‑81 (MQ=255)         atTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTg                                                              >  1:599941/1‑81 (MQ=255)          ttttGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTg                                                              <  1:599870/80‑1 (MQ=255)            ttGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGa                                                            <  2:599871/80‑1 (MQ=255)              gAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATc                                                          <  1:599872/80‑1 (MQ=255)               aaCAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCAt                                                        <  1:599873/81‑1 (MQ=255)                 cAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCAt                                                        >  1:599926/1‑79 (MQ=255)                  aGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATgag                                                     <  1:599874/81‑1 (MQ=255)                    gTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATgagga                                                   <  1:599875/81‑1 (MQ=255)                     ttACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATgaggag                                                  <  1:599877/81‑1 (MQ=255)                     ttACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATgaggag                                                  <  2:599876/81‑1 (MQ=255)                        ccTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGAt                                                >  1:599936/1‑80 (MQ=255)                             gAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAg                                           >  2:599922/1‑80 (MQ=255)                             gAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAg                                           >  2:599921/1‑80 (MQ=255)                                   aTCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTaa                                    >  2:599932/1‑81 (MQ=255)                                      ggCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAAcc                                  <  1:599878/80‑1 (MQ=255)                                      ggCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAAcc                                  >  2:599920/1‑80 (MQ=255)                                       gCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAAcc                                  >  2:599938/1‑79 (MQ=255)                                        cTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCgag                               <  2:599879/81‑1 (MQ=255)                                             gCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTg                          <  2:599880/81‑1 (MQ=255)                                              cGTGCAGCGCCCACTCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGa                         <  2:599881/81‑1 (MQ=255)                                              agagCAGCGCCCAAGCCATGAATCTCATTGAGAGGCAAACGCTGGATCATGAGGAGATGAAAGCACTTACCCGAGAGGTGa                         <  2:599882/78‑1 (MQ=255)                                               gTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGAt                        >  1:599930/1‑81 (MQ=255)                                                   aGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGaa                    <  2:599883/81‑1 (MQ=255)                                                   aGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGaa                    >  2:599937/1‑81 (MQ=255)                                                    gcgcCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAAt                   <  2:599884/81‑1 (MQ=255)                                                            gCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCaaa           <  2:599886/81‑1 (MQ=255)                                                            gCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCaaa           <  1:599885/81‑1 (MQ=255)                                                            gCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCa             >  2:599943/1‑79 (MQ=255)                                                               cTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAagag        >  1:599948/1‑81 (MQ=255)                                                                  aaTCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGc       <  1:599887/79‑1 (MQ=255)                                                                   aTCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATg    <  2:599888/81‑1 (MQ=255)                                                                        aTTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTc  <  1:599889/78‑1 (MQ=255)                                                                               |                                                                       CGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTC  >  NC_000913/3219483‑3219631

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