Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,377,088	A→G	Y86C (TAC→TGC)	ycjS →	putative NADH‑binding oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,377,088	0	A	G	100.0%	142.5 / NA	38	Y86C (TAC→TGC)	ycjS	putative NADH‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/19);  total (19/19)

TGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTA  >  NC_000913/1376954‑1377214                                                                                                                                       |                                                                                                                               tgACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGaa                                                                                                                            <  2:253546/139‑1 (MQ=255) tgACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGaa                                                                                                                            <  2:372092/139‑1 (MQ=255)       ccgccTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATAcc                                                                                                                      >  1:487065/1‑139 (MQ=255)            tttCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGAt                                                                                                                 <  1:92775/139‑1 (MQ=255)                           cgcTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  <  2:690718/139‑1 (MQ=255)                           cgcTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  <  2:441500/139‑1 (MQ=255)                                     aaaaTACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCatgtg                                                                                        <  2:1010252/139‑1 (MQ=255)                                                aTGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGaa                                                                                                       <  2:728222/113‑1 (MQ=255)                                                aTGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGaa                                                                                                       >  1:728222/1‑113 (MQ=255)                                                 tGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGaaaaa                                                                            >  2:1057089/1‑139 (MQ=255)                                                   cATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAAcc                                                                          >  1:1143112/1‑139 (MQ=255)                                                    aTCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGaa                                                                                                                            >  1:810435/1‑88 (MQ=255)                                                    aTCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGaa                                                                                                                            <  2:810435/88‑1 (MQ=255)                                                         tgtgGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAAccgcccgc                                                                    >  1:527392/1‑139 (MQ=255)                                                            gggACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCAt                                                                 <  2:487065/139‑1 (MQ=255)                                                            gggACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCAt                                                                 <  2:1005117/139‑1 (MQ=255)                                                                  aTCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGcc                                                           >  2:143552/1‑139 (MQ=255)                                                                   tCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  >  2:94200/1‑99 (MQ=255)                                                                   tCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  <  1:94200/99‑1 (MQ=255)                                                                   tCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCa                                                          >  2:1010440/1‑139 (MQ=255)                                                                   tCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCa                                                          >  2:214068/1‑139 (MQ=255)                                                                      gCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGaa                                                       <  1:143552/139‑1 (MQ=255)                                                                          gCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCg                                                                                                    >  2:846082/1‑90 (MQ=255)                                                                          gCCATGCTGCTGGAGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCg                                                                                                    <  1:846082/90‑1 (MQ=255)                                                                              tgctgcTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  <  2:47028/88‑1 (MQ=255)                                                                              tgctgcTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  >  1:47028/1‑88 (MQ=255)                                                                               gctgctGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCgg                                              >  1:597747/1‑139 (MQ=255)                                                                                ctgctgGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCggg                                             <  1:1010440/139‑1 (MQ=255)                                                                                    tgGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  <  2:795836/82‑1 (MQ=255)                                                                                    tgGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCggg                                                                                                  >  1:795836/1‑82 (MQ=255)                                                                                                         ttAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATgg                                                                                                               >  1:517539/1‑48 (MQ=255)                                                                                                         ttAGCGTCTGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATgg                                                                                                               <  2:517539/48‑1 (MQ=255)                                                                                                                 tGCTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGTGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAgg            >  2:231516/1‑139 (MQ=255)                                                                                                                   cTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCAt                                                                 >  1:159873/1‑84 (MQ=255)                                                                                                                   cTCACCTAACCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCAt                                                                 <  2:159873/84‑1 (MQ=255)                                                                                                                           aCCGTTTTCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTa  >  1:512181/1‑139 (MQ=255)                                                                                                                                tttCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGTAAAACCGCCAGCCATGACGCCAGAACAgg                                                   <  2:786450/85‑1 (MQ=255)                                                                                                                                tttCATTGCGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAgg                                                   >  1:786450/1‑85 (MQ=255)                                                                                                                                       |                                                                                                                               TGACAGCCGCCTTTCCCAGGCGCAGGCGCTGGCAGAAAAATACGGGAATGCATCCGTGTGGGACGATCCGCAGGCCATGCTGCTGGCGGTGAAACCTGATGTGGTTAGCGTCTGCTCACCTAACCGTTTTCATTACGAACATACCCTGATGGCACTGGAAGCGGGCTGCCATGTGATGTGCGAAAAACCGCCCGCCATGACGCCAGAACAGGCGCGGGAAATGTGCGATACCGCGCGCAAACTGGGCAAGGTGCTGGCCTA  >  NC_000913/1376954‑1377214

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