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%	78.9 / NA	24	K23R (AAG→AGG)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/11);  total (13/11)

CGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTC  >  NC_000913/3219483‑3219631                                                                               |                                                                       cGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGc                                                                      <  2:538586/81‑1 (MQ=255)  ggAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCg                                                                     <  2:538587/81‑1 (MQ=255)     gCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAAc                                                                  >  2:538617/1‑81 (MQ=255)       atatTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGc                                                                <  2:538588/81‑1 (MQ=255)        tatTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCt                                                               >  2:538589/1‑81 (MQ=255)          ttttGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTg                                                              >  2:538611/1‑80 (MQ=255)              gAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCa                                                         >  2:538623/1‑81 (MQ=255)                       aCCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGAt                                                >  1:538620/1‑81 (MQ=255)                        ccTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATg                                               <  2:538590/81‑1 (MQ=255)                           tCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGa                                              >  2:538625/1‑79 (MQ=255)                             gAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAg                                           <  2:538591/80‑1 (MQ=255)                             gAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAg                                           >  2:538649/1‑80 (MQ=255)                             gAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAg                                           >  1:538617/1‑80 (MQ=255)                                     cGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAAcc                                  >  1:538621/1‑81 (MQ=255)                                      ggCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCg                                 <  2:538592/81‑1 (MQ=255)                                        cTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCgag                               <  2:538593/81‑1 (MQ=255)                                           tAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAgg                            >  2:538618/1‑81 (MQ=255)                                                  cAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTg                      <  1:538595/80‑1 (MQ=255)                                                  cAGCGCCCACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTg                      <  2:538595/80‑1 (MQ=255)                                                         cACGCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTc              <  1:538597/81‑1 (MQ=255)                                                            gCCCTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCa             >  1:538635/1‑79 (MQ=255)                                                              ccTGAATCTCATTGAGAGGCGAACGCTGGGTCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCaaa           >  2:538639/1‑79 (MQ=255)                                                              ccTGAATCTCATTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCaaa           >  2:538630/1‑79 (MQ=255)                                                                        aTTGAGAGGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTc  <  1:538598/78‑1 (MQ=255)                                                                               |                                                                       CGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTC  >  NC_000913/3219483‑3219631

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