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

CGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACAAA  >  NC_000913/3220059‑3220316                                                                                                                                    |                                                                                                                               cGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGt                                                                                                                         <  1:485768/139‑1 (MQ=255)        ggCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTc                                                                                                                  >  1:339374/1‑139 (MQ=255)                                                   ggCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCgg                                                                       >  1:531523/1‑139 (MQ=255)                                                               gtgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAgg                                                                                              <  2:801399/104‑1 (MQ=255)                                                               gtgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAgg                                                                                              >  1:801399/1‑104 (MQ=255)                                                                tgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgc                                                          >  2:702219/1‑139 (MQ=255)                                                                tgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgc                                                          <  2:339374/139‑1 (MQ=255)                                                                tgCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgc                                                          <  1:702219/139‑1 (MQ=255)                                                                  ccGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCAccgccg                                                        <  1:351250/139‑1 (MQ=255)                                                                     tttGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                   <  1:357351/77‑1 (MQ=255)                                                                     tttGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCt                                                                                                                   >  2:357351/1‑77 (MQ=255)                                                                       tGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacgt                                                   <  2:799369/139‑1 (MQ=255)                                                                          cAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGAAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAg                                                >  1:725833/1‑139 (MQ=255)                                                                                   aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAgg                                                                                              >  2:928541/1‑84 (MQ=255)                                                                                   aGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAgg                                                                                              <  1:928541/84‑1 (MQ=255)                                                                                        tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATTTGACTGCGGTGATCCTCGAACCGATTc                                                                                                        <  2:559509/69‑1 (MQ=255)                                                                                        tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTc                                                                                                        >  1:559509/1‑69 (MQ=255)                                                                                        tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                  <  1:412793/139‑1 (MQ=255)                                                                                        tGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                  <  1:427968/139‑1 (MQ=255)                                                                                           gcACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCgg                               <  1:81067/139‑1 (MQ=255)                                                                                             aCGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                  >  1:422184/1‑134 (MQ=255)                                                                                             aCGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                  <  2:422184/134‑1 (MQ=255)                                                                                               ggCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                         >  1:731505/1‑61 (MQ=255)                                                                                               ggCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGAtt                                                                                                         <  2:731505/61‑1 (MQ=255)                                                                                                   ctTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                  >  1:876443/1‑128 (MQ=255)                                                                                                   ctTAACGAGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGtt                                  <  2:876443/128‑1 (MQ=255)                                                                                                          aGTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCACCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCt                >  1:552228/1‑139 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAggg                                                                                                    >  2:546534/1‑54 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAggg                                                                                                    <  1:546534/54‑1 (MQ=255)                                                                                                           gTGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTc               >  2:743459/1‑139 (MQ=255)                                                                                                            tGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacgt                                                   <  2:232204/102‑1 (MQ=255)                                                                                                            tGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCcgtacgt                                                   >  1:232204/1‑102 (MQ=255)                                                                                                            tGCAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCg              >  1:162429/1‑139 (MQ=255)                                                                                                              cAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGAt            <  2:351914/139‑1 (MQ=255)                                                                                                              cAAAAAAACCGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGAt            <  1:743459/139‑1 (MQ=255)                                                                                                                      ccGGTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACa    >  1:454915/1‑139 (MQ=255)                                                                                                                        ggTGATGATGTGACTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACaaa  <  2:725833/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                               CGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCCGTTACTGCCGGGCTTCCGTCATGTGCCGTTTGGCAATATCGAAGCCATGCGCACGGCTCTTAACGAGTGCAAAAAAACCGGTGATGATGTGGCTGCGGTGATCCTCGAACCGATTCAGGGTGAAGGTGGCGTAATTCTGCCGCCGCCGGGCTATCTCACCGCCGTACGTAAGCTATGCGATGAGTTCGGCGCACTGATGATCCTCGATGAAGTACAAA  >  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: ‑