Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,217,003	A→G	intergenic (+7/‑325)	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,003	0	A	G	100.0%	216.6 / NA	57	intergenic (+7/‑325)	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 G (26/31);  total (26/31)

CAGGCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCC  >  NC_000913/1216856‑1217147                                                                                                                                                    |                                                                                                                                                 cAGGCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTa                                                                                                                                               <  1:459798/151‑1 (MQ=255)   ggCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAg                                                                                                                                             <  1:575576/151‑1 (MQ=255)        tACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGttt                                                                                                                                         <  2:437002/150‑1 (MQ=255)        tACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGttt                                                                                                                                         <  1:345138/150‑1 (MQ=255)               gCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATTTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATa                                                                                                                                  <  2:583242/150‑1 (MQ=255)                cGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATaa                                                                                                                                 >  2:413368/1‑150 (MQ=255)                cGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATaa                                                                                                                                 >  2:561713/1‑150 (MQ=255)                cGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATCAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATaa                                                                                                                                 >  2:196764/1‑150 (MQ=255)                    tAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACAt                                                                                                                             <  1:519281/150‑1 (MQ=255)                    tAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATc                                                                                                                            >  2:341870/1‑151 (MQ=255)                     aaaaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATc                                                                                                                            <  2:763352/150‑1 (MQ=255)                      aaaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATc                                                                                                                            >  2:320512/1‑149 (MQ=255)                            gAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTAt                                                                                                                    >  1:269881/1‑151 (MQ=255)                              aaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATc                                                                                                                   >  2:54749/1‑150 (MQ=255)                                   gtgtGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATg                                                                                                             <  2:568504/151‑1 (MQ=255)                                          tctttctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATg                                                                                                       <  1:291185/150‑1 (MQ=255)                                           ctttctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGc                                                                                                      >  1:476842/1‑150 (MQ=255)                                                 gATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTg                                                                                                <  2:748821/150‑1 (MQ=255)                                                 gATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTg                                                                                                <  2:10901/150‑1 (MQ=255)                                                       tATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAttt                                                                                         >  2:332375/1‑151 (MQ=255)                                                        aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAttt                                                                                         <  1:443422/150‑1 (MQ=255)                                                            gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGc                                                                                    >  2:55689/1‑151 (MQ=255)                                                               aaTATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                                                 <  1:499301/151‑1 (MQ=255)                                                                     gCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCtt                                                                           >  1:428278/1‑151 (MQ=255)                                                                      cAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCttt                                                                          <  1:574628/151‑1 (MQ=255)                                                                       aGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTc                                                                         <  1:6639/151‑1 (MQ=255)                                                                          aCATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCtt                                                                       <  2:221711/150‑1 (MQ=255)                                                                          aCATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCtt                                                                       <  1:346689/150‑1 (MQ=255)                                                                            aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatt                                                                    <  1:37448/151‑1 (MQ=255)                                                                            aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatt                                                                    >  1:306902/1‑151 (MQ=255)                                                                              tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatta                                                                   <  2:346482/150‑1 (MQ=255)                                                                               gAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTtattat                                                                  <  1:122772/150‑1 (MQ=255)                                                                               gAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTtattat                                                                  <  2:736405/150‑1 (MQ=255)                                                                                ataTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGCGCATTTCAGGCGTTATGCTTTCTTATTATGt                                                                >  1:770037/3‑151 (MQ=255)                                                                                 aaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTc                                                               >  1:557313/1‑151 (MQ=255)                                                                                 aaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTc                                                               <  1:554277/151‑1 (MQ=255)                                                                                  aTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                                                              <  2:574302/151‑1 (MQ=255)                                                                                     cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCa                                                           <  2:176252/151‑1 (MQ=255)                                                                                      aTAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCa                                                           >  2:94202/1‑150 (MQ=255)                                                                                         aGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAatat                                                       <  1:359351/151‑1 (MQ=255)                                                                                                   aaCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa                                              <  1:184261/150‑1 (MQ=255)                                                                                                    aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAg                                            <  1:650256/151‑1 (MQ=255)                                                                                                       gATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAAt                                         >  1:715537/1‑151 (MQ=255)                                                                                                              ttattaATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGt                                  >  1:294833/1‑151 (MQ=255)                                                                                                               tattaATGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGtt                                 >  2:697237/1‑151 (MQ=255)                                                                                                                   aaTGACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg                             >  1:542064/1‑151 (MQ=255)                                                                                                                      gACAGGGAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctc                          >  1:623443/1‑151 (MQ=255)                                                                                                                          gggAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCa                      <  1:76930/151‑1 (MQ=255)                                                                                                                           ggAATCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAt                     >  2:488263/1‑151 (MQ=255)                                                                                                                              aTCCGTGCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                  >  2:10094/1‑151 (MQ=255)                                                                                                                                    gCTCTCTTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGc            >  2:74273/1‑151 (MQ=255)                                                                                                                                      tctctTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt          >  1:459075/1‑151 (MQ=255)                                                                                                                                       ctctTAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGcc           <  1:209088/149‑1 (MQ=255)                                                                                                                                           tAGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat       <  2:579445/149‑1 (MQ=255)                                                                                                                                            aGAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCTTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGatat     <  1:28630/150‑1 (MQ=255)                                                                                                                                             gAACTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATcc   >  2:697324/1‑151 (MQ=255)                                                                                                                                              aaCTTGGGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATccc  >  1:642334/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CAGGCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCC  >  NC_000913/1216856‑1217147

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