Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,245,600	G→A	*147* (TAG→TAA)	ycgY →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,245,600	0	G	A	100.0%	119.7 / NA	38	*147* (TAG→TAA)	ycgY	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/20);  total (18/20)

TAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTAAC  >  NC_000913/1245465‑1245732                                                                                                                                        |                                                                                                                                     taGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaa                                                                                                                                   <  2:336997/139‑1 (MQ=255)   gCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAaaaaaa                                                                                                                                 <  2:694987/139‑1 (MQ=255)    ccAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTaaaaaaaa                                                                                                                                >  2:972333/1‑139 (MQ=255)     cAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAAc                                                                                                                               <  1:824396/139‑1 (MQ=255)       aCGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                             >  1:246843/1‑139 (MQ=255)       aCGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                             <  1:972333/139‑1 (MQ=255)       aCGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCa                                                                                                                             >  1:359196/1‑139 (MQ=255)          tCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAg                                                                                                                          <  1:140561/139‑1 (MQ=255)            ggCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        <  1:85889/139‑1 (MQ=255)              cGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTAt                                                                                                                      >  1:1002471/1‑139 (MQ=255)                                               aCCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCAt                                                                                              <  2:973266/130‑1 (MQ=255)                                               aCCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCAt                                                                                              >  1:973266/1‑130 (MQ=255)                                                       aaaTGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGa                                                                             >  2:211923/1‑139 (MQ=255)                                                              gCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                                  >  1:496331/1‑111 (MQ=255)                                                              gCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACg                                                                                                  <  2:496331/111‑1 (MQ=255)                                                                 aTCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGa                                                                   >  1:83225/1‑139 (MQ=255)                                                                 aTCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGa                                                                   >  1:972212/1‑139 (MQ=255)                                                                  tCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        <  2:244609/85‑1 (MQ=255)                                                                  tCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        <  1:417045/85‑1 (MQ=255)                                                                  tCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        >  1:556521/1‑85 (MQ=255)                                                                  tCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        >  1:244609/1‑85 (MQ=255)                                                                  tCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        <  2:556521/85‑1 (MQ=255)                                                                  tCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        >  2:417045/1‑85 (MQ=255)                                                                       ggCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                             <  1:452750/139‑1 (MQ=255)                                                                         cAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgcgc                                                           <  2:246843/139‑1 (MQ=255)                                                                                            tttGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTtgt                                        <  1:684402/139‑1 (MQ=255)                                                                                                     tACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTggg                                            >  2:976778/1‑126 (MQ=255)                                                                                                     tACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTggg                                            <  1:976778/126‑1 (MQ=255)                                                                                                              ccAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGt                      <  2:83225/139‑1 (MQ=255)                                                                                                               cAGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGtt                     >  1:239709/1‑139 (MQ=255)                                                                                                                aGGACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTg                    <  2:239709/139‑1 (MQ=255)                                                                                                                 ggACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                             >  2:378617/1‑97 (MQ=255)                                                                                                                 ggACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAgc                                                             <  1:378617/97‑1 (MQ=255)                                                                                                                  gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        >  1:841537/1‑37 (MQ=255)                                                                                                                  gACTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCt                                                                                                                        <  2:841537/37‑1 (MQ=255)                                                                                                                   aCTTGAAAAACCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGcgt                 <  2:972212/139‑1 (MQ=255)                                                                                                                            aCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGc           >  2:512022/1‑136 (MQ=255)                                                                                                                            aCCGTCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGc           <  1:512022/136‑1 (MQ=255)                                                                                                                                tCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGaa                                                                                      >  2:343298/1‑57 (MQ=255)                                                                                                                                tCTGTGTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGaa                                                                                      <  1:343298/57‑1 (MQ=255)                                                                                                                                  tgtgTAAAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTAAc  >  1:527994/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     TAGCCAACGTCGGCGTAAATCGCCGGGTGATGAAAGTGAAACACTTACCCGTAGAAATGACGCTATCTTCGGCAACCACGTATGGCAAACATTTGCGCAATACTTCCCGCCAGGACTTGAAAAACCGTCTGTGTAGAAAAAACCATAGCTATGGTGTTATCTGCCGACAACGCCATCCGGAGAAGCCGAATGACGTTGTCGGAGTAAGCGCTGGTTAAGGTGTGGGTTGTGCCTCTTTGGTTGAGGGTTGCGTCGTTGCTGACTTAAC  >  NC_000913/1245465‑1245732

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