Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,586,281	C→T	G69D (GGC→GAC)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,586,281	0	C	T	100.0%	187.9 / NA	58	G69D (GGC→GAC)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/28);  total (30/28)

GGCTGAGTGAGTCGCCCCGCAGCCTCAAGCTCGTGGTCTCCCCAGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTT  >  NC_000913/1586152‑1586408                                                                                                                                  |                                                                                                                                ggCTGAGTGAGTCGCCCCGCAGCCTCAAGCTCGTGGTCTCCCCAGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTcac                                                                                                                        >  1:31436/1‑139 (MQ=255)                  cGCAGCCTCAAGCTCGTGGTCTCCCCAGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       >  2:491012/1‑139 (MQ=255)                  cGCAGCCTCAAGCTCGTGGTCTCCCCAGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       <  1:1431121/139‑1 (MQ=255)                                              gtaagtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACtgt                                                                                                   >  1:1193911/1‑115 (MQ=255)                                              gtaagtaaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACtgt                                                                                                   <  2:1193911/115‑1 (MQ=255)                                                   taaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGa                                                                      >  2:711100/1‑139 (MQ=255)                                                   taaTGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGa                                                                      >  1:297168/1‑139 (MQ=255)                                                       gATTGAACCGTATTCTCAGTAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGtt                                                                                                                             >  1:1041660/1‑80 (MQ=255)                                                       gATTGAACCGTATTCTCAGTAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGtt                                                                                                                             <  2:1041660/80‑1 (MQ=255)                                                        aTTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGAc                                                                                                          >  2:1502658/1‑98 (MQ=255)                                                        aTTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGAc                                                                                                          <  1:1502658/98‑1 (MQ=255)                                                        aTTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       >  2:720719/1‑101 (MQ=255)                                                        aTTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       <  1:720719/101‑1 (MQ=255)                                                        aTTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                 >  1:1386619/1‑139 (MQ=255)                                                        aTTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                 >  2:1012213/1‑139 (MQ=255)                                                             aCCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGc                                                            >  1:8542/1‑139 (MQ=255)                                                                gTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGCCAAAGCCct                                                         <  1:1012213/139‑1 (MQ=255)                                                                        aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaa                                                               <  2:218934/125‑1 (MQ=255)                                                                        aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaa                                                               >  1:218934/1‑125 (MQ=255)                                                                         gCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTAt                                                >  2:16014/1‑139 (MQ=255)                                                                          cAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGc                                                                                >  2:569588/1‑106 (MQ=255)                                                                          cAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGc                                                                                <  1:569588/106‑1 (MQ=255)                                                                             agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt                                            >  2:673379/1‑139 (MQ=255)                                                                                    ggCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATg                                                                                                            <  1:550757/68‑1 (MQ=255)                                                                                    ggCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATg                                                                                                            >  2:550757/1‑68 (MQ=255)                                                                                    ggCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAAc                                     <  1:785881/139‑1 (MQ=255)                                                                                         ttACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGc                                >  1:949420/1‑139 (MQ=255)                                                                                               gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAAc                          >  1:606718/1‑139 (MQ=255)                                                                                                  tATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATc                       <  1:16014/139‑1 (MQ=255)                                                                                                  tATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATc                       <  2:318419/139‑1 (MQ=255)                                                                                                  tATCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATc                       <  2:949420/139‑1 (MQ=255)                                                                                                   aTCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACt                                                                                                     >  1:933917/1‑60 (MQ=255)                                                                                                   aTCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACt                                                                                                     <  2:933917/60‑1 (MQ=255)                                                                                                   aTCCGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATcc                      >  2:566301/1‑139 (MQ=255)                                                                                                     ccGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       <  2:1042815/56‑1 (MQ=255)                                                                                                     ccGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       >  1:1042815/1‑56 (MQ=255)                                                                                                      cGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGAc                                                                                                          <  2:199426/52‑1 (MQ=255)                                                                                                      cGTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGAc                                                                                                          >  1:199426/1‑52 (MQ=255)                                                                                                       gTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGc                  <  2:606718/139‑1 (MQ=255)                                                                                                       gTGACTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGc                  <  2:297168/139‑1 (MQ=255)                                                                                                          aCTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTGGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCATCTGGCTTATTCATGTCAAACATGGCAATACAATCCTGTCGta               >  2:1438791/1‑139 (MQ=255)                                                                                                          aCTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGa                                                                                                               >  2:1388074/1‑43 (MQ=255)                                                                                                          aCTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGa                                                                                                               <  1:1388074/43‑1 (MQ=255)                                                                                                          aCTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGa                                                                      >  2:1262184/1‑84 (MQ=255)                                                                                                          aCTTCCCAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGa                                                                      <  1:1262184/84‑1 (MQ=255)                                                                                                              cccAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       <  2:1422839/47‑1 (MQ=255)                                                                                                              cccAGGCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCa                                                                                                       >  1:1422839/1‑47 (MQ=255)                                                                                                                  ggtgATTGCTTTTGCGTCGTTTTCACAAAAGTCGAATTACGCACTTTGCCTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATc       <  1:1438791/136‑1 (MQ=255)                                                                                                                  ggCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATc       <  2:31436/139‑1 (MQ=255)                                                                                                                  ggCGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATc       <  1:673379/139‑1 (MQ=255)                                                                                                                    cGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCAc                                                                                                      <  1:795972/42‑1 (MQ=255)                                                                                                                    cGATTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCAc                                                                                                      >  2:795972/1‑42 (MQ=255)                                                                                                                      aTTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaaa                                                              <  2:1392211/80‑1 (MQ=255)                                                                                                                      aTTGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaaa                                                              >  1:1392211/1‑80 (MQ=255)                                                                                                                       ttGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCtt  >  2:28779/1‑139 (MQ=255)                                                                                                                       ttGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCtt  <  2:8542/139‑1 (MQ=255)                                                                                                                        tGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCct                                                         >  2:1239852/1‑83 (MQ=255)                                                                                                                        tGCTTTTGCGTCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCct                                                         <  1:1239852/83‑1 (MQ=255)                                                                                                                                  |                                                                                                                                GGCTGAGTGAGTCGCCCCGCAGCCTCAAGCTCGTGGTCTCCCCAGGTAAGTAATGATTGAACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTT  >  NC_000913/1586152‑1586408

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