Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	341,378	A→G	Q85R (CAG→CGG)	yahJ →	putative metallo‑dependent hydrolase domain deaminase

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

CTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTTTACGGCGGG  >  NC_000913/341231‑341517                                                                                                                                                    |                                                                                                                                            cTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCggg                                                                                                                                           >  2:812615/1‑150 (MQ=255)   ggCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGa                                                                                                                                          <  2:72464/149‑1 (MQ=255)   ggCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGAc                                                                                                                                         <  1:179477/150‑1 (MQ=255)    gCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACgg                                                                                                                                       <  2:279170/151‑1 (MQ=255)                   aagcgaaCAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCt                                                                                                                        <  1:503719/151‑1 (MQ=255)                           aacaCCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTgcgc                                                                                                                    >  1:66758/1‑147 (MQ=255)                           aacaCCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCgaga                                                                                                                <  1:401125/151‑1 (MQ=255)                              acCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACa                                                                                                             <  2:74784/151‑1 (MQ=255)                                   gAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCt                                                                                                         <  1:434468/150‑1 (MQ=255)                                       aTCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGTCGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATc                                                                                                    >  2:289373/1‑151 (MQ=255)                                             gACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGGCGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGAc                                                                                               >  2:582623/1‑150 (MQ=255)                                             gACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGAc                                                                                               <  2:195106/150‑1 (MQ=255)                                             gACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGAc                                                                                               <  2:78434/150‑1 (MQ=255)                                                 cGCATTACTATCTCGATAACGCGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCa                                                                                           >  1:328013/1‑150 (MQ=255)                                                    aTTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCt                                                                                       <  1:15377/151‑1 (MQ=255)                                                          aTCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCg                                                                                   <  1:148863/149‑1 (MQ=255)                                                                  aaCGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGAc                                                                          >  2:292177/1‑150 (MQ=255)                                                                  aaCGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACCGCGCGCCAGAACGGGGGGATTCGGGACGGCAAAATTGTTGCCCTGCGCGAGAAAAAGCTGCATCCGGACGCCAAGCTTCCCCACTATggc                                                                          >  1:736364/1‑148 (MQ=255)                                                                     gTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTg                                                                      >  2:646920/1‑151 (MQ=255)                                                                       gctgctGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTgg                                                                     >  2:507009/1‑150 (MQ=255)                                                                             ggAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTa                                                                <  2:521890/149‑1 (MQ=255)                                                                                  ccGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCATGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGAt                                                         <  2:196003/151‑1 (MQ=255)                                                                                     gTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTTGCGGTAAGCTGAtgct                                                      <  1:766508/151‑1 (MQ=255)                                                                                     gTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGAtgct                                                      <  2:827043/151‑1 (MQ=255)                                                                                       tttGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGAACGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGGAAGc                                                             >  1:16016/1‑142 (MQ=255)                                                                                         tGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAACTGATGCCGccc                                                  >  1:88752/1‑151 (MQ=255)                                                                                           aCTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCcca                                                 >  2:100342/1‑150 (MQ=255)                                                                                            cTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCccacc                                               >  2:755711/1‑151 (MQ=255)                                                                                            cTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGCGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCCGACGCCACGCTGCCGCACTATGACGCTGGCGGGAAGCTGATGCTTCccacc                                               >  1:189775/1‑151 (MQ=255)                                                                                                 aaaaTGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCAccc                                           >  2:824811/1‑150 (MQ=255)                                                                                                 aaaaTGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCCGCCCACCAccc                                           >  2:628632/1‑150 (MQ=255)                                                                                                     tggcgtggcgGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGa                                       >  2:227389/1‑150 (MQ=255)                                                                                                      ggcgtggcgGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGAc                                      >  2:855226/1‑150 (MQ=255)                                                                                                             cgGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCatat                              <  2:323122/151‑1 (MQ=255)                                                                                                              ggTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCCTAtt                             >  2:698787/1‑151 (MQ=255)                                                                                                              ggTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATAtt                             <  2:498273/151‑1 (MQ=255)                                                                                                                  cAGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAtc                         >  1:368280/1‑151 (MQ=255)                                                                                                                   aGACCCGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCAtct                        >  2:674857/1‑151 (MQ=255)                                                                                                                      cccGCACGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATTTCg                      >  1:239295/1‑150 (MQ=255)                                                                                                                            cGGCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACaaaa                <  1:358628/150‑1 (MQ=255)                                                                                                                             ggCGCGCCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAAcc              <  2:777911/151‑1 (MQ=255)                                                                                                                                  gcCAGACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCTCGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAACCCttt           <  2:186648/149‑1 (MQ=255)                                                                                                                                      gACCGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTTTAcggc     <  1:451629/151‑1 (MQ=255)                                                                                                                                        ccGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTTTAcggcgg   <  1:26964/151‑1 (MQ=255)                                                                                                                                         cGTGGAGATTCGGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTTTACGGCggg  <  2:49666/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                            CTGGCACCCTAAACTGCGAAGCGAACAACACCATGAAAATCACTGACCCGCATTACTATCTCGATAACGTGCTGCTGGAAACCGGTTTTGACTACGAAAATGGCGTGGCGGTGCAGACCCGCACGGCGCGCCAGACCGTGGAGATTCAGGACGGCAAAATTGTCGCCCTGCGCGAGAACAAGCTGCATCCGGACGCCACGCTGCCGCACTATGACGCTGGCGGTAAGCTGATGCTGCCCACCACCCGCGACATGCATATTCATCTCGACAAAACCTTTTACGGCGGG  >  NC_000913/341231‑341517

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