Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,031,019	G→A	intergenic (‑48/‑142)	fadB ← / → pepQ	fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase/proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,031,019	0	G	A	100.0%	174.5 / NA	46	intergenic (‑48/‑142)	fadB/pepQ	fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase/proline dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/27);  total (19/27)

AGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGA  >  NC_000913/4030873‑4031165                                                                                                                                                   |                                                                                                                                                   aGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAg                                                                                                                                                 >  1:627253/1‑150 (MQ=255)  gtgtCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAgtg                                                                                                                                               <  2:178824/151‑1 (MQ=255)    gtCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAgtgtg                                                                                                                                             >  2:468332/1‑151 (MQ=255)      cGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTa                                                                                                                                           <  2:224677/151‑1 (MQ=255)         gTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTAttt                                                                                                                                        <  2:112267/151‑1 (MQ=255)            tattaACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTtgt                                                                                                                                     >  2:444729/1‑151 (MQ=255)            tattaACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTtgt                                                                                                                                     >  1:473402/1‑151 (MQ=255)            tattaACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTtg                                                                                                                                      <  2:254611/150‑1 (MQ=255)             attaACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTtgtg                                                                                                                                    <  2:29775/151‑1 (MQ=255)              ttaACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTgtgt                                                                                                                                   >  2:410316/1‑151 (MQ=255)              ttaACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTgtgt                                                                                                                                   >  1:325765/1‑151 (MQ=255)                    gAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTaaaaa                                                                                                                             <  1:549271/151‑1 (MQ=255)                    gAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTaaaa                                                                                                                              >  2:182653/1‑150 (MQ=255)                    gAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTaaa                                                                                                                               <  2:148912/149‑1 (MQ=255)                       ccTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAatat                                                                                                                          <  2:117128/151‑1 (MQ=255)                          ggggCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGca                                                                                                                       <  1:26231/151‑1 (MQ=255)                                       ccAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTa                                                                                                           <  2:225834/150‑1 (MQ=255)                                         aGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATCTTGTGTTAAAAATATGCAAACAAAAGATTAAAg                                                                                                        <  2:362742/151‑1 (MQ=255)                                            tCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGaa                                                                                                      <  1:50347/150‑1 (MQ=255)                                            tCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGaa                                                                                                      <  1:272446/150‑1 (MQ=255)                                               gCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGc                                                                                                  <  1:67263/151‑1 (MQ=255)                                                  aTGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGa                                                                                               <  2:318636/151‑1 (MQ=255)                                                      cATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCGTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTCGTGTTAAAAAGATGCAAACAAAAGATAAAAGAAATGGCGATCt                                                                                            <  2:506911/150‑1 (MQ=255)                                                           tCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGa                                                                                          <  2:503555/147‑1 (MQ=255)                                                            ccagccagTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACg                                                                                     >  1:28353/1‑151 (MQ=255)                                                                           aCAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCa                                                                       <  1:638487/150‑1 (MQ=255)                                                                             aGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGc                                                                    <  1:653039/151‑1 (MQ=255)                                                                                    cGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGAt                                                             <  1:67495/151‑1 (MQ=255)                                                                                      ccTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGa                                                           <  1:519646/151‑1 (MQ=255)                                                                                           gTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGAt                                                      >  2:462014/1‑151 (MQ=255)                                                                                             aaaGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGAtg                                                     >  1:130850/1‑150 (MQ=255)                                                                                                     tCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTa                                                >  2:50250/1‑147 (MQ=255)                                                                                                      cAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATg                                           <  1:704728/151‑1 (MQ=255)                                                                                                           tcCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCgg                                        >  2:530724/1‑149 (MQ=255)                                                                                                                      cAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTc                            >  1:132589/1‑150 (MQ=255)                                                                                                                      cAAGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCa                           <  2:59060/151‑1 (MQ=255)                                                                                                                       aaGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCa                           >  2:648591/1‑150 (MQ=255)                                                                                                                        aGGTGATCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCaaa                         <  1:719578/151‑1 (MQ=255)                                                                                                                              tCTGGTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAAc                    >  1:440274/1‑150 (MQ=255)                                                                                                                                  gTCGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAg               <  1:310849/151‑1 (MQ=255)                                                                                                                                    cGTACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAggg             >  2:344077/1‑151 (MQ=255)                                                                                                                                      tACCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGt            >  2:379894/1‑150 (MQ=255)                                                                                                                                       aCCAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTa           >  1:102734/1‑150 (MQ=255)                                                                                                                                        ccAGATGAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTaa          >  2:491589/1‑150 (MQ=255)                                                                                                                                              gAGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATg    <  2:714357/150‑1 (MQ=255)                                                                                                                                               aGTCAAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGa  <  2:197339/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGTGTCGAGTTTATTAACTGAACCTGGGGCATCAAATACCAGTTCGGCAATGCCATCTTCCAGCCAGTCAAGGTACAGGGTGTCGCCTTTGTAAAGCATGTCAGTCTCCTGAATCCGCAAGGTGATCTGGTCGTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAATGGA  >  NC_000913/4030873‑4031165

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