Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,061,157	G→A	*237* (TAG→TAA)	yihL →	putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,157	0	G	A	100.0%	133.8 / NA	42	*237* (TAG→TAA)	yihL	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/23);  total (19/23)

TCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTG  >  NC_000913/4061028‑4061288                                                                                                                                  |                                                                                                                                    tCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTAt                                                                                                                            >  2:215773/1‑139 (MQ=255)   ttGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:303543/139‑1 (MQ=255)    tGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATATGTAAGTGTTTTATGGt                                                                                                                         <  2:751561/139‑1 (MQ=255)                      gCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                                       >  2:27943/1‑139 (MQ=255)                      gCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                                       <  1:391215/139‑1 (MQ=255)                        cAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                     <  2:360065/139‑1 (MQ=255)                         aGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                                    <  2:430767/139‑1 (MQ=255)                                       aaaaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAAAgtg                                                                                      >  2:578602/1‑139 (MQ=255)                                           gCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                                  <  1:27943/139‑1 (MQ=255)                                                   ccGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                          <  1:669194/139‑1 (MQ=255)                                                       ccAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCctct                                                                      <  1:7407/139‑1 (MQ=255)                                                                  atatTGTCCAAATAGATATTGAATATTGGCGATTTGCGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                                    >  2:641070/1‑98 (MQ=255)                                                                  atatTGTCCAAATAGATATTGAATATTGGCGATTTGCGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                                    <  1:641070/98‑1 (MQ=255)                                                                         ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  2:438658/1‑69 (MQ=255)                                                                         ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  1:438658/69‑1 (MQ=255)                                                                             aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          >  1:270465/1‑65 (MQ=255)                                                                             aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                          <  2:270465/65‑1 (MQ=255)                                                                             aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                               >  2:499417/1‑124 (MQ=255)                                                                             aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                               <  1:499417/124‑1 (MQ=255)                                                                                  tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCt                                                         <  1:494221/125‑1 (MQ=255)                                                                                  tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCt                                                         >  2:494221/1‑125 (MQ=255)                                                                                       aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAata                                      <  1:215773/139‑1 (MQ=255)                                                                                              gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAAc                                                                                                            >  1:103899/1‑62 (MQ=255)                                                                                              gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAAc                                                                                                            <  2:103899/62‑1 (MQ=255)                                                                                              gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCt                               >  2:296801/1‑139 (MQ=255)                                                                                              gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCt                               <  1:578602/139‑1 (MQ=255)                                                                                                   ttGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                              <  2:515099/103‑1 (MQ=255)                                                                                                   ttGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                    >  1:608207/1‑129 (MQ=255)                                                                                                   ttGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                    <  2:608207/129‑1 (MQ=255)                                                                                                   atGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATAGTGACAATAAATAACGCAAGAAAGATTCTACTACGTGTCGACACTCTTCCTCTTTATTTAc                                                              >  1:515099/2‑103 (MQ=255)                                                                                                            tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                          <  1:729309/82‑1 (MQ=255)                                                                                                            tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                          >  2:729309/1‑82 (MQ=255)                                                                                                                    atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAAc                                                                                         <  2:524359/59‑1 (MQ=255)                                                                                                                    atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACActc                                                                             <  2:164634/71‑1 (MQ=255)                                                                                                                    atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACActc                                                                             >  1:164634/1‑71 (MQ=255)                                                                                                                    atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt         >  2:467913/1‑139 (MQ=255)                                                                                                                    atcatTAATCTGTAAGTATTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAAc                                                                                         >  1:524359/1‑59 (MQ=255)                                                                                                                      catTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtgt                                                                                     >  2:225555/1‑61 (MQ=255)                                                                                                                      catTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtgt                                                                                     <  1:225555/61‑1 (MQ=255)                                                                                                                         tAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgc             >  2:596985/1‑130 (MQ=255)                                                                                                                         tAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgc             <  1:596985/130‑1 (MQ=255)                                                                                                                           aTCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTg  >  2:312897/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                    TCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTG  >  NC_000913/4061028‑4061288

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