Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,190	G→A	A233T (GCT→ACT)	patA →	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,190	0	G	A	100.0%	183.1 / NA	57	A233T (GCT→ACT)	patA	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/32);  total (25/32)

CGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACAAA  >  NC_000913/3220059‑3220316                                                                                                                                    |                                                                                                                                 cGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                                           <  1:797710/139‑1 (MQ=255) cGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                                           <  2:1258695/139‑1 (MQ=255)        ggCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTc                                                                                                                    >  1:74795/1‑139 (MQ=255)        ggCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTc                                                                                                                    >  2:1077771/1‑139 (MQ=255)         gCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCg                                                                                                                   <  1:1087644/139‑1 (MQ=255)          cGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGa                                                                                                                  >  1:1008303/1‑139 (MQ=255)                 gACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                           <  1:810853/139‑1 (MQ=255)                 gACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                           <  1:165322/139‑1 (MQ=255)                 gACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCAAACCGAtt                                                                                                           <  1:1077771/139‑1 (MQ=255)                     ttCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAgg                                                                                                       <  1:544914/139‑1 (MQ=255)                               cGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTgg                                                                                             <  2:74795/139‑1 (MQ=255)                                gTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                     >  1:731326/1‑114 (MQ=255)                                gTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                     <  2:731326/114‑1 (MQ=255)                                 tttATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCg                                                                                           >  2:1143190/1‑139 (MQ=255)                                 tttATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCg                                                                                           >  2:780759/1‑139 (MQ=255)                                                ccGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                                            <  1:1143190/139‑1 (MQ=255)                                                   ggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCgg                                                                         <  2:1008303/139‑1 (MQ=255)                                                    gCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCggg                                                                        >  2:168483/1‑139 (MQ=255)                                                     cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAACCCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGc                                                                       <  1:1233357/139‑1 (MQ=255)                                                     cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGc                                                                       <  2:495425/139‑1 (MQ=255)                                                     cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGc                                                                       <  1:780759/139‑1 (MQ=255)                                                                tgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccgc                                                          <  1:168483/139‑1 (MQ=255)                                                                tgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccgc                                                          <  2:708117/139‑1 (MQ=255)                                                                  ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAgg                                                                                                >  2:643317/1‑101 (MQ=255)                                                                  ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAgg                                                                                                <  1:643317/101‑1 (MQ=255)                                                                       tGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCTGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccg                                                           <  2:846427/131‑1 (MQ=255)                                                                       tGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGc                                                                       <  1:486360/121‑1 (MQ=255)                                                                       tGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGc                                                                       >  2:486360/1‑121 (MQ=255)                                                                       tGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccg                                                           >  1:846427/1‑131 (MQ=255)                                                                                cGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCTCACCGCCGTACGTAAGCTAt                                            >  1:963046/1‑139 (MQ=255)                                                                                  aaGCAATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTTACCGCCGTACGTAAGATATGCGa                                        >  1:302942/1‑139 (MQ=255)                                                                                       aTGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGt                                   >  2:312110/1‑139 (MQ=255)                                                                                        tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                     >  1:372967/1‑58 (MQ=255)                                                                                        tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                     <  2:372967/58‑1 (MQ=255)                                                                                         gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGa                                                                                                   <  1:594563/75‑1 (MQ=255)                                                                                         gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGa                                                                                                   >  2:594563/1‑75 (MQ=255)                                                                                         gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTc                                 >  1:118191/1‑139 (MQ=255)                                                                                           gcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCgg                               <  1:312110/139‑1 (MQ=255)                                                                                           gcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCgg                               <  1:127626/139‑1 (MQ=255)                                                                                                gCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCAc                          <  2:1121130/139‑1 (MQ=255)                                                                                                 ctctTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTg                                                                                                    <  1:35283/66‑1 (MQ=255)                                                                                                 ctctTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTg                                                                                                    >  2:35283/1‑66 (MQ=255)                                                                                                   ctTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACtga                       >  1:606142/1‑139 (MQ=255)                                                                                                     tAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACtga                       <  2:963046/139‑1 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccgccg                                                        >  2:1232838/1‑98 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccgccg                                                        <  1:1232838/98‑1 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTc               >  2:1332851/1‑139 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTc               >  1:919095/1‑139 (MQ=255)                                                                                                             ggacaagggtcGGTGATGATGTGACTGCGGTGATGCTCGACCCGATTCTAGGCGGAGGTGGGGTAATCCTGCCGCGGCTGGGCTA‑‑CCTCACCGCCGTACGTAAGCTATGTGATTAGTTCGGCGCACTGTTGATCCTCGa             <  2:302942/129‑1 (MQ=255)                                                                                                                      ccGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTAAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACa    >  1:689841/1‑139 (MQ=255)                                                                                                                        ggTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACaaa  <  2:606142/139‑1 (MQ=255)                                                                                                                        ggTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACaaa  <  1:1332851/139‑1 (MQ=255)                                                                                                                         gTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACtgatg                     >  2:1040377/1‑119 (MQ=255)                                                                                                                         gTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACtgatg                     <  1:1040377/119‑1 (MQ=255)                                                                                                                              gatgTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccgc                                                          >  2:583729/1‑77 (MQ=255)                                                                                                                              gatgTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCAccgc                                                          <  1:583729/77‑1 (MQ=255)                                                                                                                               atgTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCgg                               >  2:1275852/1‑103 (MQ=255)                                                                                                                               atgTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCgg                               <  1:1275852/103‑1 (MQ=255)                                                                                                                                    |                                                                                                                                 CGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTA‑‑TCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACAAA  >  NC_000913/3220059‑3220316

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