Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,216,996	Δ2 bp		[ymgC]	[ymgC]

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,216,995	0	A	.	100.0%	276.1 / NA	64	coding (248/249 nt)	ymgC	Blue light, low temperature and stress induced protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (30/34);  total (30/34)
*	NC_000913	1,216,996	0	G	.	100.0%	276.1 / NA	64	coding (249/249 nt)	ymgC	Blue light, low temperature and stress induced protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (30/34);  total (30/34)

TGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCG  >  NC_000913/1216869‑1217120                                                                                                                               ||                                                                                                                             tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                                 >  2:608447/1‑139 (MQ=255)             aTGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATaaa                                                                                                     >  1:563612/1‑139 (MQ=255)              tGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                    <  1:187989/123‑1 (MQ=255)              tGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                    >  2:187989/1‑123 (MQ=255)               gAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACa                                                                                                   <  2:404081/139‑1 (MQ=255)                aaaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                  <  2:480090/139‑1 (MQ=255)                 aaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTCAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGACTAAGGTTTTAGAATAAACacc                                                                                                 <  2:808826/139‑3 (MQ=255)                      gtgtGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                            <  1:1194982/139‑1 (MQ=255)                                     aTCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                               >  1:1100461/1‑121 (MQ=255)                                     aTCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                               <  2:1100461/121‑1 (MQ=255)                                     aTCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGAt                                                                             <  1:418618/139‑1 (MQ=255)                                      tctcTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATg                                                                            >  1:887162/1‑139 (MQ=255)                                          tATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                                <  1:238752/131‑1 (MQ=255)                                          tATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                                >  2:238752/1‑131 (MQ=255)                                               gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   >  1:212985/1‑139 (MQ=255)                                               gCAAATAATTCAGGAAAATGAAATGCATAAGACCAACCCAACGGATATTATTAATAAAGACAGGGAATCCCTGCTCTCTT‑‑AACTTGAGTTAGGTTTTTGAAAAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  1:608447/139‑1 (MQ=255)                                                cAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGc                                                                  <  2:1341540/139‑1 (MQ=255)                                                    tatTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATg                                                                                  >  2:466940/1‑119 (MQ=255)                                                    tatTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATg                                                                                  <  1:466940/119‑1 (MQ=255)                                                      ttGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTaa                                                                                                                   >  2:694857/1‑84 (MQ=255)                                                      ttGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTaa                                                                                                                   <  1:694857/84‑1 (MQ=255)                                                               aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  2:1206602/139‑1 (MQ=255)                                                               aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  1:1161379/139‑1 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                    >  1:226180/1‑72 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                    <  2:226180/72‑1 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                 >  2:1206912/1‑139 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                 <  2:212985/139‑1 (MQ=255)                                                                   aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCttt                                               >  2:820781/1‑139 (MQ=255)                                                                        cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttat                                          >  2:1213906/1‑139 (MQ=255)                                                                         aTAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatt                                         <  2:887162/139‑1 (MQ=255)                                                                           aaGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTtattat                                       <  1:1206912/139‑1 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATg                                                                                  >  2:323430/1‑89 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATg                                                                                  <  1:323430/89‑1 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGa                                                                    >  2:420583/1‑103 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGa                                                                    <  1:420583/103‑1 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   >  1:262554/1‑99 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  2:262554/99‑1 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATc                           >  2:3442/1‑139 (MQ=255)                                                                                        cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                    <  2:998890/49‑1 (MQ=255)                                                                                        cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                    >  1:998890/1‑49 (MQ=255)                                                                                        cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCa                          >  1:876544/1‑139 (MQ=255)                                                                                           atatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgt                       <  1:820781/139‑1 (MQ=255)                                                                                             atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                     <  1:236129/139‑1 (MQ=255)                                                                                               cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  2:83040/107‑1 (MQ=255)                                                                                               cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   >  1:83040/1‑107 (MQ=255)                                                                                               cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa                   >  1:572858/1‑139 (MQ=255)                                                                                                 ttattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATa                                                                                                       >  2:963726/1‑53 (MQ=255)                                                                                                 ttattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATa                                                                                                       <  1:963726/53‑1 (MQ=255)                                                                                                          aCAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                  >  1:827559/1‑49 (MQ=255)                                                                                                          aCAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                  <  2:827559/49‑1 (MQ=255)                                                                                                          aCAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagag        >  1:1222300/1‑139 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTAt                                                                                         <  1:469845/56‑1 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTAt                                                                                         >  2:469845/1‑56 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTc                                                             >  2:863252/1‑84 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTc                                                             <  1:863252/84‑1 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCg                                  <  1:456533/111‑1 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCg                                  >  2:456533/1‑111 (MQ=255)                                                                                                              ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                                >  1:338014/1‑63 (MQ=255)                                                                                                              ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                                <  2:338014/63‑1 (MQ=255)                                                                                                              ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCttt                                               <  1:1232750/96‑1 (MQ=255)                                                                                                              ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCttt                                               >  2:1232750/1‑96 (MQ=255)                                                                                                                aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg  <  2:572858/139‑1 (MQ=255)                                                                                                                      tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   >  1:967381/1‑84 (MQ=255)                                                                                                                      tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  2:967381/84‑1 (MQ=255)                                                                                                                               ||                                                                                                                             TGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCG  >  NC_000913/1216869‑1217120

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