Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,276,646	T→G	F392L (TTT→TTG)	yjcC →	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,276,646	0	T	G	100.0%	186.9 / NA	48	F392L (TTT→TTG)	yjcC	putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/33);  total (15/33)

CGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTC  >  NC_000913/4276500‑4276787                                                                                                                                                   |                                                                                                                                              cgcgATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAg                                                                                                                                           <  1:496015/151‑1 (MQ=255)            cGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAgcat                                                                                                                                <  2:415193/151‑1 (MQ=255)            cGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAgcat                                                                                                                                <  1:434067/151‑1 (MQ=255)             gATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAgcatg                                                                                                                               >  2:893940/1‑151 (MQ=255)              aTTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAgcatg                                                                                                                               <  1:506919/150‑1 (MQ=255)                tACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAgcatgcat                                                                                                                            >  1:387788/1‑151 (MQ=255)                           acaTGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGtt                                                                                                                 >  2:176674/1‑151 (MQ=255)                                cAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACa                                                                                                             >  1:869046/1‑150 (MQ=255)                                  gATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACaaaa                                                                                                          <  2:557676/151‑1 (MQ=255)                                     cGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATg                                                                                                        <  1:394427/150‑1 (MQ=255)                                     cGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATg                                                                                                        >  1:97648/1‑150 (MQ=255)                                     cGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGa                                                                                                       <  2:275265/151‑1 (MQ=255)                                        tATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGc                                                                                                    >  1:705487/1‑151 (MQ=255)                                         aTGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACg                                                                                                     <  1:663561/149‑1 (MQ=255)                                         aTGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGcc                                                                                                   >  1:893740/1‑151 (MQ=255)                                           gTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGa                                                                                                 <  2:88115/151‑1 (MQ=255)                                              tCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGatt                                                                                               <  2:618547/150‑1 (MQ=255)                                                 aTTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGatta                                                                                              <  2:393980/148‑1 (MQ=255)                                                      ccTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCt                                                                                      <  1:541710/151‑1 (MQ=255)                                                         gTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCtt                                                                                     <  2:817663/149‑1 (MQ=255)                                                         gTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTc                                                                                   <  1:180010/151‑1 (MQ=255)                                                              ccTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGcc                                                                               >  1:708373/1‑150 (MQ=255)                                                                   gATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCcaggc                                                                           >  2:69235/1‑149 (MQ=255)                                                                        tCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTa                                                                     <  2:231058/150‑1 (MQ=255)                                                                              gTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTg                                                              <  2:680356/151‑1 (MQ=255)                                                                                 aCGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGc                                                            >  2:370541/1‑150 (MQ=255)                                                                                  cGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCa                                                           >  2:717122/1‑150 (MQ=255)                                                                                       gATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATtgat                                                     <  2:444205/151‑1 (MQ=255)                                                                                        aTAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATtgatg                                                    >  2:592253/1‑151 (MQ=255)                                                                                             gcgAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGAtttt                                               <  2:387699/151‑1 (MQ=255)                                                                                             gcgAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGAtttt                                               <  1:622328/151‑1 (MQ=255)                                                                                               gAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTgg                                             <  1:474471/151‑1 (MQ=255)                                                                                                 atcaatcaGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGt                                            >  1:298925/1‑150 (MQ=255)                                                                                                     atcaGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTgg                                       <  1:808755/151‑1 (MQ=255)                                                                                                           aaaCAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTActc                                 >  2:122802/1‑151 (MQ=255)                                                                                                                       aCGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATaa                     <  2:79705/151‑1 (MQ=255)                                                                                                                       aCGCGGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATa                      <  2:286786/150‑1 (MQ=255)                                                                                                                        cgcgGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAAc                    <  2:893720/151‑1 (MQ=255)                                                                                                                        cgcgGTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAAc                    >  2:522504/1‑151 (MQ=255)                                                                                                                           ggTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCtt                 <  1:140418/151‑1 (MQ=255)                                                                                                                           ggTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCtt                 <  2:690449/151‑1 (MQ=255)                                                                                                                            gTGCGTCCGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTa                <  2:10172/151‑1 (MQ=255)                                                                                                                                  ccGCAGCAAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCAt          <  1:256858/151‑1 (MQ=255)                                                                                                                                        caAATTAAATTGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATg    <  1:148119/151‑1 (MQ=255)                                                                                                                                          aattaaattGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTc  <  2:77088/151‑1 (MQ=255)                                                                                                                                          aattaaattGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTc  <  2:716633/151‑1 (MQ=255)                                                                                                                                          aattaaattGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTc  <  2:499006/151‑1 (MQ=255)                                                                                                                                          aattaaattGGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTc  <  2:484423/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              CGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTTCTTGATGTTGACAAAATGACGCCGATTATTCTGGCTTTCCGCCAGGCAGGTTACGAAGTGGCAATTGATGATTTTGGTATTGGCTACTCTAACTTGCATAACCTTAAATCATTGAATGTC  >  NC_000913/4276500‑4276787

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