Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,634,127	A→G	T111A (ACC→GCC)	ytjC →	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,634,127	0	A	G	100.0%	151.5 / NA	41	T111A (ACC→GCC)	ytjC	putative phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/18);  total (23/18)

CCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGG  >  NC_000913/4633996‑4634259                                                                                                                                    |                                                                                                                                     cccAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGAc                                                                                                                              <  2:1071005/140‑1 (MQ=255)      gCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCg                                                                                                                         >  1:991633/1‑140 (MQ=255)            ggCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTcc                                                                                                                   >  1:522011/1‑140 (MQ=255)                gtgACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGaa                                                                                                               >  1:749696/1‑140 (MQ=255)                           tttGATTCTCGCCTGCGTGAATTAAACATGGGGGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAAt                                                                                                    >  1:593369/1‑140 (MQ=255)                               aTTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAg                                                                                                >  1:827493/1‑140 (MQ=255)                                ttCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAgg                                                                                               <  2:522011/140‑1 (MQ=255)                                  ctcGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAg                                                                                             >  2:262139/1‑140 (MQ=255)                                    cGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCt                                                                                           >  1:393136/1‑140 (MQ=255)                                                   aaCATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCaa                                                                            <  2:130961/140‑1 (MQ=255)                                                        gggTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCgg                                                                       <  2:827493/140‑1 (MQ=255)                                                         ggTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGc                                                                      >  1:27794/1‑140 (MQ=255)                                                             tgCTGGAAAAAAGACATATCTATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTg                                                                  >  2:465636/1‑140 (MQ=255)                                                             tgCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTg                                                                  >  1:1128357/1‑140 (MQ=255)                                                             tgCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTg                                                                  >  1:867240/1‑140 (MQ=255)                                                                 ggAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAAt                                                              <  2:867240/140‑1 (MQ=255)                                                                      aaaGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGc                                                         <  2:570724/140‑1 (MQ=255)                                                                            atatCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCAGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGAc                                                   >  1:793738/1‑140 (MQ=255)                                                                                      tctGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAgg                                         <  2:393136/140‑1 (MQ=255)                                                                                          aCCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAg                                     >  1:577980/1‑140 (MQ=255)                                                                                             gaagaagaaGAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCg                                  >  1:209522/1‑140 (MQ=255)                                                                                              aagaagaagAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGa                                 >  2:775820/1‑140 (MQ=255)                                                                                              aagaagaagAGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGa                                 >  1:315305/1‑140 (MQ=255)                                                                                               agaagaagaGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGAc                                <  2:690550/140‑1 (MQ=255)                                                                                                   gaagaGAACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCgct                            <  1:465636/140‑1 (MQ=255)                                                                                                       agaACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGctgctg                        <  2:577980/140‑1 (MQ=255)                                                                                                        gaACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTgg                       <  1:775820/140‑1 (MQ=255)                                                                                                        gaACTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTgg                       <  2:749696/140‑1 (MQ=255)                                                                                                         aaCTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGt                      >  1:612514/1‑140 (MQ=255)                                                                                                         aaCTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGt                      <  2:27794/140‑1 (MQ=255)                                                                                                           cTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGa                    >  1:157511/1‑140 (MQ=255)                                                                                                           cTGGCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGa                    >  1:287819/1‑140 (MQ=255)                                                                                                           cTGGCGTCGGCAGCTGGTCAATAGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGa                    >  1:576399/1‑140 (MQ=255)                                                                                                            tGGCGTCGGCAGCTGGTCAATAGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAg                   <  2:576399/140‑1 (MQ=255)                                                                                                             ggCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGt                  <  2:612514/140‑1 (MQ=255)                                                                                                              gCGTCGGCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTc                 >  1:467603/1‑140 (MQ=255)                                                                                                                   ggCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGt            >  1:817906/1‑140 (MQ=255)                                                                                                                    gCAGCTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTa           <  2:655688/140‑1 (MQ=255)                                                                                                                        cTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCa                <  2:39323/131‑1 (MQ=255)                                                                                                                        cTGGTCAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCa                >  1:39323/1‑131 (MQ=255)                                                                                                                             cAATGGCGCCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTgg  <  2:315305/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     CCCAGGCCTGCGGCTGTGACATCATCTTTGATTCTCGCCTGCGTGAATTAAACATGGGTGTGCTGGAAAAAAGACATATCGATTCTCTGACCGAAGAAGAAGAGAACTGGCGTCGGCAGCTGGTCAATGGCACCGTTGACGGGCGTATTCCTGAAGGCGAGTCAATGCAGGAGCTCAGCGATCGCGTCAATGCGGCACTGGAATCTTGCCGGGACTTACCGCAGGGAAGCCGACCGCTGCTGGTGAGTCATGGTATTGCACTGG  >  NC_000913/4633996‑4634259

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