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

GTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACAAAC  >  NC_000913/3220057‑3220317                                                                                                                                      |                                                                                                                                gTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCg                                                                                                                            >  2:514843/1‑139 (MQ=255)         cGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                    <  1:483986/139‑1 (MQ=255)            cGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                    <  1:869685/136‑1 (MQ=255)            cGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                    >  2:869685/1‑136 (MQ=255)            cGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGa                                                                                                                 >  2:758970/1‑139 (MQ=255)             gAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGaa                                                                                                                <  1:737656/139‑1 (MQ=255)               aaTCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAAcc                                                                                                              <  2:848720/139‑1 (MQ=255)                   gACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                          <  1:76501/139‑1 (MQ=255)                   gACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                          >  1:973790/1‑139 (MQ=255)                          cGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTg                                                                                                   >  2:620135/1‑139 (MQ=255)                               aCCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGt                                                                                                    >  1:981539/1‑133 (MQ=255)                               aCCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGt                                                                                                    <  2:981539/133‑1 (MQ=255)                                           gTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTg                                                                                  >  2:988329/1‑139 (MQ=255)                                                  ccGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                                           <  2:973790/139‑1 (MQ=255)                                                  ccGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                                           <  1:758970/139‑1 (MQ=255)                                                       cTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGc                                                                      <  1:583515/139‑1 (MQ=255)                                                                  tgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgc                                                           <  1:195322/139‑1 (MQ=255)                                                                    ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                             >  1:713030/1‑119 (MQ=255)                                                                    ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                             <  2:713030/119‑1 (MQ=255)                                                                    ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCCCAccgc                                                           <  1:495716/137‑1 (MQ=255)                                                                    ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCCCAccgc                                                           >  2:495716/1‑137 (MQ=255)                                                                       tttGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAc                                                              >  1:358422/1‑131 (MQ=255)                                                                       tttGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGCGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAc                                                              <  2:358422/131‑1 (MQ=255)                                                                            cAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAg                                                 >  2:1094132/1‑139 (MQ=255)                                                                                    aaGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                                          >  1:544603/1‑58 (MQ=255)                                                                                    aaGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                                          <  2:544603/58‑1 (MQ=255)                                                                                    aaGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                          >  2:239514/1‑74 (MQ=255)                                                                                    aaGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                          <  1:239514/74‑1 (MQ=255)                                                                                     aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGc                                                <  1:367105/131‑1 (MQ=255)                                                                                     aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGc                                                >  2:367105/1‑131 (MQ=255)                                                                                     aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGa                                         <  1:398278/138‑1 (MQ=255)                                                                                     aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGa                                         >  2:398278/1‑138 (MQ=255)                                                                                          tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                   >  2:477080/1‑139 (MQ=255)                                                                                          tGCGCACGGCTCGTAACGAGTGCAAAAAAACCGGAGATGATGTCACTGCGGTGATCCTCGAACCGATTCAGGGAGAAGGTGGCGTAATTCTTCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                   >  2:495284/1‑139 (MQ=255)                                                                                           gcgcACGGCTCTTAACGAGTGCAAAAAAACAGGTGATGATGTGACCGCGGTGATCTTCGAACCGATTTAGGGTGATGGTGGCGTAATTATGCCGCCGCCTGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTc                                  >  2:583868/1‑139 (MQ=255)                                                                                             gcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCgg                                <  1:620135/139‑1 (MQ=255)                                                                                                cGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGgcgc                             <  2:474666/139‑1 (MQ=255)                                                                                                cGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGgcgc                             <  2:338934/139‑1 (MQ=255)                                                                                                 ggCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCa                            <  1:477080/139‑1 (MQ=255)                                                                                                        aaCGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                             >  2:215466/1‑83 (MQ=255)                                                                                                        aaCGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgcc                                                                             <  1:215466/83‑1 (MQ=255)                                                                                                             gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                                           >  1:1042208/1‑80 (MQ=255)                                                                                                             gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTc                >  2:40632/1‑139 (MQ=255)                                                                                                             gTGCAAAAAAACCGGTGATGATATGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTgccgccgc                                                                           <  2:1042208/80‑1 (MQ=255)                                                                                                                cAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGAt             <  2:512159/139‑1 (MQ=255)                                                                                                                          ggTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCC‑CGATGAAGTACAAAc  >  2:397648/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                GTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACAAAC  >  NC_000913/3220057‑3220317

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