Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,233,669	G→A	A669T (GCG→ACG)	fadH →	2,4‑dienoyl‑CoA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,233,669	0	G	A	100.0%	109.3 / NA	32	A669T (GCG→ACG)	fadH	2,4‑dienoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/14);  total (18/14)

TCTGCGCAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTC  >  NC_000913/3233532‑3233803                                                                                                                                          |                                                                                                                                       tCTGCGCAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGacg                                                                                                                                      >  1:919096/1‑140 (MQ=255)   tGCGCAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCt                                                                                                                                    >  2:803182/1‑140 (MQ=255)   tGCGCAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCt                                                                                                                                    >  2:501739/1‑140 (MQ=255)    gcgcAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTg                                                                                                                                   >  2:970914/1‑140 (MQ=255)    gcgcAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTg                                                                                                                                   <  1:991227/140‑1 (MQ=255)     cgcAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGTAGCTGGACCCACGACGGGCAATTTCCCCGGGAACACGGCTGACGCTgg                                                                                                                                  >  2:525953/1‑140 (MQ=255)     cgcAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTgg                                                                                                                                  <  2:919096/140‑1 (MQ=255)       cAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGgag                                                                                                                                <  2:888215/140‑1 (MQ=255)              agagCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTaaa                                                                                                                         <  1:310389/140‑1 (MQ=255)              agagCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTaaa                                                                                                                         <  2:216264/140‑1 (MQ=255)              agagCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTaaa                                                                                                                         <  2:1022460/140‑1 (MQ=255)               gagCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAAt                                                                                                                        >  1:543201/1‑140 (MQ=255)               gagCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAAt                                                                                                                        >  1:58909/1‑140 (MQ=255)                            cgcTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAgc                                                                                                           >  2:177533/1‑140 (MQ=255)                                 gcgcAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCgg                                                                                                      >  2:1083874/1‑140 (MQ=255)                                        cGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAAc                                                                                               <  2:800896/140‑1 (MQ=255)                                                 atAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGc                                                                                      >  1:669003/1‑140 (MQ=255)                                                       ggAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTcg                                                                                >  2:478905/1‑140 (MQ=255)                                                       ggAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTcg                                                                                <  2:58909/140‑1 (MQ=255)                                                           aaCGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGtt                                                                            <  1:478905/140‑1 (MQ=255)                                                            aCGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTa                                                                           <  2:669003/140‑1 (MQ=255)                                                                  cATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAAt                                                                     >  2:894658/1‑140 (MQ=255)                                                                              ggcTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACaa                                                               <  2:81175/134‑1 (MQ=255)                                                                              ggcTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACaa                                                               >  1:81175/1‑134 (MQ=255)                                                                                               ggAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcgcg                                        >  2:1036342/1‑140 (MQ=255)                                                                                                  gCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTg                                     >  1:147779/1‑140 (MQ=255)                                                                                                     ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCa                                  >  1:147755/1‑140 (MQ=255)                                                                                                           acgacgGTCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttttt                            >  2:645903/1‑140 (MQ=255)                                                                                                           acgacgGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttttt                            >  2:781518/1‑140 (MQ=255)                                                                                                             gacgGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTa                          <  1:781518/140‑1 (MQ=255)                                                                                                               cgGGCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                        <  1:645903/140‑1 (MQ=255)                                                                                                                gggCAATTGCCCAGGGAACACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCg                       <  2:147779/140‑1 (MQ=255)                                                                                                                                aaCACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAAtt       >  2:316515/1‑140 (MQ=255)                                                                                                                                aaCACGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAAtt       >  1:126729/1‑140 (MQ=255)                                                                                                                                   acGGCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc    >  1:922764/1‑140 (MQ=255)                                                                                                                                     ggCTGACGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCctc  >  2:396561/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       TCTGCGCAGGGCAAGAGCCAAACCGCGCGCTGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTC  >  NC_000913/3233532‑3233803

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