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%	213.9 / NA	66	*254* (TAG→TAA)	yohF	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/33);  total (33/33)

GCGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGC  >  NC_000913/2226374‑2226644                                                                                                                                        |                                                                                                                                        gCGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAt                                                                                                                                      >  1:1337514/1‑139 (MQ=255)  cGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                     <  2:1774444/139‑1 (MQ=255)      aaaCCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTgg                                                                                                                                 >  1:1906193/1‑139 (MQ=255)       aaCCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTggg                                                                                                                                <  2:101082/139‑1 (MQ=255)        aCCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGt                                                                                                                               <  1:1144873/139‑1 (MQ=255)          cATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                             <  2:1917586/139‑1 (MQ=255)                   ggATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                    <  1:968947/139‑1 (MQ=255)                         ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                     >  1:553253/1‑116 (MQ=255)                         ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                     <  2:553253/116‑1 (MQ=255)                         ggCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                              >  1:405966/1‑139 (MQ=255)                              tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACAt                                                                                                         >  1:1254992/1‑139 (MQ=255)                              tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACAt                                                                                                         >  1:346002/1‑139 (MQ=255)                               cgcgCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGt                                                                                                                               >  2:1992038/1‑116 (MQ=255)                               cgcgCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGt                                                                                                                               <  1:1992038/116‑1 (MQ=255)                                    cgcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCA‑CTCACTGCGTTATTCTggg                                                                                                                                <  1:1196787/109‑1 (MQ=255)                                    cgcATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCA‑CTCACTGCGTTATTCTggg                                                                                                                                >  2:1196787/1‑109 (MQ=255)                                          cGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATc                                                                                             <  2:405966/139‑1 (MQ=255)                                            aCATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                           <  2:346002/139‑1 (MQ=255)                                                    aGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAAc                                                                                   >  2:479108/1‑139 (MQ=255)                                                           tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGc                                                                                                               >  2:247893/1‑104 (MQ=255)                                                           tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGc                                                                                                               <  1:247893/104‑1 (MQ=255)                                                           tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGa                                                                                 <  1:1398577/134‑1 (MQ=255)                                                           tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCAACTATCAACGa                                                                                 >  2:1398577/1‑134 (MQ=255)                                                              cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                           <  2:1481822/121‑1 (MQ=255)                                                              cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                           >  1:1481822/1‑121 (MQ=255)                                                              cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCgg                                                                         >  2:52369/1‑139 (MQ=255)                                                              cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATC‑ACGACTGCCCg                                                                          >  1:161069/1‑137 (MQ=255)                                                              cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATC‑ACGACTGCCCg                                                                          <  2:161069/137‑1 (MQ=255)                                                                      acgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTaa                                                                 >  2:135947/1‑139 (MQ=255)                                                                      acgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTaa                                                                 >  2:1833109/1‑139 (MQ=255)                                                                       cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                                                     >  2:1385697/1‑102 (MQ=255)                                                                       cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAg                                                                                                     <  1:1385697/102‑1 (MQ=255)                                                                            gtgCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgcg                                                           <  1:52369/139‑1 (MQ=255)                                                                              gCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGcc                                                         <  1:479108/139‑1 (MQ=255)                                                                                gTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCt                                                       <  1:1833109/139‑1 (MQ=255)                                                                                 tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCt                                                                                                                                   >  1:1404398/1‑62 (MQ=255)                                                                                 tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCt                                                                                                                                   <  2:1404398/62‑1 (MQ=255)                                                                                     ttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACt                                                                                                                         <  2:1058260/68‑1 (MQ=255)                                                                                     ttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACt                                                                                                                         >  1:1058260/1‑68 (MQ=255)                                                                                                         aCCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCt                              >  1:1278313/1‑139 (MQ=255)                                                                                                          cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                                  <  2:1873375/118‑1 (MQ=255)                                                                                                          cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                                  >  1:1873375/1‑118 (MQ=255)                                                                                                          cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                          >  2:822747/1‑126 (MQ=255)                                                                                                          cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                          <  1:822747/126‑1 (MQ=255)                                                                                                             aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                                  >  2:284990/1‑67 (MQ=255)                                                                                                             aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAgccg                                                                                                  <  1:284990/67‑1 (MQ=255)                                                                                                             aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                                 >  1:1701499/1‑116 (MQ=255)                                                                                                             aaTGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                                 <  2:1701499/116‑1 (MQ=255)                                                                                                               tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                          <  1:1668727/89‑1 (MQ=255)                                                                                                               tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCg                                                                          >  2:1668727/1‑89 (MQ=255)                                                                                                               tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAt                        >  1:997663/1‑139 (MQ=255)                                                                                                                        cACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGt               <  2:997663/139‑1 (MQ=255)                                                                                                                         aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                    <  2:1162257/37‑1 (MQ=255)                                                                                                                         aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                    >  1:1162257/1‑37 (MQ=255)                                                                                                                         aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgc                                                            <  2:1521224/93‑1 (MQ=255)                                                                                                                         aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTgc                                                            >  1:1521224/1‑93 (MQ=255)                                                                                                                         aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                               >  2:1050159/1‑122 (MQ=255)                                                                                                                         aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                               <  1:1050159/122‑1 (MQ=255)                                                                                                                           gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                           <  2:286925/60‑1 (MQ=255)                                                                                                                           gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                           >  1:286925/1‑60 (MQ=255)                                                                                                                           gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGgtggtg                                                                    >  2:233014/1‑83 (MQ=255)                                                                                                                           gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGgtggtg                                                                    <  1:233014/83‑1 (MQ=255)                                                                                                                                 cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAAc                                                                                   >  2:257887/1‑62 (MQ=255)                                                                                                                                 cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAAc                                                                                   <  1:257887/62‑1 (MQ=255)                                                                                                                                    tGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                          >  2:1597463/1‑100 (MQ=255)                                                                                                                                    tGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                          <  1:1597463/100‑1 (MQ=255)                                                                                                                                     gcgtTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAAcgc  <  2:1254992/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        GCGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGC  >  NC_000913/2226374‑2226644

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