Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,219,963	C→T	A157V (GCG→GTG)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,219,963	0	C	T	100.0%	188.4 / NA	53	A157V (GCG→GTG)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/23);  total (30/23)

CGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCC  >  NC_000913/3219828‑3220098                                                                                                                                        |                                                                                                                                        cGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTg                                                                                                                                     >  1:108060/1‑140 (MQ=255)         ccGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGc                                                                                                                             <  2:687757/140‑1 (MQ=255)               cacaGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGa                                                                                                                           <  1:1123550/136‑1 (MQ=255)               cacaGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGa                                                                                                                           >  2:1123550/1‑136 (MQ=255)                acaGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCt                                                                                                                      >  2:1574186/1‑140 (MQ=255)                acaGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCt                                                                                                                      >  2:321604/1‑140 (MQ=255)                    ccAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTAcc                                                                                                                  >  2:1605485/1‑140 (MQ=255)                    ccAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTAcc                                                                                                                  <  1:1003443/140‑1 (MQ=255)                         aGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCa                                                                                                             >  2:1262654/1‑140 (MQ=255)                          gctgctCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCAc                                                                                                            <  1:1605485/140‑1 (MQ=255)                              ctCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACcgcg                                                                                                        >  1:209035/1‑140 (MQ=255)                              ctCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACcgcg                                                                                                        >  2:1006401/1‑140 (MQ=255)                              ctCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACcgcg                                                                                                        >  2:1217739/1‑140 (MQ=255)                              ctCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACcgcg                                                                                                        >  2:556499/1‑140 (MQ=255)                                 gATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCgg                                                                                                     <  1:1217739/140‑1 (MQ=255)                                  aTCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGc                                                                                                    <  1:321604/140‑1 (MQ=255)                                       ttACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGtt                                                                                               >  1:821602/1‑140 (MQ=255)                                       ttACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGtt                                                                                               <  1:590961/140‑1 (MQ=255)                                           gggCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACt                                                                                           >  1:540783/1‑140 (MQ=255)                                             gCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACttt                                                                                         >  1:872616/1‑140 (MQ=255)                                                        gAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGC‑GCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAg                                                                             <  1:1584390/140‑1 (MQ=255)                                                             cccTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAgcggc                                                                         <  2:209035/140‑1 (MQ=255)                                                             cccTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAgcggc                                                                         >  2:1422404/1‑140 (MQ=255)                                                                ttGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGc                                                                                   <  2:396195/127‑1 (MQ=255)                                                                ttGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGc                                                                                   >  1:396195/1‑127 (MQ=255)                                                                 tgctgcGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGt                                                                     >  1:728882/1‑140 (MQ=255)                                                                      cgcTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCAc                                                                <  1:556499/140‑1 (MQ=255)                                                                       gcTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACg                                                               >  2:522485/1‑140 (MQ=255)                                                                            cGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTaaa                                                          <  2:368110/140‑1 (MQ=255)                                                                                cGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCAc                                                      <  2:728882/140‑1 (MQ=255)                                                                                cGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCAc                                                      >  1:554143/1‑140 (MQ=255)                                                                                      aCTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCGCCAGCGGCGCGTTCCACGGTAAATCACTTGgcg                                                <  2:1597282/140‑1 (MQ=255)                                                                                      aCTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGgcg                                                <  2:1173483/140‑1 (MQ=255)                                                                                        tGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGgcgcg                                              <  1:572580/140‑1 (MQ=255)                                                                                                  cttcttCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCa                                    >  2:447343/1‑140 (MQ=255)                                                                                                  cttcttCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCa                                    >  2:1153525/1‑140 (MQ=255)                                                                                                    tcttctGTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACg                                  >  2:510258/1‑140 (MQ=255)                                                                                                          gTAATAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGCCACCCCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCCCGCTGTCGGCCACGGCGaaa                            >  2:601058/1‑140 (MQ=255)                                                                                                             aTAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCg                         >  2:572017/1‑140 (MQ=255)                                                                                                              tAGCGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGa                        >  1:702459/1‑140 (MQ=255)                                                                                                                 cGGCACCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCt                     <  2:872616/140‑1 (MQ=255)                                                                                                                     aCCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAGTCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCg                 >  1:114933/1‑140 (MQ=255)                                                                                                                     aCCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCg                 >  1:295997/1‑140 (MQ=255)                                                                                                                     aCCGAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCg                 >  2:841174/1‑140 (MQ=255)                                                                                                                        gAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCaa              <  2:429363/140‑1 (MQ=255)                                                                                                                        gAGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCaa              <  1:1262654/140‑1 (MQ=255)                                                                                                                         aGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGAc                       >  1:1472033/1‑130 (MQ=255)                                                                                                                         aGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGAc                       <  2:1472033/130‑1 (MQ=255)                                                                                                                         aGTCCGTTGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCaaa             <  2:540783/140‑1 (MQ=255)                                                                                                                               ttGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGAc                       <  2:584485/124‑1 (MQ=255)                                                                                                                               ttGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGAc                       >  1:584485/1‑124 (MQ=255)                                                                                                                                tGAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTa      >  1:54249/1‑140 (MQ=255)                                                                                                                                 gAAGCAGTGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTAt     <  2:108060/140‑1 (MQ=255)                                                                                                                                   agcagtGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGc   >  2:1225519/1‑140 (MQ=255)                                                                                                                                    gcagtGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGcc  >  2:1570200/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        CGAAACAACCGCTGCACAGCCAGGAGCTGCTCGATCCGTTACGGGCGATGTTGGCGAAAACCCTTGCTGCGCTAACGCCCGGTAAACTGAAATACAGCTTCTTCTGTAATAGCGGCACCGAGTCCGTTGAAGCAGCGCTGAAGCTGGCGAAAGCTTACCAGTCACCGCGCGGCAAGTTTACTTTTATTGCCACCAGCGGCGCGTTCCACGGTAAATCACTTGGCGCGCTGTCGGCCACGGCGAAATCGACCTTCCGCAAACCGTTTATGCC  >  NC_000913/3219828‑3220098

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