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%	285.0 / NA	68	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 . (33/35);  total (33/35)
*	NC_000913	1,216,996	0	G	.	100.0%	285.0 / NA	68	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 . (33/35);  total (33/35)

GCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCC  >  NC_000913/1216859‑1217126                                                                                                                                         ||                                                                                                                                   gCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                                   <  1:246603/139‑3 (MQ=255)      aCGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaacttgag                                                                                                                              <  1:704036/139‑8 (MQ=255)        gCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          <  2:1210387/139‑1 (MQ=255)           tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                                       >  1:645002/1‑139 (MQ=255)              gTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGttt                                                                                                                    <  1:812292/139‑1 (MQ=255)                ttAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTa                                                                                                                  <  2:1110298/139‑1 (MQ=255)                                 tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                      <  1:140149/134‑1 (MQ=255)                                 tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                      >  2:140149/1‑134 (MQ=255)                                     aTTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCa                                                                                             >  2:1282027/1‑139 (MQ=255)                                         tttctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCAt                                                                                         >  2:686221/1‑139 (MQ=255)                                           tctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGc                                                                                       >  2:1090327/1‑139 (MQ=255)                                                         gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         <  1:1282027/139‑1 (MQ=255)                                                         gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCACTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     >  2:106697/1‑111 (MQ=255)                                                         gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCACTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     <  1:106697/111‑1 (MQ=255)                                                                     gTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGt                                                             >  2:596556/1‑139 (MQ=255)                                                                       aCATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTa                                                           <  2:645002/139‑1 (MQ=255)                                                                        cATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTAt                                                          <  1:686221/139‑1 (MQ=255)                                                                         aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                         >  2:933280/1‑139 (MQ=255)                                                                           tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         >  2:195553/1‑121 (MQ=255)                                                                           tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         <  1:1044834/121‑1 (MQ=255)                                                                           tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         >  2:1044834/1‑121 (MQ=255)                                                                           tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         <  1:195553/121‑1 (MQ=255)                                                                                tGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                                      <  1:1046768/103‑1 (MQ=255)                                                                                tGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                                      >  2:1046768/1‑103 (MQ=255)                                                                                  cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          >  2:449165/1‑65 (MQ=255)                                                                                  cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          <  1:449165/65‑1 (MQ=255)                                                                                       gACTACAGTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGTTTTTAGaa                                                                                                               >  1:894807/1‑71 (MQ=255)                                                                                       gACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGaa                                                                                                               <  2:894807/71‑1 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGATTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         >  2:924966/1‑104 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         <  1:924966/104‑1 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCa                                      >  1:645875/1‑139 (MQ=255)                                                                                                aaCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAatat                                  <  1:1347270/139‑1 (MQ=255)                                                                                                aaCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAatat                                  <  2:645875/139‑1 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     >  2:638751/1‑71 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     <  1:638751/71‑1 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATc                                 >  1:79329/1‑139 (MQ=255)                                                                                                  cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          >  1:427928/1‑49 (MQ=255)                                                                                                  cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          <  2:427928/49‑1 (MQ=255)                                                                                                  cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCa                                                                                          <  2:1070026/81‑1 (MQ=255)                                                                                                  cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCa                                                                                          >  1:1070026/1‑81 (MQ=255)                                                                                                  cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCa                                >  2:1038270/1‑139 (MQ=255)                                                                                                   ggATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAg                               <  1:596556/139‑1 (MQ=255)                                                                                                      tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          <  2:716689/45‑1 (MQ=255)                                                                                                      tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                          >  1:716689/1‑45 (MQ=255)                                                                                                      tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     <  1:571161/66‑1 (MQ=255)                                                                                                      tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     >  2:571161/1‑66 (MQ=255)                                                                                                       atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                            >  1:346859/1‑106 (MQ=255)                                                                                                       atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                            <  2:346859/106‑1 (MQ=255)                                                                                                       atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                           <  2:1231655/139‑1 (MQ=255)                                                                                                       atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                           <  1:1038270/139‑1 (MQ=255)                                                                                                         cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAAc                                                                                                          <  1:362352/58‑1 (MQ=255)                                                                                                         cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAAc                                                                                                          >  2:362352/1‑58 (MQ=255)                                                                                                               taATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                         <  1:250498/101‑1 (MQ=255)                                                                                                               taATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                         >  2:250498/1‑101 (MQ=255)                                                                                                                aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGTTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATgg                  <  1:662559/139‑1 (MQ=255)                                                                                                                   gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga               >  2:1297788/1‑139 (MQ=255)                                                                                                                    aCAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                         <  2:979187/80‑1 (MQ=255)                                                                                                                    aCAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGCGAg                                                                         >  1:979187/1‑80 (MQ=255)                                                                                                                      aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCtt                                                      <  2:418866/97‑1 (MQ=255)                                                                                                                      aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCtt                                                      >  1:418866/1‑97 (MQ=255)                                                                                                                      aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGtt            >  2:10812/1‑139 (MQ=255)                                                                                                                        ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     >  2:884086/1‑48 (MQ=255)                                                                                                                        ggAATCCGTGCTCGCTT‑‑AATTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                     <  1:884086/48‑1 (MQ=255)                                                                                                                          aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg        <  2:79329/139‑1 (MQ=255)                                                                                                                            tCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                              >  1:1306366/1‑83 (MQ=255)                                                                                                                            tCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                              <  2:1306366/83‑1 (MQ=255)                                                                                                                            tCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                              <  2:816025/83‑1 (MQ=255)                                                                                                                            tCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                              >  1:816025/1‑83 (MQ=255)                                                                                                                              cGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctct    <  1:1090327/139‑1 (MQ=255)                                                                                                                                tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTcc  >  2:717391/1‑139 (MQ=255)                                                                                                                                         ||                                                                                                                                   GCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCC  >  NC_000913/1216859‑1217126

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