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%	233.9 / NA	56	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 . (27/29);  total (27/29)
*	NC_000913	1,216,996	0	G	.	100.0%	233.9 / NA	56	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 . (27/29);  total (27/29)

GCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTC  >  NC_000913/1216859‑1217123                                                                                                                                         ||                                                                                                                                gCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                                >  1:734379/1‑137 (MQ=255) gCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                                <  1:313223/139‑3 (MQ=255) gCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                                >  1:229133/1‑137 (MQ=255)           tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                                    >  1:1192086/1‑139 (MQ=255)           tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                                    >  1:1343837/1‑139 (MQ=255)                    aaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAAt                                                                                                           <  1:80024/139‑1 (MQ=255)                        tGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAAc                                                                                                       >  1:702850/1‑139 (MQ=255)                             ccAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                                  <  2:454762/139‑1 (MQ=255)                                 tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTaa                                                                                                                      <  1:689716/115‑1 (MQ=255)                                 tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTaa                                                                                                                      >  2:689716/1‑115 (MQ=255)                                 tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                    <  1:757340/133‑1 (MQ=255)                                 tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                    >  2:757340/1‑133 (MQ=255)                                     aTTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCa                                                                                          >  1:60101/1‑139 (MQ=255)                                               aTCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGAt                                                                                <  2:734379/139‑1 (MQ=255)                                                     aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaacttgag                                                                                                                           <  2:1631822/92‑8 (MQ=255)                                                     aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaacttgag                                                                                                                           >  1:1631822/1‑85 (MQ=255)                                                         gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                      >  1:214844/1‑139 (MQ=255)                                                         gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                      <  2:702850/139‑1 (MQ=255)                                                          cAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTaa                                                                                                                      >  2:688679/1‑90 (MQ=255)                                                          cAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTaa                                                                                                                      <  1:688679/90‑1 (MQ=255)                                                                  gCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAgg                                                             >  2:1781633/1‑139 (MQ=255)                                                                    aGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                           <  2:1192086/139‑1 (MQ=255)                                                                      tACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                     >  2:1259355/1‑95 (MQ=255)                                                                      tACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                     <  1:1259355/95‑1 (MQ=255)                                                                      tACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                      >  1:562084/1‑126 (MQ=255)                                                                      tACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                      <  2:562084/126‑1 (MQ=255)                                                                         aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      <  2:1343837/139‑1 (MQ=255)                                                                         aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      <  2:236414/139‑1 (MQ=255)                                                                           tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                    <  2:60101/139‑1 (MQ=255)                                                                             aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                       >  1:861162/1‑70 (MQ=255)                                                                             aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                       <  2:861162/70‑1 (MQ=255)                                                                                tGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCtt                                               >  2:1811500/1‑139 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                                <  1:798389/48‑3 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                                >  2:798389/1‑46 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      <  2:383272/120‑1 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      >  1:383272/1‑120 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                    >  2:868731/1‑122 (MQ=255)                                                                                            cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                                    <  1:868731/122‑1 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                       <  2:852532/50‑1 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                       >  1:852532/1‑50 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      <  1:1690971/115‑1 (MQ=255)                                                                                                 aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      >  2:1690971/1‑115 (MQ=255)                                                                                                   ggATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAg                            <  2:229133/139‑1 (MQ=255)                                                                                                      tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtg                         >  2:243046/1‑139 (MQ=255)                                                                                                       atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAgtgtc                        >  1:1486989/1‑135 (MQ=255)                                                                                                       atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAgtgtc                        >  1:1595254/1‑135 (MQ=255)                                                                                                         cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa                      >  2:1083373/1‑139 (MQ=255)                                                                                                         cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa                      <  1:1811500/139‑1 (MQ=255)                                                                                                              ttaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                   <  2:1574575/57‑1 (MQ=255)                                                                                                              ttaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                   >  1:1574575/1‑57 (MQ=255)                                                                                                                aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                         <  2:794946/97‑1 (MQ=255)                                                                                                                aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      >  2:1165800/1‑100 (MQ=255)                                                                                                                aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                      <  1:1165800/100‑1 (MQ=255)                                                                                                                aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGCTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGtt                                                         >  1:794946/1‑97 (MQ=255)                                                                                                                 aTGACAGGTAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                        <  2:1232548/129‑1 (MQ=255)                                                                                                                 aTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                        >  1:1232548/1‑129 (MQ=255)                                                                                                                  tGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgag             >  1:1458387/1‑139 (MQ=255)                                                                                                                   gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga              <  1:612102/137‑1 (MQ=255)                                                                                                                   gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga              >  2:612102/1‑137 (MQ=255)                                                                                                                        ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCA‑GTGTCAAGAATGGAGAGTtctc      <  2:1595254/139‑1 (MQ=255)                                                                                                                        ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCA‑GTGTCAAGAATGGAGAGTtctc      <  2:1486989/139‑1 (MQ=255)                                                                                                                          aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg     <  2:1458387/139‑1 (MQ=255)                                                                                                                             ccGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctc  <  1:1781633/139‑1 (MQ=255)                                                                                                                                         ||                                                                                                                                GCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTC  >  NC_000913/1216859‑1217123

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