Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	679,874	G→A	G123R (GGG→AGG)	ybeU →	DUF1266 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	679,874	0	G	A	100.0%	127.5 / NA	41	G123R (GGG→AGG)	ybeU	DUF1266 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/20);  total (21/20)

CGCTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCAC  >  NC_000913/679739‑679987                                                                                                                                        |                                                                                                                  cGCTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGt                                                                                                                <  1:68639/139‑1 (MQ=255)              aaGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGTTTATGTTCGAATGAGGTTTCTGTCACGGGt                                                                                                   >  2:404152/1‑139 (MQ=255)              aaGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGt                                                                                                   >  2:579107/1‑139 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                              >  2:360334/1‑139 (MQ=255)                   tGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGgtg                                                                                              >  1:743828/1‑139 (MQ=255)                        aGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTc                                                                                         >  1:20547/1‑139 (MQ=255)                        aGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTc                                                                                         >  2:86703/1‑139 (MQ=255)                             tAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTcaacaa                                                                                    >  1:216055/1‑139 (MQ=255)                                     ggtggCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGa                                                                            <  1:469687/139‑1 (MQ=255)                                       tggCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACg                                                                          <  1:538007/139‑1 (MQ=255)                                                   agGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACg                                                                                                      >  2:809078/1‑99 (MQ=255)                                                   agGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGGGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACg                                                                                                      <  1:809078/99‑1 (MQ=255)                                                                    aaTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                             >  2:425034/1‑139 (MQ=255)                                                                     aTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcgc                                            <  1:594764/139‑1 (MQ=255)                                                                     aTTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcgc                                            >  1:282435/1‑139 (MQ=255)                                                                      ttCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                             >  2:619799/1‑137 (MQ=255)                                                                      ttCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                             <  1:619799/137‑1 (MQ=255)                                                                        cGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTc                                         <  2:282435/139‑1 (MQ=255)                                                                              ccAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTa                                   <  1:254346/139‑1 (MQ=255)                                                                                    cGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCt                             <  1:425034/139‑1 (MQ=255)                                                                                       tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                           >  2:37238/1‑138 (MQ=255)                                                                                       tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCa                           <  1:37238/138‑1 (MQ=255)                                                                                       tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAc                          <  1:438328/139‑1 (MQ=255)                                                                                       tATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCCCAc                          >  2:827718/1‑139 (MQ=255)                                                                                        aTGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                         <  1:65334/139‑1 (MQ=255)                                                                                          gtgCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCAt                       <  1:827718/139‑1 (MQ=255)                                                                                              tgtgGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACt                   >  1:81692/1‑139 (MQ=255)                                                                                                     aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGc                              >  1:787965/1‑121 (MQ=255)                                                                                                     aaGGGGGAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGc                              <  2:787965/121‑1 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGttt                                                            <  1:421237/86‑1 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGttt                                                            >  2:421237/1‑86 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                         <  1:372188/121‑1 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCAcc                         >  2:372188/1‑121 (MQ=255)                                                                                                          ggAATCAAAGCCTGGGATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGa       <  2:216055/139‑1 (MQ=255)                                                                                                               cAAAGCCTGGGATTAAGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACCAGTGGCTGACGGAAGAAGTTAGTTTATGGCTGCAATCGCGCGTATATCTACGAGCTCACCATTACTATAACAGCTAGATGCAc  >  2:183067/1‑139 (MQ=255)                                                                                                                        ggATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCt                                                                                          >  1:837639/1‑42 (MQ=255)                                                                                                                        ggATTATGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCt                                                                                          <  2:837639/42‑1 (MQ=255)                                                                                                                              tGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                             >  2:765278/1‑81 (MQ=255)                                                                                                                              tGTTCGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATcgcg                                             <  1:765278/81‑1 (MQ=255)                                                                                                                                  cGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTAt                                                          >  1:183291/1‑64 (MQ=255)                                                                                                                                  cGAATGAGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTAt                                                          <  2:183291/64‑1 (MQ=255)                                                                                                                                        |                                                                                                                  CGCTCTTCCCCTGAAGAGTGGAAAGCATTAATTGCTGGTGGCTCGGAGAGAGGGTTGATTTATACCCAATTCGTTGCCCAAACCGCTATGTGCTGTGGGGAAGGGGGAATCAAAGCCTGGGATTATGTTCGAATGGGGTTTCTGTCACGGGTCGGTGTGCTCAACAAGTGGCTGACGGAAGAAGAGAGTTTATGGCTGCAATCGCGCGTCTATGTACGAGCTCACCATTACTATCACAGCTGGATGCAC  >  NC_000913/679739‑679987

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