Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,379	C→T	G6D (GGC→GAC)	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,379	0	C	T	100.0%	163.4 / NA	42	G6D (GGC→GAC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/23);  total (19/23)

GTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGAT  >  NC_000913/4128234‑4128526                                                                                                                                                  |                                                                                                                                                    gTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCt                                                                                                                                                 <  2:235503/150‑1 (MQ=255)    ggTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCAt                                                                                                                                             >  2:504791/1‑151 (MQ=255)     gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCAtt                                                                                                                                            <  1:404529/151‑1 (MQ=255)       aGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCa                                                                                                                                          <  2:553367/151‑1 (MQ=255)               gCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATgag                                                                                                                                  <  2:20022/151‑1 (MQ=255)                   gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATAc                                                                                                                              <  1:312050/151‑1 (MQ=255)                      gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACtt                                                                                                                            <  2:759961/150‑1 (MQ=255)                      gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTa                                                                                                                           <  2:79989/151‑1 (MQ=255)                        tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTaa                                                                                                                          <  2:705859/150‑1 (MQ=255)                        tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTaa                                                                                                                          <  1:91821/150‑1 (MQ=255)                          aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATcc                                                                                                                       >  2:260917/1‑151 (MQ=255)                           ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCct                                                                                                                      <  2:620/151‑1 (MQ=255)                             tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCctct                                                                                                                    >  1:321164/1‑151 (MQ=255)                               acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGTGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTc                                                                                                                  <  2:2710/151‑1 (MQ=255)                                          gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                                                                                       <  2:341634/151‑1 (MQ=255)                                            tCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                                                                                      >  2:866873/1‑150 (MQ=255)                                                cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                                                                                 >  2:130261/1‑151 (MQ=255)                                                   tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                                                                                 >  1:788113/1‑148 (MQ=255)                                                     aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                                                                                            <  2:437293/151‑1 (MQ=255)                                                     aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                                                                                            <  2:195833/151‑1 (MQ=255)                                                      ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                                                                                            >  2:545710/1‑150 (MQ=255)                                                              aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATAATTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg                                                                                   >  1:311523/1‑151 (MQ=255)                                                                  ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTACACAGCGTATGGGCTTATATACTCGTCGCAACAAACAGCCATGAGAGACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgt                                                                               <  1:636448/151‑1 (MQ=255)                                                                     tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg                                                                            >  2:671413/1‑151 (MQ=255)                                                                            aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAc                                                                     >  1:17234/1‑151 (MQ=255)                                                                               ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                                                  >  2:861291/1‑151 (MQ=255)                                                                                gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                                                  >  2:470749/1‑150 (MQ=255)                                                                                         aTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCAt                                                        >  1:310150/1‑151 (MQ=255)                                                                                             ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCATCGATCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACAAGACCACACAGTTGATGTGGGTACTGACCGTAAAACCGCATa                                                       <  1:220115/148‑1 (MQ=255)                                                                                                        tCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACa                                           <  2:926768/149‑1 (MQ=255)                                                                                                           ctctTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCGTTCAGCCATGAGATACTTACTCCTCGTCGTCAATAAATTGAGACTTGACCAAAGAGTTGATGTGGGTACTGTCCGTAAACCCACATATTTTACCGTACAggaa                                       <  1:456254/150‑3 (MQ=255)                                                                                                             cttcttGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGt                                      <  1:66579/149‑1 (MQ=255)                                                                                                                 ttGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGt                                >  2:146061/1‑151 (MQ=255)                                                                                                                  tGCCGTGCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTAcc                                  <  1:967785/148‑1 (MQ=255)                                                                                                                        gCTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATc                          <  2:942871/150‑1 (MQ=255)                                                                                                                         cTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgt                        >  2:683215/1‑151 (MQ=255)                                                                                                                         cTCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgt                        >  1:402365/1‑151 (MQ=255)                                                                                                                          tCAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgtg                       >  1:21856/1‑151 (MQ=255)                                                                                                                           cAGCGTATGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgtg                       <  2:823185/150‑1 (MQ=255)                                                                                                                            aGCGTATGGGCTGATATATTCGTAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgtgt                      >  1:225316/1‑150 (MQ=255)                                                                                                                                  tGGGCTGATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCAc               >  1:47049/1‑151 (MQ=255)                                                                                                                                        gATATATTCGTCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCg         <  2:310192/151‑1 (MQ=255)                                                                                                                                               tcgtcgCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGAt  >  1:86000/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGAT  >  NC_000913/4128234‑4128526

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: ‑