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

CGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAA  >  NC_000913/2226375‑2226641                                                                                                                                       |                                                                                                                                     cGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTAtt                                                                                                                                  <  2:643759/139‑1 (MQ=255)          aTTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGa                                                                                                                         >  1:412438/1‑139 (MQ=255)           ttGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGtt                                                                                                                           >  2:641574/1‑136 (MQ=255)           ttGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGtt                                                                                                                           <  1:641574/136‑1 (MQ=255)                  ggATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                 <  1:1245271/139‑1 (MQ=255)                             tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACAt                                                                                                      >  2:934505/1‑139 (MQ=255)                             tCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGCGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATTTCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACAt                                                                                                      >  1:910381/1‑139 (MQ=255)                                           aCATCCGGAGTACCATGCCTGATACGTCGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        <  2:910381/139‑1 (MQ=255)                                               ccGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                          >  2:31557/1‑101 (MQ=255)                                               ccGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTg                                                                                                                          <  1:31557/101‑1 (MQ=255)                                               ccGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTAt                                                                                    >  1:541449/1‑139 (MQ=255)                                                       ccATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACt                                                                            >  1:981784/1‑139 (MQ=255)                                                          tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCCAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGACCTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                                         >  1:852502/1‑139 (MQ=255)                                                          tGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                                         <  1:659895/139‑1 (MQ=255)                                                             cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTggg                                                                                                                             <  1:1119416/84‑1 (MQ=255)                                                             cTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTggg                                                                                                                             >  2:1119416/1‑84 (MQ=255)                                                                aTACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggtgg                                                                   >  2:1129460/1‑139 (MQ=255)                                                                  acgacgCTTGTCCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGgtggtg                                                                 <  1:934505/139‑1 (MQ=255)                                                                    gacgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTa                                                               >  1:896695/1‑139 (MQ=255)                                                                      cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTTCCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        <  2:905538/112‑1 (MQ=255)                                                                      cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTggg                                                                                                                             >  2:391938/1‑75 (MQ=255)                                                                      cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTggg                                                                                                                             <  1:391938/75‑1 (MQ=255)                                                                      cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        >  1:905538/1‑112 (MQ=255)                                                                      cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGgtggtg                                                                 >  1:582793/1‑135 (MQ=255)                                                                      cgCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGgtggtg                                                                 <  2:582793/135‑1 (MQ=255)                                                                              cGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                           >  1:1251404/1‑85 (MQ=255)                                                                              cGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGcc                                                                                                           <  2:1251404/85‑1 (MQ=255)                                                                                tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                                         >  1:212305/1‑117 (MQ=255)                                                                                tCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGcc                                                                         <  2:212305/117‑1 (MQ=255)                                                                                    ttCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                               <  2:981784/139‑1 (MQ=255)                                                                                            gACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCca                                       <  2:412438/139‑1 (MQ=255)                                                                                                       tACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                 <  2:78894/54‑1 (MQ=255)                                                                                                       tACCCAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGa                                                                                                                 >  1:78894/1‑54 (MQ=255)                                                                                                         cccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTg                          >  2:152538/1‑139 (MQ=255)                                                                                                          ccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTgg                         <  2:735910/139‑1 (MQ=255)                                                                                                          ccAATGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTgg                         <  2:704740/139‑1 (MQ=255)                                                                                                              tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggt                                                                     <  2:1013821/91‑1 (MQ=255)                                                                                                              tGAAATATCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCggt                                                                     >  1:1013821/1‑91 (MQ=255)                                                                                                                    atCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                               <  1:972179/107‑1 (MQ=255)                                                                                                                    atCCACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                               >  2:972179/1‑107 (MQ=255)                                                                                                                      ccACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATc                                                                                          <  2:1358530/62‑1 (MQ=255)                                                                                                                      ccACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATc                                                                                          >  1:1358530/1‑62 (MQ=255)                                                                                                                       cACGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGt            <  2:541449/139‑1 (MQ=255)                                                                                                                        aCGCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGtt           >  1:161533/1‑139 (MQ=255)                                                                                                                        aCGCACCTCACAGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        <  1:427317/62‑1 (MQ=255)                                                                                                                        aCGCACCTCACAGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCa                                                                                        >  2:427317/1‑62 (MQ=255)                                                                                                                          gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                            <  2:526720/120‑1 (MQ=255)                                                                                                                          gCACCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                            >  1:526720/1‑120 (MQ=255)                                                                                                                            aCCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAt                     >  2:571505/1‑125 (MQ=255)                                                                                                                            aCCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAt                     <  1:571505/125‑1 (MQ=255)                                                                                                                            aCCTCACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTgcgc       >  1:771057/1‑139 (MQ=255)                                                                                                                                cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                               >  2:668912/1‑95 (MQ=255)                                                                                                                                cACTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGaa                                               <  1:668912/95‑1 (MQ=255)                                                                                                                                 aCTGCGTTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCaaaa  >  1:479218/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                     CGTAAAACCATTGCCCTGGATAAGGCGTTCGCGCCGCATCCGACATCCGGAGTACCATGCCTGATACGACGCTTGTGCGTCTTTTCATGCCGACCAAAATATTACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAA  >  NC_000913/2226375‑2226641

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