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

TGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCG  >  NC_000913/1216869‑1217120                                                                                                                               ||                                                                                                                             tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                                 >  1:312940/1‑139 (MQ=255) tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                                 >  2:171149/1‑139 (MQ=255)     tttAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGtttt                                                                                                             >  2:841659/1‑139 (MQ=255)         aaaGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGaa                                                                                                         <  1:841659/139‑1 (MQ=255)               gAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACa                                                                                                   >  1:11973/1‑139 (MQ=255)                aaaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                  >  1:64826/1‑139 (MQ=255)                 aaCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                 <  2:64826/139‑1 (MQ=255)                  aCCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                <  2:358061/139‑1 (MQ=255)                      gtgtGATTCTTTCTGATCTCAATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                            >  1:1335149/1‑139 (MQ=255)                       tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGt                                                                                           >  1:1403159/1‑139 (MQ=255)                        gtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTa                                                                                          <  2:312940/139‑1 (MQ=255)                                 tctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGc                                                                                 >  1:1627212/1‑139 (MQ=255)                                           aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                             >  1:318395/1‑85 (MQ=255)                                           aTCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                             <  2:318395/87‑3 (MQ=255)                                             cTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTg                                                                     <  2:1335149/139‑1 (MQ=255)                                               gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaa                                                                                                                              <  2:1097262/82‑2 (MQ=255)                                               gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaa                                                                                                                              >  1:1097262/1‑81 (MQ=255)                                               gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  2:11973/139‑1 (MQ=255)                                                      ttGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCa                                                            >  1:480220/1‑139 (MQ=255)                                                        gCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaacttg                                                                                                                          <  1:1419848/77‑6 (MQ=255)                                                        gCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaacttg                                                                                                                          >  2:1419848/1‑72 (MQ=255)                                                               aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  1:171149/139‑1 (MQ=255)                                                               aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  2:1627212/139‑1 (MQ=255)                                                               aTTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  2:480220/139‑1 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACTGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                               <  2:1036755/93‑1 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                               >  1:1036755/1‑93 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  2:839642/137‑1 (MQ=255)                                                                 tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   >  1:839642/1‑137 (MQ=255)                                                                        cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCaaa                                                                        >  1:1052985/1‑109 (MQ=255)                                                                        cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCaaa                                                                        <  2:1052985/109‑1 (MQ=255)                                                                          tAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  1:926500/128‑1 (MQ=255)                                                                          tAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   >  2:926500/1‑128 (MQ=255)                                                                             gACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  1:1468454/125‑1 (MQ=255)                                                                             gACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   >  2:1468454/1‑125 (MQ=255)                                                                                tACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                        >  2:232501/1‑117 (MQ=255)                                                                                tACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCg                                                        <  1:232501/117‑1 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGc                                                                                              <  2:406342/77‑1 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGc                                                                                              >  1:406342/1‑77 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   >  2:944179/1‑120 (MQ=255)                                                                                  cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTAGCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                                   <  1:944179/120‑1 (MQ=255)                                                                                     tAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   >  2:693101/1‑101 (MQ=255)                                                                                     tAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  1:693101/101‑1 (MQ=255)                                                                                      aaCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAatat                            <  2:1403159/139‑1 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  1:1271453/99‑1 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  1:1257314/99‑1 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   >  2:1257314/1‑99 (MQ=255)                                                                                       aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   >  2:1271453/1‑99 (MQ=255)                                                                                        cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                 <  2:72938/68‑1 (MQ=255)                                                                                        cGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                 >  1:72938/1‑68 (MQ=255)                                                                                            tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGAt                                                                             <  1:389557/84‑1 (MQ=255)                                                                                            tatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGAt                                                                             >  2:389557/1‑84 (MQ=255)                                                                                                  tattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCAc                                                                                      <  2:787906/69‑1 (MQ=255)                                                                                                  tattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCAc                                                                                      >  1:787906/1‑69 (MQ=255)                                                                                                         gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                <  1:736004/52‑1 (MQ=255)                                                                                                         gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                                >  2:736004/1‑52 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                 >  1:1659635/1‑48 (MQ=255)                                                                                                            aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                                 <  2:1659635/48‑1 (MQ=255)                                                                                                              ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTa                                                                                          >  1:1196131/1‑53 (MQ=255)                                                                                                              ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTa                                                                                          <  2:1196131/53‑1 (MQ=255)                                                                                                                  tCCGTGCCCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   >  2:1125022/1‑72 (MQ=255)                                                                                                                  tCCGTGCCCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                                   <  1:1125022/72‑1 (MQ=255)                                                                                                                      tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                  <  1:29698/37‑1 (MQ=255)                                                                                                                      tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACAt                                                                                                  >  2:29698/1‑37 (MQ=255)                                                                                                                      tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg  >  1:1228943/1‑133 (MQ=255)                                                                                                                      tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCg  <  2:1228943/133‑1 (MQ=255)                                                                                                                               ||                                                                                                                             TGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCG  >  NC_000913/1216869‑1217120

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