Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	341,276	A→G	D51G (GAC→GGC)	yahJ →	putative deaminase with metallo‑dependent hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	341,276	0	A	G	100.0%	155.4 / NA	44	D51G (GAC→GGC)	yahJ	putative deaminase with metallo‑dependent hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/21);  total (23/21)

TAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAG  >  NC_000913/341139‑341412                                                                                                                                          |                                                                                                                                         tAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGGTAAGCGGACAACACCATGAAAATCACTGGcc                                                                                                                                        >  2:640656/1‑140 (MQ=255)  aGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGTGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGccc                                                                                                                                       <  2:434351/140‑1 (MQ=255)  aGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGccc                                                                                                                                       <  2:562157/140‑1 (MQ=255)   gccgccGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCg                                                                                                                                      <  2:579053/140‑1 (MQ=255)   gccgccGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCg                                                                                                                                      >  2:802323/1‑140 (MQ=255)      gccgTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCAt                                                                                                                                   >  1:243783/1‑140 (MQ=255)       ccgTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCAtt                                                                                                                                  <  2:222204/140‑1 (MQ=255)       ccgTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCAtt                                                                                                                                  <  1:802323/140‑1 (MQ=255)        cgTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTa                                                                                                                                 >  1:677725/1‑140 (MQ=255)           gAATTTCTGAGCCAGAGCGGTAAGTAGGACACCGCCGCCGCGCTGTTTGGTACCTATGTGCCGCACGCCCATGCGGCGGTAGCTGGCTCCCTAAACTGCGAAGCGAACAACACCAGGAAAATCACTGGCCCGCATTACTa                                                                                                                              <  2:416006/140‑1 (MQ=255)            aaTTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAGAATCACTGGCCCGCATTACTAt                                                                                                                             <  2:285423/140‑1 (MQ=255)                        agagCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGtg                                                                                                                 >  2:552661/1‑140 (MQ=255)                          agCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGtgct                                                                                                               >  2:748956/1‑140 (MQ=255)                               tAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTgg                                                                                                          >  2:724026/1‑140 (MQ=255)                               tAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTgg                                                                                                          >  1:335170/1‑140 (MQ=255)                                aaGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGa                                                                                                         <  1:640656/140‑1 (MQ=255)                                  gATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGaaa                                                                                                       <  2:506770/140‑1 (MQ=255)                                             ccgccgcGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACGACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGAc                                                                                            >  2:121957/1‑140 (MQ=255)                                             ccgccgcGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGAc                                                                                            <  1:552661/140‑1 (MQ=255)                                                        tttGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAAtgg                                                                                 <  1:556531/140‑1 (MQ=255)                                                        tttGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAAtgg                                                                                 <  1:399810/140‑1 (MQ=255)                                                          tGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAAtggcg                                                                               <  2:109923/140‑1 (MQ=255)                                                          tGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAAtggcg                                                                               <  2:243783/140‑1 (MQ=255)                                                             tACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAAtggcgtgg                                                                            <  1:748956/140‑1 (MQ=255)                                                                ctctGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACCCCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAACCCGGTTTTGCCTACAAAAATGGCGTGGCgg                                                                         <  2:335170/140‑1 (MQ=255)                                                                  ctGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTg                                                                       >  1:876349/1‑140 (MQ=255)                                                                  ctGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAACCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTg                                                                       >  1:764561/1‑140 (MQ=255)                                                                      gCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGa                                                                   >  2:823846/1‑140 (MQ=255)                                                                       ccGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGAc                                                                  >  1:596597/1‑140 (MQ=255)                                                                          cTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCg                                                               >  1:317922/1‑140 (MQ=255)                                                                              cccATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACg                                                           <  1:823846/140‑1 (MQ=255)                                                                               ccATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACgg                                                          >  2:865961/1‑140 (MQ=255)                                                                                     ggcggTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGcc                                                    <  1:757757/140‑1 (MQ=255)                                                                                      gcggTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAACCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCa                                                   <  2:764561/140‑1 (MQ=255)                                                                                                 cACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGGGCAGACCCCCCCGGCGCCCCAGACCggggggg                                        >  2:445169/1‑137 (MQ=255)                                                                                                             cgaagcgaaCAACACCATGAAAATCACTGGCCCGCATTACTATTTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCa                            >  2:551202/1‑140 (MQ=255)                                                                                                             cgaagcgaaCAACACCATGAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCa                            >  2:724215/1‑140 (MQ=255)                                                                                                                               gAAAATCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGcgcg          >  2:904794/1‑140 (MQ=255)                                                                                                                                 aaaTCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCgag        >  1:791697/1‑140 (MQ=255)                                                                                                                                   aTCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGaa      >  1:550102/1‑140 (MQ=255)                                                                                                                                   aTCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGaa      <  1:445169/140‑1 (MQ=255)                                                                                                                                    tCACTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAAc     <  2:550102/140‑1 (MQ=255)                                                                                                                                      aCTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACaa   >  2:582450/1‑140 (MQ=255)                                                                                                                                      aCTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACaa   >  2:754084/1‑140 (MQ=255)                                                                                                                                       cTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAg  <  1:754084/140‑1 (MQ=255)                                                                                                                                       cTGGCCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAg  <  1:724026/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         TAGCCGCCGTGAATTTCTGAGCCAGAGCGGTAAGATGGTCACCGCCGCCGCGCTGTTTGGTACCTCTGTGCCGCTCGCCCATGCGGCGGTAGCTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAG  >  NC_000913/341139‑341412

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