Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	341,434	C→G	P104A (CCG→GCG)	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,434	0	C	G	100.0%	132.5 / NA	38	P104A (CCG→GCG)	yahJ	putative deaminase with metallo‑dependent hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (14/24);  total (14/24)

GCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTT  >  NC_000913/341360‑341508                                                                           |                                                                           gcCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCAc                                                                       >  2:70633/1‑80 (MQ=255) gcCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACt                                                                      <  2:70584/81‑1 (MQ=255)   cAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTAt                                                                    <  1:70585/81‑1 (MQ=255)   cAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTAt                                                                    >  1:70625/1‑81 (MQ=255)        cGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGc                                                               <  2:70586/81‑1 (MQ=255)         gTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCt                                                              <  2:70587/81‑1 (MQ=255)          tGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTg                                                             >  2:70643/1‑81 (MQ=255)          tGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTg                                                             <  2:70588/81‑1 (MQ=255)            gagaTTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGc                                                           >  2:70631/1‑81 (MQ=255)              gaTTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCgg                                                         <  1:70589/81‑1 (MQ=255)              gaTTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCgg                                                         >  1:70655/1‑81 (MQ=255)               aTTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGt                                                        <  2:70590/81‑1 (MQ=255)                   aGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGc                                                    <  1:70591/81‑1 (MQ=255)                   aGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGc                                                    <  1:70592/81‑1 (MQ=255)                       cGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAt                                                >  1:70648/1‑81 (MQ=255)                       cGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAt                                                <  1:70593/81‑1 (MQ=255)                                 gTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCccac                                       <  2:70594/80‑1 (MQ=255)                                  tCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCccacc                                      <  2:70595/80‑1 (MQ=255)                                     cccTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCAccc                                  >  2:70642/1‑81 (MQ=255)                                     cccTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCAccc                                  <  1:70596/81‑1 (MQ=255)                                      ccTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCcg                                 <  1:70597/81‑1 (MQ=255)                                          cgcgAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGa                              >  1:70640/1‑80 (MQ=255)                                            cgAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACAt                           <  2:70599/81‑1 (MQ=255)                                            cgAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACAt                           <  1:70598/81‑1 (MQ=255)                                             gagaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATg                          <  1:70600/81‑1 (MQ=255)                                             gagaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATg                          >  1:70647/1‑81 (MQ=255)                                             gagaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATg                          >  2:70645/1‑81 (MQ=255)                                               gaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCa                        >  2:70654/1‑81 (MQ=255)                                               gaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCa                        >  1:70646/1‑81 (MQ=255)                                                   aaGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATAtt                    <  1:70601/81‑1 (MQ=255)                                                   aaGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATAtt                    <  1:70602/81‑1 (MQ=255)                                                    aGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATAtt                    >  2:70651/1‑80 (MQ=255)                                                             cGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGTGACATGCATATTCATCTCGAc           <  2:70603/80‑1 (MQ=255)                                                             cGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACa          <  1:70604/81‑1 (MQ=255)                                                              ggACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACaa         <  2:70605/81‑1 (MQ=255)                                                                   ccACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCt    <  2:70606/81‑1 (MQ=255)                                                                   ccACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCt    <  1:70607/81‑1 (MQ=255)                                                                     aCGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCttt  >  1:70644/1‑81 (MQ=255)                                                                           |                                                                           GCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTT  >  NC_000913/341360‑341508

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