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%	153.6 / NA	42	P104A (CCG→GCG)	yahJ	putative deaminase with metallo‑dependent hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (20/22);  total (20/22)

CGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAA  >  NC_000913/341357‑341501                                                                              |                                                                    cgcgcCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCg                                                                   >  1:66742/1‑80 (MQ=255)       aGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATg                                                            <  2:66677/81‑1 (MQ=255)          ccGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACg                                                         <  1:66678/81‑1 (MQ=255)             tGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTg                                                      <  2:66679/81‑1 (MQ=255)              ggAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTgg                                                     >  1:66741/1‑81 (MQ=255)               gagaTTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGc                                                    <  2:66680/81‑1 (MQ=255)                agaTTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCg                                                   <  1:66681/81‑1 (MQ=255)                     cAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTaa                                               <  1:66682/80‑1 (MQ=255)                     cAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTaa                                               <  1:66683/80‑1 (MQ=255)                           ggCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtg                                        >  1:66748/1‑81 (MQ=255)                           ggCAAAATTCTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtg                                        >  1:66734/1‑81 (MQ=255)                             cAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtgct                                      <  2:66684/81‑1 (MQ=255)                             cAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtgct                                      >  2:66727/1‑81 (MQ=255)                               aaaTTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtgctgc                                    >  2:66731/1‑81 (MQ=255)                               aaaTTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtgctgc                                    <  2:66685/81‑1 (MQ=255)                               aaaTTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGAtgctgc                                    <  1:66686/81‑1 (MQ=255)                                aaTTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGcc                                   <  2:66687/81‑1 (MQ=255)                                    gTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCccac                                <  2:66688/80‑1 (MQ=255)                                       gCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCcaccac                             <  1:66689/80‑1 (MQ=255)                                       gCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCaccacc                            >  1:66737/1‑81 (MQ=255)                                       gCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCaccacc                            >  2:66730/1‑81 (MQ=255)                                       gCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCaccacc                            >  1:66732/1‑81 (MQ=255)                                        cccTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCAccc                           <  2:66690/81‑1 (MQ=255)                                          cTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCcgc                         >  1:66749/1‑81 (MQ=255)                                           tGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCcgc                         >  2:66745/1‑80 (MQ=255)                                           tGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCcgc                         >  2:66744/1‑80 (MQ=255)                                            gcgcGAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGa                       >  2:66733/1‑81 (MQ=255)                                             cgcgAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCcgcg                        <  2:66691/79‑1 (MQ=255)                                              gcgAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACa                     >  1:66735/1‑81 (MQ=255)                                               cgAGAACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACAt                    <  1:66692/81‑1 (MQ=255)                                                gagaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATg                   <  2:66693/81‑1 (MQ=255)                                                gagaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATg                   <  2:66694/81‑1 (MQ=255)                                                 agaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGc                  <  2:66695/81‑1 (MQ=255)                                                 agaACAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGc                  <  2:66696/81‑1 (MQ=255)                                                   aaCAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCGGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCat                <  1:66697/81‑1 (MQ=255)                                                    aCAAGCTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCata               <  2:66698/81‑1 (MQ=255)                                                         cTGCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAt          >  1:66750/1‑81 (MQ=255)                                                           gCATCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAtct        >  1:66736/1‑81 (MQ=255)                                                             aTCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAtctc       >  2:66761/1‑80 (MQ=255)                                                             aTCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAtctc       >  1:66743/1‑80 (MQ=255)                                                             aTCCGGACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAtctc       >  1:66728/1‑80 (MQ=255)                                                                  gACGCCACGCTGGCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACaa  >  1:66758/1‑80 (MQ=255)                                                                              |                                                                    CGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAA  >  NC_000913/341357‑341501

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