Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,951,576	G→A	C61Y (TGT→TAT)	yecE →	DUF72 domain‑containing protein YecE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,951,576	0	G	A	100.0%	79.8 / NA	25	C61Y (TGT→TAT)	yecE	DUF72 domain‑containing protein YecE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/9);  total (16/9)

ACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCA  >  NC_000913/1951337‑1951809                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACAACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTc                                                                                                                                                                                                                                >  1:194100/1‑251 (MQ=255)                                         agagATCCTCAACGCGTTATGATTTCCATCGGTCTACCGCAATGGACGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCg                                                                                                                                                                                         >  1:49153/1‑250 (MQ=255)                                                           atgttttccagcGGTCTCCCGCAATGTTCGCACCATAAATGGGTGGGGTTGGGGATCCCCAGGCTTGAAGAGTATTCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAAtt                                                                                                                                                                      <  1:183794/240‑1 (MQ=255)                                                                                            ccTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCAt                                                                                                                                     >  2:281493/1‑251 (MQ=255)                                                                                                  tGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGc                                                                                                                                >  1:309574/1‑250 (MQ=255)                                                                                                     gTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTACCTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAAt                                                                                                                             >  1:265671/1‑250 (MQ=255)                                                                                                        cGGTTGGGGATCGCGAGCCTTGCAGAGTATGCCCGCCACTTTAAGTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTGGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGAATTATGGCATTGCGTTGGTTTAGTGACTGAATTTTTGACCCGCGTGGCACCGGTGTGTCCGCGCATTGGGCAGTACTg                                                                                                                         >  2:269929/1‑251 (MQ=255)                                                                                                          gtgggggcacACCAGCCTTGAAGATTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGc                                                                                                                       <  1:193911/242‑1 (MQ=255)                                                                                                                 aTCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGACAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAAct                                                                                                                >  1:170405/1‑251 (MQ=255)                                                                                                                              aaGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTc                                                                                                   >  1:392598/1‑251 (MQ=255)                                                                                                                                    tgncccGCCACGTTGACTGCGTGGAGGGCAAGACCACGCTTTACGCCCTGCCGAAACCCGAGGGTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCGGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGAGGATTTAGTGACGGAATTTTTGTCCCGCATGTCACCGTAGGCTCCGCGCATGGGGCAATAGTGGCTGCAACTGCCTGCCACATTTGGccc                                                                                              <  2:379719/247‑1 (MQ=255)                                                                                                                                        ccGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTACCTGCCACATTCGGCCCACggg                                                                                         >  1:404841/1‑251 (MQ=255)                                                                                                                                         cgcctatttAACTGCGTGGAGGGCCACACCACGATTTACGCCCTGCCGAAACCCGAGTTTGTCCTGCGCTGGCGTGAGCTGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGa                                                                                        <  1:211649/245‑1 (MQ=255)                                                                                                                                            cACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGct                                                                                     >  2:307440/1‑251 (MQ=255)                                                                                                                                                     tGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCt                                                                            >  1:73888/1‑251 (MQ=255)                                                                                                                                                           gAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCa                                                                      <  1:307440/251‑1 (MQ=255)                                                                                                                                                           gAGGGCAACACCACGCTGTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGTCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTTACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCa                                                                      <  2:136694/251‑1 (MQ=255)                                                                                                                                                             gggCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCAtt                                                                    >  1:340290/1‑251 (MQ=255)                                                                                                                                                               gCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCAGTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCAtttt                                                                  <  2:47002/251‑1 (MQ=255)                                                                                                                                                                     ccacGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGAGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGAATTTTCTCGAt                                                            >  1:169517/1‑251 (MQ=255)                                                                                                                                                                      cacGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGAtt                                                           >  2:174260/1‑251 (MQ=255)                                                                                                                                                                       acGCTTTACGCCCTCCCGAAACCCGAGGTTGTAATGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGACTTTTTGACCCGCATGTCACCGTTGGCTCCGAGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtc                                                          >  1:128916/1‑251 (MQ=255)                                                                                                                                                                                     gCCGAAAGCCGAGGTTGTCATGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACTCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATGCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAt                                            <  1:161129/251‑1 (MQ=255)                                                                                                                                                                                             ccGAGGTTGGCCTTCGCTGGCGTGAGCAGACCACAGAGGACTTCCGCTTCTATTGTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCAGTGCGAGGATTTAGTGACTGAATTTTGGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCATATTCGGCCCACGGGAGCTGTCTGCGCTTTGGTATTTTCTCGATTCGCGTCCTGGCGAATTTAATTAt                                    <  2:170038/251‑1 (MQ=255)                                                                                                                                                                                                                                gATGACTTCCGCTTCTATTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCa  >  1:361736/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          ACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCACTTTAACTGCGTGGAGGGCAACACCACGCTTTACGCCCTGCCGAAACCCGAGGTTGTCCTGCGCTGGCGTGAGCAGACCACAGATGACTTCCGCTTCTGTTTTAAGTTTCCGGCGACCATTTCGCATCAGGCAGCATTACGGCATTGCGATGATTTAGTGACTGAATTTTTGACCCGCATGTCACCGTTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCA  >  NC_000913/1951337‑1951809

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