Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,031,161	(A)6→7	coding (1/1332 nt)	pepQ →	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,031,155	1	.	A	100.0%	194.5 / NA	48	intergenic (‑184/‑6)	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 . (0/0);  new base A (27/21);  total (27/21)

GTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGT‑AAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCC  >  NC_000913/4031005‑4031300                                                                                                                                                        |                                                                                                                                                  gTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGt                                                                                                                                                    <  1:463148/151‑1 (MQ=255)       gATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGT‑aaaaaa                                                                                                                                             <  2:324140/151‑1 (MQ=255)       gATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGT‑aaaaaa                                                                                                                                             <  2:244200/151‑1 (MQ=255)                  gtgtgtATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACt                                                                                                                                   >  2:426151/1‑151 (MQ=255)                  gtgtgtATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACt                                                                                                                                   >  2:505264/1‑151 (MQ=255)                        aTTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCt                                                                                                                              <  2:280966/150‑1 (MQ=255)                        aTTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTc                                                                                                                             >  2:467036/1‑151 (MQ=255)                            tgtgTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGctc                                                                                                                         >  2:517216/1‑151 (MQ=255)                            tgtgTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGctc                                                                                                                         >  2:536864/1‑151 (MQ=255)                               gtTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCta                                                                                                                       >  1:471282/1‑150 (MQ=255)                                 tAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCtat                                                                                                                      >  1:188788/1‑149 (MQ=255)                                  aaaaaTATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATaaa                                                                                                                   <  1:591726/151‑1 (MQ=255)                                  aaaaaTATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATaa                                                                                                                    <  2:473105/150‑1 (MQ=255)                                  aaaaaTATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATaa                                                                                                                    >  1:205110/1‑150 (MQ=255)                                    aaaTATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATaaaaa                                                                                                                 <  1:145291/151‑1 (MQ=255)                                    aaaTATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATaaaa                                                                                                                  >  2:464157/1‑150 (MQ=255)                                        atGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCa                                                                                                              <  1:567989/150‑1 (MQ=255)                                            aaacaaaAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCatat                                                                                                           >  2:589843/1‑149 (MQ=255)                                            aaacaaaAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAg                                                                                                         <  1:464163/151‑1 (MQ=255)                                                aaaaGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCGCGCTCTATAAAAATCATATAGCTa                                                                                                      <  1:26719/150‑1 (MQ=255)                                                       tAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACa                                                                                               <  2:492025/150‑1 (MQ=255)                                                       tAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACa                                                                                               <  1:51071/150‑1 (MQ=255)                                                               tGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGAc                                                                                      >  1:47297/1‑151 (MQ=255)                                                               tGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGAc                                                                                      >  1:280672/1‑151 (MQ=255)                                                                 ccGCTCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTc                                                                                    >  1:499909/1‑151 (MQ=255)                                                                  cGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTcg                                                                                   >  2:164120/1‑151 (MQ=255)                                                                   gATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTc                                                                                    >  1:378941/1‑149 (MQ=255)                                                                      cTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGAt                                                                               >  2:639216/1‑151 (MQ=255)                                                                        gATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATgc                                                                             >  1:204840/1‑151 (MQ=255)                                                                           cACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATgcgc                                                                           <  1:490671/150‑1 (MQ=255)                                                                            aCGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTg                                                                         <  2:94292/151‑1 (MQ=255)                                                                             cGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTgg                                                                        >  2:564282/1‑151 (MQ=255)                                                                                tCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGgcgc                                                                     >  2:26200/1‑151 (MQ=255)                                                                                    cAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCtt                                                                 >  1:3262/1‑151 (MQ=255)                                                                                              gCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGAt                                                       >  1:44634/1‑151 (MQ=255)                                                                                                tctGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGc                                                     >  1:439617/1‑151 (MQ=255)                                                                                                 ctGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATg                                                      <  2:223254/149‑1 (MQ=255)                                                                                                   ggATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGtt                                                  <  1:257080/151‑1 (MQ=255)                                                                                                       gAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACtt                                              >  2:9482/1‑151 (MQ=255)                                                                                                       gAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACt                                               <  1:427258/150‑1 (MQ=255)                                                                                                             tgtgCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCAc                                        >  1:146723/1‑151 (MQ=255)                                                                                                                     gATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGa                                <  2:683941/151‑1 (MQ=255)                                                                                                                        gCGGAGACTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCt                             <  1:235855/151‑1 (MQ=255)                                                                                                                         cgtagaCTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCATCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTg                            >  1:339861/4‑151 (MQ=255)                                                                                                                               cTTAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAAt                      >  2:440122/1‑151 (MQ=255)                                                                                                                                 tAAGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATg                     <  2:227198/150‑1 (MQ=255)                                                                                                                                  aaGGTCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGtt                   <  2:47297/151‑1 (MQ=255)                                                                                                                                      tCAAAAAAACAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTtc               >  2:519005/1‑151 (MQ=255)                                                                                                                                         aaaaaaCAGAAGGTTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGTATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTtctc             <  1:173113/150‑1 (MQ=255)                                                                                                                                         aaaaaaCAGAAGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTcg            <  1:492157/151‑1 (MQ=255)                                                                                                                                                   aGGGTAAAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATcc  <  1:539123/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  GTACCAGATGAGTCGAAGTGTGTATTTTGTGTTAAAAATATGCAAACAAAAGATTAAAGAAATGCCGATCTGATCACGCTCGGCAGAAATCACGCTCTGGATGAACGATGTGCTAAGATGCGGAGACTTAAGGTCAAAAAAACAGAAGGGT‑AAAAAATGGAATCACTGGCCTCGCTCTATAAAAATCATATAGCTACCTTACAAGAACGGACTCGCGATGCGCTGGCGCGCTTCAAGCTGGATGCGTTACTTATTCACTCCGGCGAGCTGTTCAATGTTTTTCTCGACGATCATCC  >  NC_000913/4031005‑4031300

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