Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,571	G→A	E360K (GAG→AAG)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,571	0	G	A	100.0%	39.0 / NA	17	E360K (GAG→AAG)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/9);  total (8/9)

GCGAGCATGAGAACGTACAGCCGGATATCCTCTGCCTTGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTG  >  NC_000913/3220350‑3220816                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      gCGAGCATGAGAACGTACAGCCGGATATCCTCTGCCTTGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACg                                                                                                                                                                                                                            <  1:152030/249‑1 (MQ=255)          agaACGTACAGCCGGATATCCTCTGCCTTGCCAAAGCGCTCGGCGGCGGAGTGATGCTGATTGGCGCGCCCATCGCCCCTGAAGCGGTGCTTTTAGTTCCGTTCGTCAACCCCTTACTGCCTACCACCACCTTTGGCGGCAACCCGCTGCCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCCGAACTTACCGGCACAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCa                                                                                                                                                                                                                  >  1:237641/1‑250 (MQ=255)                       ggATATCCTCTGCCTTGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCgg                                                                                                                                                                                                        <  1:248456/247‑1 (MQ=255)                                                   ngcggcggcgTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTTGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCCTCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCCGTCTGGGACAGGAAGcgcgggg                                                                                                                                                                        >  1:82630/2‑248 (MQ=255)                                                                    gATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGGTGGCGACCATCAATGTGTTGCTGGAGCAGAACGTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGGTTCCGTCAACTGGCGCGGGAATATCCCGAGCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATgg                                                                                                                                                       <  1:253505/251‑1 (MQ=255)                                                                                cATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGAGACCATCAATGTGTTGCTGGATCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGAc                                                                                                                                                                                                                             >  1:74348/1‑169 (MQ=255)                                                                                 aTCGCCACTGAAGAGGTGTTTTCAGTGCTGTTCGACAACCCATTCCTGCATATCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCGGGCGACCATCAATGTGTTGCTTTAGCAGAACTTACCGGCTCAGGCTAAGCTAAAAGGCGGGATGTTGCTGGACGTTTTCCGTCAACTGGAGCGGTAATATCCCGATCTGGTACAGGAAGCGCGTGTTAAAGGGATGTTGATGGCGATTGAGTttgt                                                                                                                                          <  1:406311/251‑1 (MQ=255)                                                                                                                          aTTCCTGCATACCACCCCCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTcn                                                                                                 <  1:113335/251‑2 (MQ=255)                                                                                                                                                ggcggcAACGCGCTGGCCTGTGCGGCGGCGCTGACGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCCGGCTAAGCAAAAAGGCGATATGTTGATGGACGGTTTCCGTCAACTGGAGCGGGAATATTCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGAACTGGTGGCCgg                                                                           >  1:369090/1‑251 (MQ=255)                                                                                                                                                ggcggcAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGTCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACtgg                                                                                  <  1:219352/244‑1 (MQ=255)                                                                                                                                                               gCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTcc                                                                                                 <  1:144666/214‑1 (MQ=255)                                                                                                                                                               gCCTGGGCGTCGGCGCTGGCGACCATCAAGGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTTGACGGTTTCCGTCAACTGGCTCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATAGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACGTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGc                                                            <  1:297289/251‑1 (MQ=255)                                                                                                                                                                                            tgtgTTGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACGGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGGTGATGGCGATTGAGTGTGTTGATAACg                                                                                                                                  <  1:143083/152‑1 (MQ=255)                                                                                                                                                                                                 tGCTGGAGCAGAACTTACCGGCTCAGGCTAAGCAAAACGGCGATATGTTGATGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCcactgacactg                          >  1:172949/1‑251 (MQ=255)                                                                                                                                                                                                                        cAGGCTAAGCAAAAAGGCGATATGTTGCTGGACGGCTTACGTCAACTGGCGCGGGAATATCCGCATATAGCACAGGACGCGCGTGGTGAAGGGATGTTGATGGCGATTGAGTTTGTTGATAAAGACATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGACGGAACGCTCAATAACGCCAAAACGATCCCCCTTGAACCGCCACTGCCACTGACCATTGCATCGTGTGAACTGGt   >  1:386193/1‑251 (MQ=255)                                                                                                                                                                                                                          ggCTAAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTg  >  1:77481/1‑250 (MQ=255)                                                                                                                                                                                                                          ggCTAAGCAAAAAGGCGATATGTGGCTGGACGGTTTCCGAAAACTGGCGCGGGAATATCCGGATCTGGAACAGGAAGAGCGTGGTAAAGGGATGCTGATGACGACTGAGTTCGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGAAATGGCGGCCGGAACGCACAATAACGACAAAACGATCCGCATTGAACCGACACTGACAATGACCAGTGACAAGGGTGAACTGGTg  >  1:77490/1‑250 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      GCGAGCATGAGAACGTACAGCCGGATATCCTCTGCCTTGCCAAAGCGCTCGGCGGCGGCGTGATGCCGATTGGCGCGACCATCGCCACTGAAGAGGTGTTTTCAGTTCTGTTCGACAACCCATTCCTGCATACCACCACCTTTGGCGGCAACCCGCTGGCCTGTGCGGCGGCGCTGGCGACCATCAATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTG  >  NC_000913/3220350‑3220816

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