Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,303	T→C	M201T (ATG→ACG)	cmoA →	carboxy‑SAM synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,303	0	T	C	100.0%	153.7 / NA	42	M201T (ATG→ACG)	cmoA	carboxy‑SAM synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/26);  total (16/26)

AAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGA  >  NC_000913/1953156‑1953445                                                                                                                                                    |                                                                                                                                               aaGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCt                                                                                                                                             <  2:237228/151‑1 (MQ=255)  aGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTg                                                                                                                                            >  1:672815/1‑151 (MQ=255)      cTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGaa                                                                                                                                         >  1:553539/1‑150 (MQ=255)       tGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGaaa                                                                                                                                        >  1:105407/1‑150 (MQ=255)         aaCCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACg                                                                                                                                     <  1:18043/151‑1 (MQ=255)         aaCCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACg                                                                                                                                     >  1:394604/1‑151 (MQ=255)              cggcggCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATg                                                                                                                                <  2:38253/151‑1 (MQ=255)                gcggcgCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCt                                                                                                                              <  2:224222/151‑1 (MQ=255)                 cggcgCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTg                                                                                                                             >  1:219224/1‑151 (MQ=255)                     gcgcTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCg                                                                                                                         <  2:744486/151‑1 (MQ=255)                         tGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTc                                                                                                                     >  2:648355/1‑151 (MQ=255)                            tGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGt                                                                                                                  >  2:333264/1‑151 (MQ=255)                                      aaaaTTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCAt                                                                                                        >  1:431062/1‑151 (MQ=255)                                           tCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGc                                                                                                   >  1:595636/1‑151 (MQ=255)                                             aGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCAc                                                                                                 >  1:400804/1‑151 (MQ=255)                                                   gAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGc                                                                                           >  1:220332/1‑151 (MQ=255)                                                     aGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCAt                                                                                         <  2:516326/151‑1 (MQ=255)                                                     aGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCAt                                                                                         <  1:706560/151‑1 (MQ=255)                                                      gATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATa                                                                                        <  1:188668/151‑1 (MQ=255)                                                       aTGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATaa                                                                                       <  1:601333/151‑1 (MQ=255)                                                              aGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGt                                                                                >  1:239774/1‑151 (MQ=255)                                                                                 aaCATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCtgt                                                             <  1:604252/151‑1 (MQ=255)                                                                                    aTGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGt                                                          <  1:333231/151‑1 (MQ=255)                                                                                       caccacGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTcc                                                       <  2:637599/151‑1 (MQ=255)                                                                                         ccacGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTcc                                                       >  2:568030/1‑149 (MQ=255)                                                                                           acGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTg                                                   >  1:449877/1‑151 (MQ=255)                                                                                              aCTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCtt                                                <  1:352495/151‑1 (MQ=255)                                                                                                 ttAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCttt                                               <  1:691103/149‑1 (MQ=255)                                                                                                 ttAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTa                                              <  1:77807/150‑1 (MQ=255)                                                                                                        tGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGt                                      >  2:312414/1‑151 (MQ=255)                                                                                                                  tACAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACtggtg                             <  1:466672/150‑1 (MQ=255)                                                                                                                   aCAGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGc                           <  2:165891/151‑1 (MQ=255)                                                                                                                     aGCGAACTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCAt                         <  1:77209/151‑1 (MQ=255)                                                                                                                         aaCTGGAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTaaa                     <  1:125921/151‑1 (MQ=255)                                                                                                                             ggAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAg                   <  1:504881/149‑1 (MQ=255)                                                                                                                             ggAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCa                 <  2:317402/151‑1 (MQ=255)                                                                                                                             ggAGATCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCa                 <  1:256255/151‑1 (MQ=255)                                                                                                                                gaTCAGCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCagag              <  2:553591/151‑1 (MQ=255)                                                                                                                                     gCCAGAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGc         <  2:212271/151‑1 (MQ=255)                                                                                                                                         gAAGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGc      >  1:481347/1‑150 (MQ=255)                                                                                                                                           aGCGCAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATg   <  2:275156/151‑1 (MQ=255)                                                                                                                                            gcgcAGCACGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGa  <  1:416459/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               AAGGGCTGAACCCCGGCGGCGCGCTGGTGCTTTCGGAAAAATTCAGTTTCGAAGATGCCAAAGTTGGTGAACTGTTGTTCAACATGCACCACGACTTTAAACGTGCCAACGGTTACAGCGAACTGGAGATCAGCCAGAAGCGCAGCATGCTGGAAAACGTGATGCTGACCGATTCCGTGGAAACCCATAAAGCACGCCTGCATAACGCCGGTTTTGAGCATAGTGAGCTGTGGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGA  >  NC_000913/1953156‑1953445

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