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%	173.9 / NA	54	A233T (GCT→ACT)	patA	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/26);  total (28/26)

GGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGA  >  NC_000913/3220066‑3220305                                                                                                                             |                                                                                                                    ggCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTc                                                                                                       >  2:625805/1‑139 (MQ=255) ggCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTc                                                                                                       >  1:208190/1‑139 (MQ=255)                         gTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGc                                                                               <  2:545712/139‑1 (MQ=255)                                  gTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                              >  2:1410085/1‑99 (MQ=255)                                  gTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                              <  1:1410085/99‑1 (MQ=255)                                     aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                              <  2:172769/96‑1 (MQ=255)                                     aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                              >  1:172769/1‑96 (MQ=255)                                     aCTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATCATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccg                                                                   >  1:662081/1‑139 (MQ=255)                                       tgtcGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                 >  1:309052/4‑139 (MQ=255)                                         tcGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                               <  2:309052/138‑1 (MQ=255)                                         ccGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                               <  2:1388426/139‑1 (MQ=255)                                         ccGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                               <  2:662081/139‑1 (MQ=255)                                           gggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAAc                                                                                                   >  1:39631/1‑101 (MQ=255)                                           gggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAAc                                                                                                   <  2:39631/101‑1 (MQ=255)                                             gCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCggg                                                           <  2:1539345/139‑1 (MQ=255)                                                          gCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgcc                                              <  2:863243/139‑1 (MQ=255)                                                               gagacaggatCGAAGCCATGCGCACGTCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTCACTGCGGGGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacg                                         <  1:1062456/131‑1 (MQ=255)                                                                tGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacgt                                        >  2:1436113/1‑139 (MQ=255)                                                                        tCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATg                                >  2:580535/1‑139 (MQ=255)                                                                           aaGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                              >  1:1504621/1‑74 (MQ=255)                                                                           aaGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                              <  2:1504621/74‑1 (MQ=255)                                                                            aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGAt                            >  2:989112/1‑139 (MQ=255)                                                                                 tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGGTGATGTGACTACGGTGATCCTCGAACAGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGCTgtctc                       >  2:1062456/1‑135 (MQ=255)                                                                                 tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                               <  1:470913/99‑1 (MQ=255)                                                                                 tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                               >  2:470913/1‑99 (MQ=255)                                                                                 tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAg                                     <  1:435923/125‑1 (MQ=255)                                                                                 tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAg                                     >  2:435923/1‑125 (MQ=255)                                                                                 tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                       >  2:85871/1‑139 (MQ=255)                                                                                  gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAAcc                                                                                                  >  2:396307/1‑63 (MQ=255)                                                                                  gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAAcc                                                                                                  <  1:396307/63‑1 (MQ=255)                                                                                  gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgc                                               <  2:1343881/114‑1 (MQ=255)                                                                                  gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgc                                               >  1:1343881/1‑114 (MQ=255)                                                                                  gcgcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTc                      >  1:1510623/1‑139 (MQ=255)                                                                                   cgcACTGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCTCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCg                     <  1:625805/139‑1 (MQ=255)                                                                                      aCGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGgcg                  <  2:208190/139‑1 (MQ=255)                                                                                         gCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCAc               <  1:580535/139‑1 (MQ=255)                                                                                                aCGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                        <  2:998563/43‑1 (MQ=255)                                                                                                aCGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                        >  1:998563/1‑43 (MQ=255)                                                                                                    gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                              >  1:1033020/1‑49 (MQ=255)                                                                                                    gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                              <  2:1033020/49‑1 (MQ=255)                                                                                                    gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccg                                                                   <  1:1393894/76‑1 (MQ=255)                                                                                                    gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccg                                                                   >  2:1393894/1‑76 (MQ=255)                                                                                                    gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGgcg                  <  1:979188/125‑1 (MQ=255)                                                                                                    gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGgcg                  >  2:979188/1‑125 (MQ=255)                                                                                                     tGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCg   >  1:1318850/1‑139 (MQ=255)                                                                                                      gCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGa  >  2:74280/1‑139 (MQ=255)                                                                                                       cAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAt                                                                                               >  2:1087835/1‑45 (MQ=255)                                                                                                       cAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAt                                                                                               <  1:1087835/45‑1 (MQ=255)                                                                                                             aaCCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                       >  1:49540/1‑111 (MQ=255)                                                                                                             aaCCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                       <  2:49540/111‑1 (MQ=255)                                                                                                                  gTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                 >  1:1508087/1‑64 (MQ=255)                                                                                                                  gTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                 <  2:1508087/64‑1 (MQ=255)                                                                                                                      tgatgTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacg                                         >  1:1537754/1‑84 (MQ=255)                                                                                                                      tgatgTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacg                                         <  2:1537754/84‑1 (MQ=255)                                                                                                                             |                                                                                                                    GGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGA  >  NC_000913/3220066‑3220305

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: ‑