Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,226,509	C→T	*254* (TAG→TAA)	yohF ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,226,509	0	C	T	100.0%	128.6 / NA	41	*254* (TAG→TAA)	yohF	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/24);  total (17/24)

GGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCG  >  NC_000913/2226398‑2226630                                                                                                                |                                                                                                                          ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                               <  1:115350/124‑1 (MQ=255) ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                               >  2:115350/1‑124 (MQ=255)          gccgcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                                       <  2:750040/139‑1 (MQ=255)             gcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                    <  2:484737/139‑1 (MQ=255)             gcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                    <  2:128264/139‑1 (MQ=255)             gcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                    <  2:713591/139‑1 (MQ=255)                 ccGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCAt                                                                                <  1:293886/139‑1 (MQ=255)                         cGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCAGCTCACTGCGTTCTTCTGGGTTGAACTATGGATGCGCCAACATAAAGCCGCCATCCACTATc                                                                        <  1:517557/139‑1 (MQ=255)                            aGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAAc                                                                     <  2:715285/139‑1 (MQ=255)                             gTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGTCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccgc                                                                                   >  1:708417/1‑124 (MQ=255)                             gTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGTCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccgc                                                                                   <  2:708417/124‑1 (MQ=255)                             gTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGATTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACg                                                                    >  1:543636/1‑139 (MQ=255)                                 cATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTg                                                                <  2:665037/139‑1 (MQ=255)                                   tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                              <  2:256958/139‑1 (MQ=255)                                   tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                              >  1:712778/1‑139 (MQ=255)                                   tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                              <  2:712778/139‑1 (MQ=255)                                        gATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAGTATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtg                                                         >  1:494166/1‑139 (MQ=255)                                                   tgtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgc                                              >  1:515351/1‑139 (MQ=255)                                                    gtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgcg                                             <  2:167533/139‑1 (MQ=255)                                                    gtgCGTCTTTTCATGCAGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGATTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACAGCCCGGTGGTGTAATTTgcg                                             <  2:543636/139‑1 (MQ=255)                                                     tgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgcgc                                            >  2:69164/1‑139 (MQ=255)                                                           ttttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCg                                      <  1:423247/139‑1 (MQ=255)                                                            tttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                      >  2:67599/1‑74 (MQ=255)                                                            tttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                      <  1:67599/74‑1 (MQ=255)                                                                     gACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                            <  1:69164/139‑1 (MQ=255)                                                                      aCCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                           >  2:419979/1‑139 (MQ=255)                                                                                   ccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTgg              <  2:661070/139‑1 (MQ=255)                                                                                     aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCa                                                                                 >  2:791089/1‑70 (MQ=255)                                                                                     aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCa                                                                                 <  1:791089/70‑1 (MQ=255)                                                                                     aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa            >  2:207215/1‑139 (MQ=255)                                                                                       tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCaa                                                                      <  2:126738/79‑1 (MQ=255)                                                                                       tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCaa                                                                      >  1:126738/1‑79 (MQ=255)                                                                                       tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                            >  1:771128/1‑121 (MQ=255)                                                                                       tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                            <  2:771128/121‑1 (MQ=255)                                                                                       tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAt          >  1:494607/1‑139 (MQ=255)                                                                                                 aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                             >  2:402289/1‑62 (MQ=255)                                                                                                 aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                             <  1:402289/62‑1 (MQ=255)                                                                                                   gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCg  >  2:232467/1‑135 (MQ=255)                                                                                                   gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCg  <  1:232467/135‑1 (MQ=255)                                                                                                       cTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                            <  1:726916/73‑1 (MQ=255)                                                                                                       cTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                            >  2:726916/1‑73 (MQ=255)                                                                                                                |                                                                                                                          GGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCG  >  NC_000913/2226398‑2226630

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