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%	147.4 / NA	40	Y12C (TAC→TGC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/22);  total (18/22)

TCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGT  >  NC_000913/4128217‑4128504                                                                                                                                                 |                                                                                                                                                tcccacaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGCGCGATCACCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGc                                                                                                                                           >  1:111162/4‑151 (MQ=255)     acaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGat                                                                                                                                       <  2:240520/151‑1 (MQ=255)        gcagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGatat                                                                                                                                     <  1:76452/150‑1 (MQ=255)          agcTCGCTGTTGGTAGCGTGACGCAGTTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGATGAATGGAAACCGTAAACTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATAtt                                                                                                                                  <  1:635615/151‑1 (MQ=255)           gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTc                                                                                                                                 >  2:358443/1‑151 (MQ=255)           gcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTc                                                                                                                                 <  2:128684/151‑1 (MQ=255)              cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCTTGCGTTCATCGGTGATGATTTTTAACACCGTAAGAGGAATGGAAACCGTAGTCTTTTTGATTTGTTCGCTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTcgcc                                                                                                                              <  1:542524/151‑1 (MQ=255)                  gTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCt                                                                                                                           <  2:451161/150‑1 (MQ=255)                                       gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACtt                                                                                                      <  2:386334/150‑1 (MQ=255)                                                acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCtt                                                                                             <  2:387727/150‑1 (MQ=255)                                                acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTc                                                                                            <  1:706944/151‑1 (MQ=255)                                                  gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCg                                                                                           <  1:331446/150‑1 (MQ=255)                                                  gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                                                                          <  1:157570/151‑1 (MQ=255)                                                     gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCa                                                                                        >  1:450788/1‑150 (MQ=255)                                                      cgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                                      <  1:239509/151‑1 (MQ=255)                                                        tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGATGAATGGAAACCGTAATCTTTTTGTCTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATAGATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                                      <  1:352621/149‑1 (MQ=255)                                                                  ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacca                                                                          <  1:715607/151‑1 (MQ=255)                                                                  ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacca                                                                          <  2:21016/151‑1 (MQ=255)                                                                  ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                                                           >  2:449576/1‑150 (MQ=255)                                                                        gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac                                                                    >  2:351651/1‑151 (MQ=255)                                                                                acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGa                                                            <  1:386024/151‑1 (MQ=255)                                                                                 cacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt                                                           <  1:767568/151‑1 (MQ=255)                                                                                     ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCAAGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg                                                       >  2:216115/1‑151 (MQ=255)                                                                                      tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg                                                      >  2:338994/1‑151 (MQ=255)                                                                                        agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTa                                                    <  2:62969/151‑1 (MQ=255)                                                                                              aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc                                              <  2:663488/151‑1 (MQ=255)                                                                                              aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc                                              <  2:861/151‑1 (MQ=255)                                                                                                tgAACCCGTAATCTTTTTGACATGCTCACTCTTCTTGCCGTGATCAGAGCATGGGCTGATATATTCGCCGCTACATTCAGCCATGAGAGACTTAATCCTCTTCGGCAATAAATTGCGACCAGACCAAACCGTTGATGAGGGTACTGTCCGt                                            >  1:559700/2‑151 (MQ=255)                                                                                                      cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAAccc                                      >  2:260996/1‑151 (MQ=255)                                                                                                         aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGc                                    >  1:278744/1‑150 (MQ=255)                                                                                                               ttttACTTGTTCACGCTTCTGTCCATCCTCAGCGCATGGCGTGATATAGTCGACGCTCCATTCATCCATGAGATACTTATTCCTCTTCTACAATAAAATGAGACCAGACCCCACAGTAGATGTGGTTGCTGACCGTAAACCCGCATa                                 <  1:705838/147‑1 (MQ=255)                                                                                                               tttGACTTGTTCCCTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCACTAGATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCCTAgct                              >  1:779887/1‑148 (MQ=255)                                                                                                                    cTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGt                        >  2:443125/1‑151 (MQ=255)                                                                                                                          cACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGc                  >  2:245793/1‑151 (MQ=255)                                                                                                                            ctctTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGt                >  1:783506/1‑151 (MQ=255)                                                                                                                             tcttctTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGtt               >  2:672326/1‑151 (MQ=255)                                                                                                                                 cttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCg           >  1:490788/1‑151 (MQ=255)                                                                                                                                   tGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTg         >  1:154386/1‑151 (MQ=255)                                                                                                                                         gCTCAGCGCGTGGGCGGATATATTCGCTGCTACATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCg   >  1:398939/1‑151 (MQ=255)                                                                                                                                            cAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgt  <  2:637982/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                TCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGT  >  NC_000913/4128217‑4128504

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