Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,086,811	Δ1 bp	coding (796/834 nt)	fdhD →	sulfurtransferase for molybdenum cofactor sulfuration

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,086,809	0	G	.	100.0%	172.6 / NA	41	coding (794/834 nt)	fdhD	sulfurtransferase for molybdenum cofactor sulfuration
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (23/18);  total (23/18)

GTGCCAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGG  >  NC_000913/4086671‑4086941                                                                                                                                           |                                                                                                                                     gTGCCAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCgg                                                                                                                                     >  2:16512/1‑140 (MQ=255)    ccAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCggcaa                                                                                                                                  <  1:509510/140‑4 (MQ=255)    ccAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCggcaa                                                                                                                                  <  2:468490/140‑4 (MQ=255)    ccAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCggcaa                                                                                                                                  >  1:589900/1‑137 (MQ=255)    ccAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCggcaa                                                                                                                                  <  1:16512/140‑4 (MQ=255)    ccAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCggcaa                                                                                                                                  <  1:464368/140‑4 (MQ=255)    ccAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCggcaa                                                                                                                                  <  1:266067/140‑4 (MQ=255)      aGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACg                                                                                                                               >  2:348092/1‑140 (MQ=255)           tGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTta                                                                                                                          <  1:348092/140‑1 (MQ=255)             aaaTGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTtata                                                                                                                        <  2:72108/140‑1 (MQ=255)                 ggTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCa                                                                                                                    >  1:323364/1‑140 (MQ=255)                  gTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCAt                                                                                                                   >  2:942816/1‑140 (MQ=255)                  gTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCAt                                                                                                                   >  2:841458/1‑140 (MQ=255)                  gTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCAt                                                                                                                   >  1:6049/1‑140 (MQ=255)                     cAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCg                                                                                                                <  2:900801/140‑1 (MQ=255)                     cAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCg                                                                                                                >  2:484689/1‑140 (MQ=255)                     cAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCg                                                                                                                <  2:323364/140‑1 (MQ=255)                     cAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCg                                                                                                                >  1:900801/1‑140 (MQ=255)                      aaaaGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGc                                                                                                               <  1:942816/140‑1 (MQ=255)                        aaGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAg                                                                                                             >  1:683420/1‑140 (MQ=255)                             ggCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGTAGCGttt                                                                                                        >  1:242518/1‑140 (MQ=255)                                   gtgCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCaa                                                                                                  <  2:589900/140‑1 (MQ=255)                                   gtgCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCaa                                                                                                  <  2:834309/140‑1 (MQ=255)                                         cgTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGTAGCGTTTGAGCaattaatt                                                                                            <  2:242518/140‑1 (MQ=255)                                            agagATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTaaa                                                                                         >  1:712965/1‑140 (MQ=255)                                            agagATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTaaa                                                                                         >  2:532530/1‑140 (MQ=255)                                            agagATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTaaa                                                                                         >  2:128559/1‑140 (MQ=255)                                                  tttgtttgCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGt                                                                                   >  1:761868/1‑140 (MQ=255)                                                     gtttgCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCt                                                                                <  2:761868/140‑1 (MQ=255)                                                        tgCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTGAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCg                                                                             <  1:529782/140‑1 (MQ=255)                                                          cGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGaa                                                                           <  2:932828/140‑1 (MQ=255)                                                                      cgACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCt                                                                                <  1:757447/123‑1 (MQ=255)                                                                      cgACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCt                                                                                >  2:757447/1‑123 (MQ=255)                                                                        aCCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGc                                                             <  2:683420/140‑1 (MQ=255)                                                                        aCCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGc                                                             <  1:484689/140‑1 (MQ=255)                                                                        aCCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGc                                                             <  2:712965/140‑1 (MQ=255)                                                                         ccACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCa                                                            <  2:6049/140‑1 (MQ=255)                                                                         ccACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCa                                                            <  1:4923/140‑1 (MQ=255)                                                                           aCGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAgg                                                          >  2:942771/1‑140 (MQ=255)                                                                             gCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGAc                                                        >  2:363863/1‑140 (MQ=255)                                                                                            ggCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGa                                         >  1:530583/1‑140 (MQ=255)                                                                                             gCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GTCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGaa                                        >  2:429766/1‑140 (MQ=255)                                                                                                      tGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAattttattat                               >  2:98220/1‑140 (MQ=255)                                                                                                           tCTGACACTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAattttattattttat                          >  1:456846/1‑140 (MQ=255)                                                                                                                cacTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTtata                     <  2:530583/140‑1 (MQ=255)                                                                                                                  cTGGTAGGTTTTTGTAAACCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATaaa                   >  2:402211/1‑140 (MQ=255)                                                                                                                                   aCCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTgg  >  1:605431/1‑140 (MQ=255)                                                                                                                                   aCCGGGTC‑GGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTgg  >  1:403678/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                     GTGCCAGTTATGAAATGGTGCAAAAGTCGGCGATGTGCGGCGTAGAGATTTTGTTTGCGGTGTCTGCCGCGACCACGCTTGCTGTAGAAGTGGCTGAGCGCTGTAATCTGACACTGGTAGGTTTTTGTAAACCGGGTCGGGCAACGGTTTATACCCATCCGCAGCGTTTGAGCAATTAATTAAAATAAGTCCTGCGAAATATATTAAATTGCAGGACTTATTCATTTCGTGAATTTTATTATTTTATTTATAAACAAAATGCGTGTAGTGG  >  NC_000913/4086671‑4086941

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