Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,217,327	A→C	intergenic (+331/‑1)	ymgC → / → ycgG	Blue light, low temperature and stress induced protein/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	1,217,327	0	A	C	100.0%	185.9 / NA	51	intergenic (+331/‑1)	ymgC/ycgG	Blue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (23/28);  total (23/28)

AATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCAAATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATATCAATCATATACT  >  NC_000913/1217185‑1217473                                                                                                                                               |                                                                                                                                                   aattaatGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCa                                                                                                                                             <  1:385606/150‑1 (MQ=255)                     aTGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCg                                                                                                                        <  2:755415/151‑1 (MQ=255)                         agaAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGc                                                                                                                    <  2:854731/151‑1 (MQ=255)                               ttttATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATg                                                                                                              >  2:652140/1‑151 (MQ=255)                                      ttCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGttt                                                                                                       >  1:379166/1‑151 (MQ=255)                                       tCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTg                                                                                                          <  1:872836/147‑1 (MQ=255)                                       tCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGttt                                                                                                       <  1:697896/150‑1 (MQ=255)                                        cGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTtta                                                                                                      <  2:852385/150‑1 (MQ=255)                                        cGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGAAGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTttat                                                                                                     <  2:72516/151‑1 (MQ=255)                                         ggCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTttatt                                                                                                    <  2:839110/151‑1 (MQ=255)                                          gCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTttatt                                                                                                    >  1:901575/1‑150 (MQ=255)                                           cGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTAc                                                                                                  >  1:228282/1‑151 (MQ=255)                                                 tatTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGgtg                                                                                             <  2:920230/150‑1 (MQ=255)                                                  atTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGgtgt                                                                                            <  1:17331/150‑1 (MQ=255)                                                  atTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTc                                                                                           >  1:211635/1‑151 (MQ=255)                                                  atTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTc                                                                                           <  2:490871/151‑1 (MQ=255)                                                          tGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTAtttt                                                                                   >  1:891332/1‑151 (MQ=255)                                                                ctACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAAc                                                                               <  2:122312/149‑1 (MQ=255)                                                                  aCCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACAtac                                                                           <  2:488761/151‑3 (MQ=255)                                                                  aCCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATcc                                                                           <  1:129895/151‑1 (MQ=255)                                                                      tGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGct                                                                       >  1:471336/1‑151 (MQ=255)                                                                      tGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGct                                                                       >  1:650232/1‑151 (MQ=255)                                                                            ctGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGt                                                                  >  2:43641/1‑150 (MQ=255)                                                                              gTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGt                                                                  <  1:261069/148‑1 (MQ=255)                                                                              gTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTAtt                                                               >  1:921917/1‑151 (MQ=255)                                                                                aGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTcc                                                             >  1:543212/1‑151 (MQ=255)                                                                                 gATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCg                                                            <  1:115554/151‑1 (MQ=255)                                                                                     ccAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCaa                                                         >  1:398679/1‑150 (MQ=255)                                                                                         cgcgAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCa                                                    <  2:290756/151‑1 (MQ=255)                                                                                          gcgagcgaATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCa                                                    <  2:348565/150‑1 (MQ=255)                                                                                          gcgagcgaATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCa                                                    <  2:375075/150‑1 (MQ=255)                                                                                              gcgaATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGt                                                >  1:645172/1‑150 (MQ=255)                                                                                                      atgGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGTTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCgg                                       >  1:183388/1‑151 (MQ=255)                                                                                                        ggCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGa                                     <  2:564743/151‑1 (MQ=255)                                                                                                         gCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCggg                                      >  2:411323/1‑149 (MQ=255)                                                                                                             tttCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCg                                  >  1:105707/1‑149 (MQ=255)                                                                                                               tCATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCgaga                               <  2:219215/150‑1 (MQ=255)                                                                                                                cATGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGat                              <  2:450392/150‑1 (MQ=255)                                                                                                                 aTGTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGatat                            <  2:430473/151‑1 (MQ=255)                                                                                                                   gTATTAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATgc                          >  2:302111/1‑151 (MQ=255)                                                                                                                      ttAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATgccgc                       >  1:831492/1‑151 (MQ=255)                                                                                                                       tAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATgccgc                       <  1:611984/150‑1 (MQ=255)                                                                                                                       tAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATgccgc                       <  1:684089/150‑1 (MQ=255)                                                                                                                       tAAACCCTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATgccgc                       >  1:200915/1‑150 (MQ=255)                                                                                                                           cccTTTGGCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCcaac                   >  1:828188/1‑150 (MQ=255)                                                                                                                                  gCCGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATatcaa           <  2:753707/151‑1 (MQ=255)                                                                                                                                    cGGGACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATatcaatc         <  2:395663/151‑1 (MQ=255)                                                                                                                                      ggACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATATCAATCat       >  1:603852/1‑151 (MQ=255)                                                                                                                                       gACAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATATCAATCata      <  2:500156/151‑1 (MQ=255)                                                                                                                                         cAAGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATATCAATCatat     >  2:689978/1‑150 (MQ=255)                                                                                                                                           aGCACATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATATCAATCATATACt  >  2:380818/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                   AATTAATGGGTTAGATCAAGATGGAGAAAATTTTATCTTCGGCGTGGATATTTCTTATGCTCTCTACCATGCTCTCTGTAGATGCCAGCGCGAGCGAATGAATGGCATTTTCATGTATTAAACCCTTTGGCCGGGACAAGCAAATGCGCAATACACTCATACCCATCCTCGTTGCGATATGCTTGTTTATTACCGGTGTCGCTATTTTGAACATCCAGCTCTGGTATTCCGCAAAAGCAGAGTACCTGGCGGGAGCGAGATATGCCGCCAACAATATCAATCATATACT  >  NC_000913/1217185‑1217473

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