Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,361	T→C	Y12C (TAC→TGC)	metJ ←	transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,361	0	T	C	100.0%	217.2 / NA	57	Y12C (TAC→TGC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/23);  total (34/23)

CGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATC  >  NC_000913/4128213‑4128502                                                                                                                                                     |                                                                                                                                              cGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCAt                                                                                                                                             >  1:117899/1‑151 (MQ=255)     tCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGc                                                                                                                                         >  2:488064/1‑151 (MQ=255)         acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGat                                                                                                                                     <  2:139805/151‑1 (MQ=255)         acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGat                                                                                                                                     >  1:731493/1‑151 (MQ=255)             cagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGatatat                                                                                                                                 >  1:156027/1‑151 (MQ=255)              agcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATAtt                                                                                                                                <  2:625878/151‑1 (MQ=255)              accTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATAtt                                                                                                                                >  2:146032/3‑151 (MQ=255)               gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTc                                                                                                                               <  2:354794/151‑1 (MQ=255)               gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTc                                                                                                                               <  2:674148/151‑1 (MQ=255)               gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTc                                                                                                                               <  2:113917/151‑1 (MQ=255)                cTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTc                                                                                                                               >  1:201515/1‑150 (MQ=255)                 tCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTcgc                                                                                                                             >  2:306263/1‑151 (MQ=255)                         ggTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCAt                                                                                                                     >  1:208883/1‑151 (MQ=255)                          gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCAtt                                                                                                                    <  1:709882/151‑1 (MQ=255)                               gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                                                                               >  2:97196/1‑151 (MQ=255)                               gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                                                                               >  1:18514/1‑151 (MQ=255)                               gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGc                                                                                                                <  2:380583/150‑1 (MQ=255)                                   cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                                                                            >  1:283088/1‑150 (MQ=255)                                   cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATg                                                                                                            >  2:299600/1‑150 (MQ=255)                                             tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                                                                  <  2:356339/150‑1 (MQ=255)                                               aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATc                                                                                                <  2:466278/150‑1 (MQ=255)                                                      gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                                                                        <  1:59556/151‑1 (MQ=255)                                                      gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                                                                        >  1:166147/1‑151 (MQ=255)                                                      gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                                                                        >  1:430356/1‑151 (MQ=255)                                                        cgcgTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                                                                      <  2:382146/151‑1 (MQ=255)                                                         gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCaa                                                                                     <  2:659184/151‑1 (MQ=255)                                                           gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATa                                                                                   <  1:592334/151‑1 (MQ=255)                                                            tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATa                                                                                   >  2:383546/1‑150 (MQ=255)                                                               gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                                                               >  1:304128/1‑151 (MQ=255)                                                                    tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                                                                          >  2:235753/1‑151 (MQ=255)                                                                         gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga                                                                      <  2:206083/150‑1 (MQ=255)                                                                            gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac                                                                  >  2:678557/1‑151 (MQ=255)                                                                              tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                                                                >  1:530985/1‑151 (MQ=255)                                                                                 ttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt                                                             >  1:389530/1‑151 (MQ=255)                                                                                   aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt                                                            <  1:202584/150‑1 (MQ=255)                                                                                    acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg                                                           >  2:616954/1‑150 (MQ=255)                                                                                        cTTAAGAGGAATGGAAACCGTTATCTTTTTGACTTGTTCCCTCCTCCTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                                                        <  2:270659/149‑1 (MQ=255)                                                                                                gAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTg                                               >  1:107441/1‑150 (MQ=255)                                                                                                       aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTaa                                        >  2:22292/1‑150 (MQ=255)                                                                                                        aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAc                                      <  1:143350/151‑1 (MQ=255)                                                                                                          cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAccc                                    <  1:265759/151‑1 (MQ=255)                                                                                                          cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAccc                                    <  1:59227/151‑1 (MQ=255)                                                                                                               tCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGc                                  >  2:286340/1‑148 (MQ=255)                                                                                                                cTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCAt                                >  1:426103/1‑149 (MQ=255)                                                                                                                    ttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAATCCCGCATAGttt                           <  1:668443/150‑1 (MQ=255)                                                                                                                    ttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGttt                           >  1:303350/1‑150 (MQ=255)                                                                                                                            ttCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACa                   <  2:437916/150‑1 (MQ=255)                                                                                                                              cACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGc                >  2:632785/1‑151 (MQ=255)                                                                                                                                ctctTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGt              >  2:8606/1‑151 (MQ=255)                                                                                                                                  cttcttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGtt             >  1:363087/1‑150 (MQ=255)                                                                                                                                  cttcttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGtt             >  2:31348/1‑150 (MQ=255)                                                                                                                                   ttcttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGt              >  2:635382/1‑148 (MQ=255)                                                                                                                                    tcttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGtt             <  2:97569/148‑1 (MQ=255)                                                                                                                                        gCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTg       <  1:389092/150‑1 (MQ=255)                                                                                                                                         ccGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGAc     >  1:196565/1‑151 (MQ=255)                                                                                                                                         ccGTGCTCAGCGCATGGGCAGATATATTCGCCGCTCCATTCAGCCATGAGAGACTTAATCCTCTTCGTCAATAAATTGAGACCAGACAACACATTTGATGTGGGTACTGACTGTAAACACGCATAGTTTACCGTACAGGCGTTACCGTGa      >  2:588796/1‑150 (MQ=255)                                                                                                                                            ggcgcAGCGCATGGGCTGATATATTCGCAGCACGACCCAGCCATGAGATACTTAATCCTCTCCGACAATAAAAGGAGACCAGACCACACAGAGGATGAGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGCTACCGTGACATc  >  2:721774/5‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                              CGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATC  >  NC_000913/4128213‑4128502

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