Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,366,953	G→A	G10R (GGA→AGA)	yhcF →	putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,366,953	0	G	A	100.0%	143.5 / NA	41	G10R (GGA→AGA)	yhcF	putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/19);  total (22/19)

CTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCGGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCA  >  NC_000913/3366831‑3367084                                                                                                                           |                                                                                                                                    cTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                     <  2:566254/139‑1 (MQ=255)     cacTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATg                                                                                                                 >  2:75574/1‑139 (MQ=255)          aCAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgc                                                                                                            >  1:591930/1‑139 (MQ=255)                gCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                      <  2:52505/139‑1 (MQ=255)                        ttttAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATa                                                                                              <  2:303852/139‑1 (MQ=255)                                   aaGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                      >  1:451625/1‑120 (MQ=255)                                   aaGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                      <  2:451625/120‑1 (MQ=255)                                                           atatTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCt                                                           >  1:228180/1‑139 (MQ=255)                                                           atatTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAACCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCt                                                           <  2:449829/139‑1 (MQ=255)                                                                aCAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGc                                                      >  2:606269/1‑139 (MQ=255)                                                                 cAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATa                                                                                              >  2:143857/1‑98 (MQ=255)                                                                 cAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATa                                                                                              <  1:143857/98‑1 (MQ=255)                                                                   ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGACGCTGCAGCTTg                                                   <  2:76082/139‑1 (MQ=255)                                                                   ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGTAGCTTg                                                   <  1:15217/139‑1 (MQ=255)                                                                   ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   <  1:606269/139‑1 (MQ=255)                                                                   ggCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   <  1:75574/139‑1 (MQ=255)                                                                     cAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTa                                                 >  2:241913/1‑139 (MQ=255)                                                                        tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTg                                                                                  <  1:82603/103‑1 (MQ=255)                                                                        tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTg                                                                                  >  2:82603/1‑103 (MQ=255)                                                                        tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   <  2:52720/134‑1 (MQ=255)                                                                        tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   >  1:52720/1‑134 (MQ=255)                                                                            aCGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAAt                                          >  2:601803/1‑139 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgcgc                                                                                                          <  2:484529/61‑1 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgcgc                                                                                                          >  1:484529/1‑61 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   <  2:568377/116‑1 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   <  1:131608/116‑1 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   >  2:131608/1‑116 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   >  1:568377/1‑116 (MQ=255)                                                                                          ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAAcggcg                            >  2:270776/1‑139 (MQ=255)                                                                                           gCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAAcggcgg                           >  2:495281/1‑139 (MQ=255)                                                                                               aaTGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGAc                                                                                <  2:228406/82‑1 (MQ=255)                                                                                               aaTGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGAc                                                                                >  1:228406/1‑82 (MQ=255)                                                                                                       aTGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACtt               >  2:200839/1‑139 (MQ=255)                                                                                                       aTGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACtt               >  1:130583/1‑139 (MQ=255)                                                                                                                  cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCaa                                                                                                        <  2:233858/39‑1 (MQ=255)                                                                                                                  cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCaa                                                                                                        >  1:233858/1‑39 (MQ=255)                                                                                                                  cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCt                                                           <  2:534688/84‑1 (MQ=255)                                                                                                                  cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCt                                                           >  1:534688/1‑84 (MQ=255)                                                                                                                    gATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   >  2:254995/1‑90 (MQ=255)                                                                                                                    gATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                   <  1:254995/90‑1 (MQ=255)                                                                                                                    gATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCa  >  1:500275/1‑139 (MQ=255)                                                                                                                           |                                                                                                                                    CTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCGGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATTCA  >  NC_000913/3366831‑3367084

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