Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	357,606	C→T	intergenic (+152/+185)	codA → / ← cynR	cytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,606	0	C	T	100.0%	70.7 / NA	31	intergenic (+152/+185)	codA/cynR	cytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/14);  total (17/14)

GCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGC  >  NC_000913/357473‑357739                                                                                                                                      |                                                                                                                                      gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGcc                                                                                                                                  >  2:1024844/1‑139 (MQ=21)          aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcgc                                                                                                                         >  1:513371/1‑139 (MQ=16)              cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                     <  2:513371/139‑1 (MQ=12)                                                  agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCcaaca                                                                                                                              >  1:190701/1‑94 (MQ=255)                                                  agCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCcaaca                                                                                                                              <  2:190701/94‑1 (MQ=255)                                                       tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGt                                                                            >  1:837035/1‑139 (MQ=255)                                                       tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATT‑CAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                          <  2:754499/92‑1 (MQ=11)                                                       tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATT‑CAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAAAAACcg                                                                                                                          >  1:754499/1‑92 (MQ=11)                                                                   ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTa                                                                >  1:351448/1‑139 (MQ=11)                                                                      gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGtt                                                             >  1:832722/1‑139 (MQ=11)                                                                              tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                     <  1:1007927/75‑1 (MQ=255)                                                                              tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                                     >  2:1007927/1‑75 (MQ=255)                                                                              tAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCgg                                                     >  2:730977/1‑139 (MQ=14)                                                                               aGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGa                                                    <  2:362167/139‑1 (MQ=14)                                                                                 cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATg                                                  <  1:27247/139‑1 (MQ=255)                                                                                 cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATAGTTTATGCCgaacg                                                  >  2:595368/1‑135 (MQ=11)                                                                                     tCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                          >  2:236414/1‑63 (MQ=14)                                                                                     tCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg                                                                                                                          <  1:236414/63‑1 (MQ=14)                                                                                      cAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcggcg                                             >  2:134580/1‑139 (MQ=17)                                                                                        ggCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                           <  1:65756/139‑1 (MQ=17)                                                                                        ggCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                                           <  2:790138/139‑1 (MQ=17)                                                                                             ttCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGc                                      <  2:986025/139‑1 (MQ=255)                                                                                                              tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa                     >  1:15402/1‑139 (MQ=18)                                                                                                                   gCGGTTACCGTGATTCGTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGCGCa                <  1:730977/139‑1 (MQ=18)                                                                                                                   gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                       <  2:475642/132‑1 (MQ=18)                                                                                                                   gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg                       >  1:475642/1‑132 (MQ=18)                                                                                                                   gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa                <  2:581310/139‑1 (MQ=21)                                                                                                                   gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa                >  2:968155/1‑139 (MQ=21)                                                                                                                   gCGGTTACCGTGATTCATTTCCGCCAACAAC‑‑CGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCCGATGCAGCGTGAACGCCTTATCCGGCCTACGTACACCACt              >  2:406810/1‑139 (MQ=18)                                                                                                                      gttACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACACCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTg             >  1:766981/3‑139 (MQ=21)                                                                                                                            gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCg       <  2:416598/139‑1 (MQ=255)                                                                                                                                 tcatttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGc  <  2:351448/139‑1 (MQ=255)                                                                                                                                 tcatttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGc  <  1:968155/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      GCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGC  >  NC_000913/357473‑357739

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