Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,829	C→A	R446S (CGT→AGT)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,829	0	C	A	100.0%	166.1 / NA	47	R446S (CGT→AGT)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (30/17);  total (30/17)

GCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCA  >  NC_000913/3220695‑3220964                                                                                                                                       |                                                                                                                                         gCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAg                                                                                                                                     >  1:1023134/1‑140 (MQ=255)     taACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGgcgc                                                                                                                                 >  1:253730/1‑140 (MQ=255)                 cGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGc                                                                                                                     >  2:689938/1‑140 (MQ=255)                           cGCCGGCGCGTTCGGGGGGCGGAACCCCTAAAAAACGCAAAAACGCCCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAgtgt                                                                                                           <  2:322547/140‑1 (MQ=255)                                gcgcGTACTGGTGGCCGGAACGCTCAATAACGCAAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCgaag                                                                                                      <  2:1014126/140‑1 (MQ=255)                                gcgcGTACGGGGGGCCGAAACGCTCAATAACGCAAAAAGAACCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCgaag                                                                                                      <  2:1023134/140‑1 (MQ=255)                                  gcgTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCgaagaa                                                                                                    >  2:203942/1‑140 (MQ=255)                                     tattgggggCGGGACCGCTAAAAAACGCAAAAAGGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCg                                                                                                 <  1:396576/137‑1 (MQ=255)                                     tACTGGTGGCCGGAACGCTCAATAACGCAAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCg                                                                                                 <  2:55579/140‑1 (MQ=255)                                      aCTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAG‑AGCGTa                                                                                               >  1:345551/1‑140 (MQ=255)                                             gCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATcc                                                                                         >  1:972764/1‑139 (MQ=255)                                                  aaCGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCg                                                                                    >  1:902031/1‑140 (MQ=255)                                                       tCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATgg                                                                               >  2:1009730/1‑140 (MQ=255)                                                                 cAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATacgacg                                                                     >  2:780260/1‑140 (MQ=255)                                                                   aaaCGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTc                                                                   >  1:757243/1‑140 (MQ=255)                                                                        aTCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCAt                                                             >  1:560842/1‑140 (MQ=255)                                                                            gCATTGACCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGa                                                         <  2:412646/140‑1 (MQ=255)                                                                              aTTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATc                                                       >  1:1063581/1‑140 (MQ=255)                                                                                 gAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCttt                                                    >  2:214121/1‑140 (MQ=255)                                                                                  aaCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCtttt                                                   >  2:1202612/1‑140 (MQ=255)                                                                                  aaCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCCGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCtttt                                                   >  1:30337/1‑140 (MQ=255)                                                                                     cGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTtat                                                >  2:765949/1‑140 (MQ=255)                                                                                        cactgacactgaCCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATacg                                             <  2:902031/140‑1 (MQ=255)                                                                                            gacactgacCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGAcgc                                         >  2:657810/1‑140 (MQ=255)                                                                                                   acCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGcc                                  >  2:712871/1‑140 (MQ=255)                                                                                                   acCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGcc                                  <  1:1202612/140‑1 (MQ=255)                                                                                                     cATTGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCg                                <  2:1063581/140‑1 (MQ=255)                                                                                                       ttGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCg                                <  2:496899/138‑1 (MQ=255)                                                                                                       ttGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCg                                >  1:496899/1‑138 (MQ=255)                                                                                                       ttGAACAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCa                              >  1:1064171/1‑140 (MQ=255)                                                                                                          aaCAGTGTGAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAG‑AGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGca                          <  2:345551/140‑1 (MQ=255)                                                                                                               tgtgAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGcatcat                      >  1:398340/1‑140 (MQ=255)                                                                                                                gtgAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCatcatc                     <  1:765949/140‑1 (MQ=255)                                                                                                                 tgAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCatcatc                     >  1:217300/1‑140 (MQ=255)                                                                                                                 tgAACTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATcc                    >  1:513183/1‑140 (MQ=255)                                                                                                                    aCTGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                 <  2:1064171/140‑1 (MQ=255)                                                                                                                      tGGTGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTcc               <  1:1009730/140‑1 (MQ=255)                                                                                                                         tGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                 >  2:633091/1‑135 (MQ=255)                                                                                                                         tGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                 <  1:633091/135‑1 (MQ=255)                                                                                                                         tGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGc            >  2:507045/1‑140 (MQ=255)                                                                                                                         tGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGc            >  2:760148/1‑140 (MQ=255)                                                                                                                         tGATCAAAGCGGCGAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGc            >  2:1043987/1‑140 (MQ=255)                                                                                                                                  cggcgAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGc   >  1:1108750/1‑140 (MQ=255)                                                                                                                                  cggcgAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGc   <  1:657810/140‑1 (MQ=255)                                                                                                                                   ggcgAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGc   <  2:217300/140‑1 (MQ=255)                                                                                                                                   ggcgAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGcc  >  2:246022/1‑139 (MQ=255)                                                                                                                                   ggcgAGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCCTTTTATACGACGCGTACGCCCGCAGGCAACATCCGCTCCGGCGTTAACAGCa  >  2:261366/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                         GCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCG‑CCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCA  >  NC_000913/3220695‑3220964

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