Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,380	C→A	G6C (GGC→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,380	0	C	A	100.0%	148.8 / NA	40	G6C (GGC→TGC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (26/14);  total (26/14)

TGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATA  >  NC_000913/4128236‑4128527                                                                                                                                                 |                                                                                                                                                    tGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTcc                                                                                                                                                >  1:603700/1‑150 (MQ=255) tGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTcc                                                                                                                                                <  2:119680/150‑1 (MQ=255)      gCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAg                                                                                                                                          <  2:345409/151‑1 (MQ=255)            cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCAt                                                                                                                                      >  2:435722/1‑149 (MQ=255)               aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGAt                                                                                                                                 <  1:131859/151‑1 (MQ=255)               aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGAt                                                                                                                                 <  2:93886/151‑1 (MQ=255)                        aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATcc                                                                                                                        >  2:521521/1‑151 (MQ=255)                         ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCct                                                                                                                       >  1:317369/1‑151 (MQ=255)                            gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCtt                                                                                                                    <  1:435781/151‑1 (MQ=255)                                          tCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                                                                                       <  1:657669/150‑1 (MQ=255)                                            aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                                                                                                    <  1:44746/151‑1 (MQ=255)                                            aTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                                                                                       >  1:677996/1‑148 (MQ=255)                                                 tGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTACGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgaga                                                                                                    >  2:211448/1‑146 (MQ=255)                                                   aGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagacc                                                                                             >  2:626382/1‑151 (MQ=255)                                                       tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGATTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                                                                                         >  1:387691/1‑151 (MQ=255)                                                         tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacaca                                                                                        <  2:286377/150‑1 (MQ=255)                                                             acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg                                                                                    >  1:195909/1‑150 (MQ=255)                                                               acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                                                                                 >  2:724272/1‑151 (MQ=255)                                                                         gAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTg                                                                        >  1:157595/1‑150 (MQ=255)                                                                           aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc                                                                     >  2:178741/1‑151 (MQ=255)                                                                             ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                                                   <  2:673759/151‑1 (MQ=255)                                                                              gAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                                                   >  2:582034/1‑150 (MQ=255)                                                                                       aTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCAt                                                         >  1:224510/1‑151 (MQ=255)                                                                                         cTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCAt                                                         >  1:356839/1‑149 (MQ=255)                                                                                           ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGt                                                      <  1:509889/150‑1 (MQ=255)                                                                                              tGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTAc                                                  >  1:319224/1‑151 (MQ=255)                                                                                                     ttCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTAc                                             >  1:786962/1‑149 (MQ=255)                                                                                                      tCACTCTTCTAGCCGTGCCCAGCGTATGGGCTGATATATTCGCAGCTCCATTCCGCCATGAGATACTTAATCCTCTTCGTCACTAAATTGAGACCAGACCACACAGTTGATGACGGTCCTGACTGTAAACCCGCATACTTTACCGTACAgg                                          >  2:349596/1‑151 (MQ=255)                                                                                                           cttcttGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATCCTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGtt                                      >  2:233968/1‑150 (MQ=255)                                                                                                              cttGCCGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCg                                  >  2:532348/1‑151 (MQ=255)                                                                                                                   cGTGCTCAGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACa                             <  1:247934/151‑1 (MQ=255)                                                                                                                         cAGCGTATGGGCTGATATATTCGTAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgtgt                       <  1:359895/151‑1 (MQ=255)                                                                                                                          aGCGTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgtgt                       >  2:324721/1‑150 (MQ=255)                                                                                                                             gTATGGGCTGATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATg                   >  1:649047/1‑151 (MQ=255)                                                                                                                                      gATATATTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCg          <  1:785928/151‑1 (MQ=255)                                                                                                                                       atatatTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCgg         <  1:352614/151‑1 (MQ=255)                                                                                                                                        tatatTCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGc        >  1:326395/1‑151 (MQ=255)                                                                                                                                            ttCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTg     >  1:28860/1‑150 (MQ=255)                                                                                                                                             tCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGAt   >  2:383185/1‑151 (MQ=255)                                                                                                                                             tCGCAGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGAt   >  2:241503/1‑151 (MQ=255)                                                                                                                                              cgcagcTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATa  <  2:604367/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATA  >  NC_000913/4128236‑4128527

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