Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,629,481	G→A	*64* (TAG→TAA)	yfgG →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,629,481	0	G	A	100.0%	84.8 / NA	28	*64* (TAG→TAA)	yfgG	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/19);  total (9/19)

TATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACT  >  NC_000913/2629346‑2629602                                                                                                                                        |                                                                                                                          tATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACgg                                                                                                                        <  1:232383/139‑1 (MQ=255) tATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACgg                                                                                                                        >  2:710517/1‑139 (MQ=255)              atcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCAc                                                                                                           <  1:162906/139‑1 (MQ=255)              atcaGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCAc                                                                                                           <  2:847359/139‑1 (MQ=255)                   cTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                      <  1:38088/139‑1 (MQ=255)                      tATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAg                                                                                                   <  2:96191/139‑1 (MQ=255)                             ttcttTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATaa                                                                                            <  2:201603/139‑1 (MQ=255)                                tttGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACgg                                                                                         >  2:577726/1‑139 (MQ=255)                                  tGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCa                                                                                       >  1:168431/1‑139 (MQ=255)                                  tGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCa                                                                                       <  1:710517/139‑1 (MQ=255)                                    gCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACt                                                                                     >  2:585582/1‑139 (MQ=255)                                                  tCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGc                                                                       >  2:818776/1‑139 (MQ=255)                                                     tccgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCttt                                                                    <  1:818776/139‑1 (MQ=255)                                                       cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                  <  1:710988/139‑1 (MQ=255)                                                       cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGATACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                  <  1:577726/139‑1 (MQ=255)                                                         tcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCt                                                                >  1:518926/1‑139 (MQ=255)                                                             tGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGttt                                                            <  2:249569/139‑1 (MQ=255)                                                                              agagCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGc                                           <  2:162130/139‑1 (MQ=255)                                                                               gagCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCa                                          <  2:518926/139‑1 (MQ=255)                                                                                            gCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                  <  2:630060/102‑1 (MQ=255)                                                                                            gCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                  >  1:630060/1‑102 (MQ=255)                                                                                                       cacaCTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTa                  <  2:495485/139‑1 (MQ=255)                                                                                                              tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCttt           >  2:569264/1‑139 (MQ=255)                                                                                                              tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCttt           <  2:406105/139‑1 (MQ=255)                                                                                                                       tCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact  <  2:463406/139‑1 (MQ=255)                                                                                                                                   cGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                  <  2:232722/63‑1 (MQ=255)                                                                                                                                   cGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                  >  1:232722/1‑63 (MQ=255)                                                                                                                                     ttAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTa                  <  1:910741/109‑1 (MQ=255)                                                                                                                                        |                                                                                                                          TATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACT  >  NC_000913/2629346‑2629602

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