Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,319,158	T→C	Y150C (TAC→TGC)	phnJ ←	carbon‑phosphorus lyase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,319,158	0	T	C	100.0%	192.0 / NA	53	Y150C (TAC→TGC)	phnJ	carbon‑phosphorus lyase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/31);  total (22/31)

GTCCATTTTTGGGTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGTACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTC  >  NC_000913/4319016‑4319304                                                                                                                                               |                                                                                                                                                   gTCCATTTTTGGGTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTcc                                                                                                                                            <  2:687078/151‑1 (MQ=255)       ttttGGGTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGcgcg                                                                                                                                      >  2:428437/1‑151 (MQ=255)         ttGGGTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGt                                                                                                                                     <  1:110779/150‑1 (MQ=255)         ttGGGTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGt                                                                                                                                     <  1:355794/150‑1 (MQ=255)           gggTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGc                                                                                                                                   <  1:239426/150‑1 (MQ=255)                       ttCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGtt                                                                                                                       <  2:194337/150‑1 (MQ=255)                       ttCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGtt                                                                                                                       <  1:74132/150‑1 (MQ=255)                            gATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGtt                                                                                                                 >  1:707400/1‑151 (MQ=255)                            gATCGGCGACGGGTCCATCACGTAGCGCCCATTTAACTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCAAGCGCGTGCATGGTGCGGGTTTCCGtt                                                                                                                 <  2:318736/151‑1 (MQ=255)                                 gCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTcgcg                                                                                                            >  2:720312/1‑151 (MQ=255)                                 gCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTcgcg                                                                                                            <  2:739781/151‑1 (MQ=255)                                  cGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTcgcgc                                                                                                           >  2:894778/1‑151 (MQ=255)                                   gACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTcgcgcg                                                                                                          >  1:904491/1‑151 (MQ=255)                                       ggTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTc                                                                                                      <  1:775199/151‑1 (MQ=255)                                          ccATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGa                                                                                                    <  2:825515/150‑1 (MQ=255)                                              cACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCCTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATaaa                                                                                                <  2:683548/150‑1 (MQ=255)                                                  tAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAgcgc                                                                                            >  1:835612/1‑150 (MQ=255)                                                    gcgcCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGc                                                                                         <  1:797327/151‑1 (MQ=255)                                                     cgcCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAg                                                                                          >  1:365655/1‑149 (MQ=255)                                                          aTTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCg                                                                                   >  2:415109/1‑151 (MQ=255)                                                                  tCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGAtt                                                                             <  1:287973/149‑1 (MQ=255)                                                                      cGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCAcc                                                                       >  2:16947/1‑151 (MQ=255)                                                                       ggATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCt                                                                      <  2:350512/151‑1 (MQ=255)                                                                        gATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTg                                                                     <  2:903155/151‑1 (MQ=255)                                                                           aGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGaa                                                                  <  2:314393/151‑1 (MQ=255)                                                                           aGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGa                                                                   <  1:683647/150‑1 (MQ=255)                                                                                tAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTAACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATa                                                              >  1:523204/1‑150 (MQ=255)                                                                                     ggtggCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATc                                                        <  1:724751/151‑1 (MQ=255)                                                                                        ggCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTg                                                      <  2:667699/150‑1 (MQ=255)                                                                                         gCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGAt                                                    <  1:304719/151‑1 (MQ=255)                                                                                         gCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGAt                                                    <  1:61970/151‑1 (MQ=255)                                                                                           gATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCt                                                  >  1:890913/1‑151 (MQ=255)                                                                                           gATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCt                                                  <  2:365629/151‑1 (MQ=255)                                                                                           gATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCt                                                  <  1:415154/151‑1 (MQ=255)                                                                                           gATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCt                                                  >  1:147501/1‑151 (MQ=255)                                                                                           gATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCt                                                  <  2:43977/151‑1 (MQ=255)                                                                                              atGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTc                                                <  1:753420/150‑1 (MQ=255)                                                                                                gACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCgg                                              <  1:16969/150‑1 (MQ=255)                                                                                                    gAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAgc                                         <  2:516806/151‑1 (MQ=255)                                                                                                        cGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAgcggcg                                     >  2:121299/1‑151 (MQ=255)                                                                                                            cGATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGttt                                  <  1:775966/150‑1 (MQ=255)                                                                                                             gATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCg                                >  1:126721/1‑151 (MQ=255)                                                                                                             gATATCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCg                                >  1:907061/1‑151 (MQ=255)                                                                                                              atatCTTCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCgg                               <  1:304081/151‑1 (MQ=255)                                                                                                                    tCATACAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGc                         >  1:178376/1‑151 (MQ=255)                                                                                                                        aCAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTg                     >  1:406159/1‑151 (MQ=255)                                                                                                                         cAGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTg                     >  1:34880/1‑150 (MQ=255)                                                                                                                          aGTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGAc                   >  2:655557/1‑151 (MQ=255)                                                                                                                           gTTTCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGAcg                  <  1:78927/151‑1 (MQ=255)                                                                                                                              tCACCTGCATCACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGt               >  2:387066/1‑151 (MQ=255)                                                                                                                                        cACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAAcc     >  2:906466/1‑151 (MQ=255)                                                                                                                                        cACGCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAAc      <  2:79902/150‑1 (MQ=255)                                                                                                                                           gCCGCACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTc  >  2:529408/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                   GTCCATTTTTGGGTTATCGAATTTCGGGATCGGCGACGGGTCCATCACGTAGCGCCCATTTACCTTCACCGGATAGGCGTAGGTGGTGGCGATATGACCGAAGCGGGCGATATCTTCATACAGTTTCACCTGCATCACGCCGTACTCTTCCAGCGCGTGCATGGTGCGGGTTTCCGTTTCGCGCGGCTCGATAAAGCGCAGCGGCTCGGGGATTGGCACCTGGAAGATAATTATCTGATCTTCGGTCAGCGGCGTTTCGGGGATGCGGTGACGCGTCTGGATAACCGTC  >  NC_000913/4319016‑4319304

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