Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,236,414	A→C	D152A (GAC→GCC)	ygjP →	putative metal dependent hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,236,414	0	A	C	100.0%	159.1 / NA	42	D152A (GAC→GCC)	ygjP	putative metal dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (18/24);  total (18/24)

ACCGTGTTTCGCAACGCGCCGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGACACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGC  >  NC_000913/3236272‑3236550                                                                                                                                               |                                                                                                                                         aCCGTGTTTCGCAACGCGCCGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGa                                                                                                                                   <  1:49604/150‑1 (MQ=255)      gtTTCGCAACGCGCCGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTgg                                                                                                                             >  2:305909/1‑151 (MQ=255)          cGCAACGGGCGGGAACCGTTGTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAAAGCACAACAAAGCGTTTTACCAGTTGTGTTGTCATAGGGAACCCCAGTGCCGCGGGCTTGAGTGCGCCACCCGACTGTGGCt                                                                                                                           <  2:580343/149‑1 (MQ=255)           gCAACGCGCCGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAAc                                                                                                                        <  1:580558/151‑1 (MQ=255)                 cgcCGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGtt                                                                                                                  <  1:660228/151‑1 (MQ=255)                    cGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTc                                                                                                               <  1:225508/151‑1 (MQ=255)                       aaCCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCt                                                                                                            >  2:681312/1‑151 (MQ=255)                                   gcATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACaa                                                                                                <  2:379428/151‑1 (MQ=255)                                     aTGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACaaaa                                                                                              <  2:478274/151‑1 (MQ=255)                                         tCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCGCCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCt                                                                                           <  2:125623/150‑1 (MQ=255)                                            tggtgCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGaa                                                                                        >  1:165051/1‑150 (MQ=255)                                              gtgCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAat                                                                                     <  2:21447/151‑1 (MQ=255)                                                    gAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTc                                                                               <  1:421330/151‑1 (MQ=255)                                                      gCTGGCGCACCTGAAGGAGAAAGAGCATAACAACGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTc                                                                               >  2:398480/1‑149 (MQ=255)                                                          gcgcACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGt                                                                         <  1:253277/151‑1 (MQ=255)                                                            gcACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTa                                                                       <  1:362520/151‑1 (MQ=255)                                                               ccTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTAtct                                                                    >  1:584899/1‑151 (MQ=255)                                                                   aaGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTAtcttttc                                                                <  1:38007/151‑1 (MQ=255)                                                                     ggAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTAtcttttctt                                                              <  1:586121/151‑1 (MQ=255)                                                                     ggAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTAtcttttctt                                                              >  2:459626/1‑151 (MQ=255)                                                                              gCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAgcg                                                     >  1:191708/1‑151 (MQ=255)                                                                               cATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAgcgc                                                    <  2:552616/151‑1 (MQ=255)                                                                                      aaGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAAt                                             <  1:452814/151‑1 (MQ=255)                                                                                        gCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCg                                           >  2:647589/1‑151 (MQ=255)                                                                                         cGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCg                                           >  2:673559/1‑150 (MQ=255)                                                                                         cGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGa                                          <  1:327851/151‑1 (MQ=255)                                                                                          gTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGAc                                         <  2:575121/151‑1 (MQ=255)                                                                                           ttttACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGAc                                         <  2:674614/150‑1 (MQ=255)                                                                                              tACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTcc                                     >  2:470148/1‑151 (MQ=255)                                                                                                ccAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTcccc                                   >  1:646645/1‑151 (MQ=255)                                                                                                  aGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTcccccc                                 >  1:144646/1‑151 (MQ=255)                                                                                                  aGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTcccccc                                 >  1:141029/1‑151 (MQ=255)                                                                                                   gttgtgttgtCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCt                                <  1:103480/151‑1 (MQ=255)                                                                                                          tgtCATATGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGtt                          >  2:498806/1‑150 (MQ=255)                                                                                                                atGGAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGAt                    >  1:449560/1‑150 (MQ=255)                                                                                                                   gAACCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATGGACCTCCCCCCTTTCGTTAAGGATAAc                 >  2:695773/1‑150 (MQ=255)                                                                                                                    aaCCCCAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACg                >  1:507698/1‑150 (MQ=255)                                                                                                                        ccAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCAt           <  1:398622/151‑1 (MQ=255)                                                                                                                        ccAGTACCACCAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCAt           <  1:114552/151‑1 (MQ=255)                                                                                                                             accaccAGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCAtgctgc      <  1:403841/151‑1 (MQ=255)                                                                                                                                 ccaGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTgcgcgc  >  2:512432/1‑151 (MQ=255)                                                                                                                                 ccaGCTTGAGTTCGCCACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTgcgcg   <  2:602753/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                         ACCGTGTTTCGCAACGCGCCGGAACCGTTTTTGCGCATGATCGTGGTGCACGAGCTGGCGCACCTGAAGGAGAAAGAGCATAACAAAGCGTTTTACCAGTTGTGTTGTCATATGGAACCCCAGTACCACCAGCTTGAGTTCGACACCCGACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGC  >  NC_000913/3236272‑3236550

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