Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,129,647	T→C	Q190R (CAG→CGG)	yghO ←	putative DNA‑binding transcriptional regulator YghO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,129,647	0	T	C	100.0%	132.8 / NA	38	Q190R (CAG→CGG)	yghO	putative DNA‑binding transcriptional regulator YghO
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/22);  total (16/22)

ATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCC  >  NC_000913/3129524‑3129774                                                                                                                            |                                                                                                                                aTCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATccacc                                                                                                                        <  1:601417/133‑1 (MQ=255)  tctcACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCTCCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATccacca                                                                                                                       <  2:593496/133‑1 (MQ=255)    tcACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACg                                                                                                                     >  1:134597/1‑133 (MQ=255)     cACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGc                                                                                                                    >  2:7365/1‑133 (MQ=255)        gCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAg                                                                                                                 <  2:514245/133‑1 (MQ=255)         cAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGt                                                                                                                <  1:7365/133‑1 (MQ=255)            gATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTaaa                                                                                                             >  2:587810/1‑133 (MQ=255)                 gCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAAt                                                                                                        <  2:134597/133‑1 (MQ=255)                  cATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAATGAGATCGGTTCTCCGCATCCACCATGCCAATAAATCAATg                                                                                                       >  2:247786/1‑133 (MQ=255)                      tctTCGGCCAACCGCCGACGATGGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCtt                                                                                                   <  1:226376/133‑1 (MQ=255)                      tctTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCtt                                                                                                   >  2:524373/1‑133 (MQ=255)                      tctTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCtt                                                                                                   >  1:501431/1‑133 (MQ=255)                      tctTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCtt                                                                                                   >  1:510893/1‑133 (MQ=255)                          cGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTtgtg                                                                                               >  2:141945/1‑133 (MQ=255)                          cGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTtgtg                                                                                               >  2:496559/1‑133 (MQ=255)                                 ccGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAAccc                                                                                        >  2:341753/1‑133 (MQ=255)                                  cGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCa                                                                                       <  2:17095/133‑1 (MQ=255)                                   gCTGACGATTAATGCTGCGAATGGCCACGTTTTTGCGCACCGGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAg                                                                                      <  2:199698/133‑1 (MQ=255)                                                 cAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATgcg                                                                        <  1:340760/133‑1 (MQ=255)                                                 cAGCGAATGGTAACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATgcg                                                                        <  2:427626/133‑1 (MQ=255)                                                   gCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCGAGATGTTCAATATGcgcg                                                                      >  1:238853/1‑133 (MQ=255)                                                   gCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGcgcg                                                                      <  1:341753/133‑1 (MQ=255)                                                            cACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGGCTTTGTGAAAACCCAGTTGTTCAATATGCGGGGCATACCAc                                                             >  2:199198/1‑133 (MQ=255)                                                             aCCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAAGCCAGTTGCTCACTAGGCGCGGAATACCACg                                                            >  1:280697/1‑133 (MQ=255)                                                                   ttGCGCACCTCATCCATCAGTTTTCTTAGCGCCGGAGAGAATGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGc                                                      <  2:204823/133‑1 (MQ=255)                                                                          ccTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTgg                                               <  1:587810/133‑1 (MQ=255)                                                                          ccTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTgg                                               <  2:547619/133‑1 (MQ=255)                                                                               tCCATCAGTTTTTTTGGCGCCGGATAGAAAGTAAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGcatca                                          <  2:257029/133‑1 (MQ=255)                                                                               tCCATCAGTTTTTTTCGCTCCGCAGAGAAATTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGcatca                                          <  1:175113/133‑1 (MQ=255)                                                                                 cATCAGTTTTTTTAGTGCTGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCatcatc                                        <  1:314211/133‑1 (MQ=255)                                                                                   tCAGTTTTTTTAGCGCCGGAGATAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGc                                      <  1:272897/133‑1 (MQ=255)                                                                                   tCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGc                                      <  1:524373/133‑1 (MQ=255)                                                                                    cAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGca                                     <  2:510893/133‑1 (MQ=255)                                                                                                    ggAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaa                     <  2:479385/133‑1 (MQ=255)                                                                                                    ggAGAGAAAGTGAGATCGGTTCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaa                     <  2:501431/133‑1 (MQ=255)                                                                                                                       ttCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATcc  >  1:426905/1‑133 (MQ=255)                                                                                                                       ttCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATcc  >  1:454014/1‑133 (MQ=255)                                                                                                                       ttCGCCGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATcc  >  1:640202/1‑133 (MQ=255)                                                                                                                            |                                                                                                                                ATCTCACGCAGGATCTGCATCTCTTCGGCAAACCGCTGACGATTGATGCAGCGAATGGTCACCTTTTTGCGCACCTGATCCATCAGTTTTTTTAGCGCCGGAGAGAAAGTGAGATCGGTTCGCTGCATCCACCACGCCAGTAAATCAATGCCTTTGTGATAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCC  >  NC_000913/3129524‑3129774

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