Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,496,196	A→G	R303R (CGA→CGG)	ydcJ →	putative metalloenzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,496,196	0	A	G	100.0%	181.6 / NA	48	R303R (CGA→CGG)	ydcJ	putative metalloenzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/31);  total (17/31)

CGAGGTACCGATTTTACTACGCCAGACCAGCTTTAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTAC  >  NC_000913/1496055‑1496337                                                                                                                                              |                                                                                                                                              cgAGGTACCGATTTTACTACGCCAGACCAGCTTTAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTat                                                                                                                                      <  2:395637/151‑1 (MQ=255)   aGGTACCGATTTTACTACGCCAGACCAGCTTTAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTatg                                                                                                                                     >  2:620242/1‑150 (MQ=255)                   acGCCAGACCAGCTTTAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTaa                                                                                                                     <  2:361184/150‑1 (MQ=255)                       cAGACCAGCTTTAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTg                                                                                                                <  2:73400/151‑1 (MQ=255)                                 ttAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCggg                                                                                                        <  1:389759/149‑1 (MQ=255)                                   aaaGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAgg                                                                                                    >  2:761176/1‑151 (MQ=255)                                   aaaGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAgg                                                                                                    <  2:592939/151‑1 (MQ=255)                                     aGCACTGGGAGAGACGGGGTTGTTTGCGGGGCAGAAACAGGGCACGCATGCCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATGAACGCCGAAAGGGCGGCATCTGCATGATGATCTTCTGCGTAACGCTGGATCCGGGCAGGAt                                                                                                  <  2:415231/151‑1 (MQ=255)                                      gCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATa                                                                                                 <  2:690205/151‑1 (MQ=255)                                        aCTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAAt                                                                                               <  1:196130/151‑1 (MQ=255)                                           ggAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAAtctc                                                                                            <  1:113969/151‑1 (MQ=255)                                           ggAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAAtctc                                                                                            <  2:248655/151‑1 (MQ=255)                                           ggAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAAtct                                                                                             <  1:619108/150‑1 (MQ=255)                                              agagACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATctcact                                                                                         >  1:610186/1‑151 (MQ=255)                                               gagaCGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATctcactc                                                                                        <  2:703439/151‑1 (MQ=255)                                                     gtgtTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCaaa                                                                                  <  2:528038/151‑1 (MQ=255)                                                       gttgttTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATg                                                                                >  1:187818/1‑151 (MQ=255)                                                            ttGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCAttt                                                                           >  2:109583/1‑151 (MQ=255)                                                             tGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCAttt                                                                           >  2:501937/1‑150 (MQ=255)                                                                   gCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGa                                                                     >  1:13135/1‑150 (MQ=255)                                                                      gAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAAcc                                                                 >  2:617327/1‑151 (MQ=255)                                                                          cAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTcc                                                             <  1:181984/151‑1 (MQ=255)                                                                                    aTACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTcc                                                   <  2:334316/151‑1 (MQ=255)                                                                                      aCCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTg                                                 >  1:173532/1‑151 (MQ=255)                                                                                      aCCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTg                                                 <  1:750770/151‑1 (MQ=255)                                                                                      aCCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTg                                                 <  1:519355/151‑1 (MQ=255)                                                                                       ccGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTg                                                 >  1:72678/1‑150 (MQ=255)                                                                                       ccGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTg                                                 >  2:223031/1‑150 (MQ=255)                                                                                       ccGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGa                                                <  2:472812/151‑1 (MQ=255)                                                                                       ccGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGa                                                <  1:606914/151‑1 (MQ=255)                                                                                        cgcgcgCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGAc                                               <  1:432615/151‑1 (MQ=255)                                                                                          cgcgcTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAg                                             <  1:702809/151‑1 (MQ=255)                                                                                            cgcTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTg                                           >  1:423834/1‑151 (MQ=255)                                                                                            cgcTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTg                                           >  1:26135/1‑151 (MQ=255)                                                                                                  ggTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGtttt                                     >  1:33413/1‑151 (MQ=255)                                                                                                     gAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTa                                   <  1:300740/150‑1 (MQ=255)                                                                                                     gAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTa                                   <  2:753/150‑1 (MQ=255)                                                                                                        aTTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGc                               <  2:59833/151‑1 (MQ=255)                                                                                                           gAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGt                             <  1:308520/150‑1 (MQ=255)                                                                                                            agcagcGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCa                           <  2:103298/151‑1 (MQ=255)                                                                                                             gcagcGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAg                          <  1:708522/151‑1 (MQ=255)                                                                                                               agcGTGGCGTGGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGc                         >  1:690923/1‑150 (MQ=255)                                                                                                                gcgtggcgtgGCATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCaa                       <  1:361504/151‑1 (MQ=255)                                                                                                                    ggcgtggcATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGt                   <  2:171941/151‑1 (MQ=255)                                                                                                                    ggcgtggcATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGt                   >  2:512212/1‑151 (MQ=255)                                                                                                                     gcgtggcATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGtt                  >  2:71777/1‑151 (MQ=255)                                                                                                                         ggcATTAACGCCGAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCa              <  2:535863/151‑1 (MQ=255)                                                                                                                                     gAAAGGGCGGCAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTAc  <  2:343220/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                              CGAGGTACCGATTTTACTACGCCAGACCAGCTTTAAAGCACTGGAAGAGACGGTGTTGTTTGCGGGGCAGAAACAGGGCACGCATACCGCGCGCTTTGGTGAAATTGAGCAGCGTGGCGTGGCATTAACGCCGAAAGGGCGACAACTGTATGATGATCTTCTGCGTAACGCTGGAACCGGGCAGGATAATCTCACTCACCAAATGCATTTACAGGAAACCTTCCGCACTTTTCCTGACAGTGAGTTTTTAATGCGTCAGCAAGGGTTGGCATGGTTCCGGTAC  >  NC_000913/1496055‑1496337

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