Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,350	Δ1 bp	coding (483/633 nt)	ycfQ ←	repressor for bhsA(ycfR)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,348	0	G	.	100.0%	145.6 / NA	35	coding (485/633 nt)	ycfQ	repressor for bhsA(ycfR)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (16/19);  total (16/19)

CGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGAT  >  NC_000913/1168233‑1168475                                                                                                                    |                                                                                                                                cGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGtt                                                                                                         <  1:58311/139‑1 (MQ=255)                        tttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                  <  2:971820/139‑1 (MQ=255)                        tttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAAGTAAGGAATTCTGCCAAATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGGTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                  <  2:771436/139‑1 (MQ=255)                                gATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTg                                                                          >  1:80690/1‑139 (MQ=255)                                 aTGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGc                                                                         <  1:534891/139‑1 (MQ=255)                                                            tCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTg                                              >  1:419244/1‑139 (MQ=255)                                                             ccccTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTgcgc                                                                               <  2:459551/105‑1 (MQ=255)                                                             ccccTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTgcgc                                                                               >  1:459551/1‑105 (MQ=255)                                                               ccTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGt                                                             <  1:680272/121‑1 (MQ=255)                                                               ccTCAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGAAAAGTGCGCTCTTGCATCGCATGGCGt                                                             >  2:680272/1‑121 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTa                                                 <  2:652965/124‑1 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTa                                                 >  1:652965/1‑124 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                  >  2:397573/1‑139 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                  >  2:935952/1‑139 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                  >  1:355240/1‑139 (MQ=255)                                                                         aCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTg                                 <  1:667809/139‑1 (MQ=255)                                                                            aTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAAt                              <  1:980178/139‑1 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCaaa                                                                                     <  2:1011717/83‑1 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCaaa                                                                                     >  1:1011717/1‑83 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATc                             >  1:1051533/1‑139 (MQ=255)                                                                              tAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCt                            <  2:383627/139‑1 (MQ=255)                                                                                     ttCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCg                                                                                                           >  2:555870/1‑53 (MQ=255)                                                                                     ttCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCg                                                                                                           <  1:555870/53‑1 (MQ=255)                                                                                     ttCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCaa                                                     <  2:173673/107‑1 (MQ=255)                                                                                     ttCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCaa                                                     >  1:173673/1‑107 (MQ=255)                                                                                        tGCCAGATGTGTCACGTGACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCTTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAggg                  <  2:767496/139‑1 (MQ=255)                                                                                          ccAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCCTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGggtg                >  1:689171/1‑139 (MQ=255)                                                                                           cAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGggtgg               <  1:58781/139‑1 (MQ=255)                                                                                             gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGgtggtg             <  1:981960/139‑1 (MQ=255)                                                                                             gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGgtggtg             <  1:397573/139‑1 (MQ=255)                                                                                                gtgtCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAg          <  2:1051533/139‑1 (MQ=255)                                                                                                gtgtCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAg          <  2:933030/139‑1 (MQ=255)                                                                                                   tcacgtcacAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAgttg       >  1:971908/1‑139 (MQ=255)                                                                                                    cacgtcacAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAgttgt      >  2:993771/1‑139 (MQ=255)                                                                                                        tcacAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGAt  >  2:7970/1‑139 (MQ=255)                                                                                                                    |                                                                                                                                CGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGAT  >  NC_000913/1168233‑1168475

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: ‑