Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,219,596	A→G	N35D (AAC→GAC)	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,219,596	0	A	G	100.0%	174.0 / NA	49	N35D (AAC→GAC)	patA	putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/30);  total (19/30)

CTTAAATCGGTATCATGTGATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGG  >  NC_000913/3219451‑3219735                                                                                                                                                  |                                                                                                                                            cTTAAATCGGTATCATGTGATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCg                                                                                                                                         >  2:289955/1‑150 (MQ=255)   tAAATCGGTATCATGTGATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCga                                                                                                                                        <  1:588767/149‑1 (MQ=255)    aaaTCGGTATCATGTGATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGagag                                                                                                                                     <  2:183611/151‑1 (MQ=255)            aTCATGTGATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTg                                                                                                                              >  1:373364/1‑150 (MQ=255)                 gtgATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATAc                                                                                                                         <  1:210144/150‑1 (MQ=255)                 gtgATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATAc                                                                                                                         <  1:412031/150‑1 (MQ=255)                 gtgATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACt                                                                                                                        <  2:37054/151‑1 (MQ=255)                      aCGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCaa                                                                                                                    <  2:373388/150‑1 (MQ=255)                      aCGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCaa                                                                                                                    >  1:766432/1‑150 (MQ=255)                         cgAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAagag                                                                                                                <  2:739251/151‑1 (MQ=255)                         cgAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAagag                                                                                                                <  2:625110/151‑1 (MQ=255)                         cgAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAagag                                                                                                                <  2:166713/151‑1 (MQ=255)                         cgAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAag                                                                                                                  >  1:401549/1‑149 (MQ=255)                           aGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGc                                                                                                               >  2:89480/1‑150 (MQ=255)                                ccGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTc                                                                                                          >  1:116979/1‑150 (MQ=255)                                ccGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGACCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTc                                                                                                          <  1:282301/150‑1 (MQ=255)                                  ggAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAAt                                                                                                       >  1:525010/1‑151 (MQ=255)                                          ttttGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGtt                                                                                               >  2:254161/1‑151 (MQ=255)                                                  aGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGt                                                                                       >  2:255614/1‑151 (MQ=255)                                                       aCCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCg                                                                                   <  2:715191/150‑1 (MQ=255)                                                          ttCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCaa                                                                                >  2:691583/1‑150 (MQ=255)                                                              aGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCt                                                                            >  1:591346/1‑150 (MQ=255)                                                                  cATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTAc                                                                       <  2:300858/151‑1 (MQ=255)                                                                     cGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTAccg                                                                     <  2:636422/150‑1 (MQ=255)                                                                                 gCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGAtt                                                         <  2:853962/150‑1 (MQ=255)                                                                                   aGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGTTGGATCATGAGGAGATGAAAGCGCTTGACCGAGAGGTGATTGAATCCTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCGGTTACCGCCGGCGGGGATTAc                                                       <  1:685366/150‑1 (MQ=255)                                                                                   aGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTAc                                                       <  2:125487/150‑1 (MQ=255)                                                                                   aGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACg                                                      >  1:346583/1‑151 (MQ=255)                                                                                    gcgcCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACgg                                                     <  2:808569/151‑1 (MQ=255)                                                                                        ccACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGcc                                                 <  2:794341/151‑1 (MQ=255)                                                                                         cACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCg                                                <  2:329529/151‑1 (MQ=255)                                                                                            gCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGt                                               >  2:722682/1‑149 (MQ=255)                                                                                               cTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTaga                                            >  1:390087/1‑149 (MQ=255)                                                                                                 gAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTgg                                        <  1:89457/151‑1 (MQ=255)                                                                                                     ctcATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCaa                                     <  2:780173/150‑1 (MQ=255)                                                                                                        aTTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCgg                                 <  1:686921/151‑1 (MQ=255)                                                                                                          tGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCggg                                <  2:138906/150‑1 (MQ=255)                                                                                                            agaAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCggg                                >  1:509218/1‑148 (MQ=255)                                                                                                              aaGCGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGt                            <  1:92031/150‑1 (MQ=255)                                                                                                                 cGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTaa                        <  1:714173/151‑1 (MQ=255)                                                                                                                 cGAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTa                         <  2:33196/150‑1 (MQ=255)                                                                                                                  gAACGCTGGATCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTaa                        >  1:545286/1‑150 (MQ=255)                                                                                                                           aTCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGt              >  2:776002/1‑151 (MQ=255)                                                                                                                            tCATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTc             >  2:243357/1‑151 (MQ=255)                                                                                                                             cATGAGGAGATGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCg            >  2:527633/1‑151 (MQ=255)                                                                                                                                 agaggaTGAAAGCACTTGACTGGGAGGTGACTGAATGCTTCAAAGAACATGTGAGCCCGGGGTTTTTAGAGTACCGCAAGTCTTTTACCGCCGCCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAGATACCCTTGTCGaca         <  1:775590/146‑1 (MQ=255)                                                                                                                                     gaTGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCAGGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAg    <  1:602898/151‑1 (MQ=255)                                                                                                                                      aTGAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAgg   <  2:263284/151‑1 (MQ=255)                                                                                                                                        gAAAGCACTTGACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAggg  <  2:546770/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                            CTTAAATCGGTATCATGTGATACGCGAGCCTCCGGAGCATATTTTGAACAGGTTACCTTCGAGCGCATCGGCTTTAGCGTGCAGCGCCCACGCCCTGAATCTCATTGAGAAGCGAACGCTGGATCATGAGGAGATGAAAGCACTTAACCGAGAGGTGATTGAATACTTCAAAGAGCATGTCAATCCGGGGTTTTTAGAGTATCGCAAATCTGTTACCGCCGGCGGGGATTACGGAGCCGTAGAGTGGCAAGCGGGAAGTTTAAATACGCTTGTCGACACCCAGGG  >  NC_000913/3219451‑3219735

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