Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	851,235	A→G	S455P (TCC→CCC)	opgE ←	phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	851,235	0	A	G	100.0%	167.0 / NA	47	S455P (TCC→CCC)	opgE	phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/19);  total (28/19)

GCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGAAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAGCTGCTGCCGCTGTT  >  NC_000913/851099‑851370                                                                                                                                         |                                                                                                                                        gCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAAt                                                                                                                                      >  2:277116/1‑140 (MQ=255)                   aaCGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTg                                                                                                                    >  2:450066/1‑140 (MQ=255)                        gTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAAt                                                                                                               >  2:80863/1‑140 (MQ=255)                        gTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGAGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAAt                                                                                                               >  1:54164/1‑140 (MQ=255)                           tttGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAAtttt                                                                                                            >  1:781073/1‑140 (MQ=255)                            ttGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTg                                                                                                           <  2:54164/140‑1 (MQ=255)                                tGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGa                                                                                                       <  2:793992/140‑1 (MQ=255)                                  ttAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCt                                                                                                         >  1:122840/1‑136 (MQ=255)                                  ttAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCt                                                                                                         <  2:122840/136‑1 (MQ=255)                                          cTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTatcat                                                                                             >  2:191334/1‑140 (MQ=255)                                           ttCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTatcatc                                                                                            <  1:155979/140‑1 (MQ=255)                                           ttCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTatcatc                                                                                            <  1:80863/140‑1 (MQ=255)                                           ttCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTatcatc                                                                                            <  1:696535/140‑1 (MQ=255)                                            tCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTAtcatca                                                                                           >  2:252889/1‑140 (MQ=255)                                                 ttAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGc                                                                                      <  2:608244/140‑1 (MQ=255)                                                         cGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGc                                                                                      <  1:283962/132‑1 (MQ=255)                                                         cGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGc                                                                                      >  2:283962/1‑132 (MQ=255)                                                          gTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTg                                                                             >  1:69465/1‑140 (MQ=255)                                                            ccACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCa                                                                           <  2:69465/140‑1 (MQ=255)                                                            ccACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCa                                                                           <  1:450066/140‑1 (MQ=255)                                                                 gggAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCt                                                                      >  1:407045/1‑140 (MQ=255)                                                                 gggAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCGGTAATTTTTCATCATGGGCAAGGAATTGCCGGTCt                                                                      >  1:728069/1‑140 (MQ=255)                                                                    agaAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTa                                                                   >  1:624707/1‑140 (MQ=255)                                                                       aaaaGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACcg                                                                <  1:323221/140‑1 (MQ=255)                                                                       aaaaGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACcg                                                                <  2:719556/140‑1 (MQ=255)                                                                        aaaGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACcgc                                                               <  2:6239/140‑1 (MQ=255)                                                                           gCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGct                                                            <  1:191334/140‑1 (MQ=255)                                                                             cTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGctct                                                          >  1:605167/1‑140 (MQ=255)                                                                                     tCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGc                                                  >  2:340813/1‑140 (MQ=255)                                                                                       aTTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCa                                                >  2:95107/1‑140 (MQ=255)                                                                                            cTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCg                                           >  2:306590/1‑140 (MQ=255)                                                                                             tGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCgg                                          >  2:41550/1‑139 (MQ=255)                                                                                              gAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTg                                         >  1:284475/1‑140 (MQ=255)                                                                                               aGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTgg                                        >  1:100060/1‑140 (MQ=255)                                                                                               aGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTgg                                        >  1:497609/1‑140 (MQ=255)                                                                                                gcggcgGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGt                                       >  2:288138/1‑140 (MQ=255)                                                                                                   gcggGCTTTAATCACACGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAg                                    <  1:258598/140‑1 (MQ=255)                                                                                                                 cacGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGc                      >  1:66269/1‑140 (MQ=255)                                                                                                                 cacGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGc                      >  1:380844/1‑140 (MQ=255)                                                                                                                 cacGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGc                      >  1:479393/1‑140 (MQ=255)                                                                                                                 cacGATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGc                      >  2:250269/1‑140 (MQ=255)                                                                                                                    gATGCGCTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCgag                   <  2:624707/140‑1 (MQ=255)                                                                                                                          cTTTATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAgctg             >  2:403979/1‑140 (MQ=255)                                                                                                                             tATCGTCGCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAGctgctg          <  2:671822/140‑1 (MQ=255)                                                                                                                                 gtcgCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAGCTGCTGCCGc      <  1:340813/140‑1 (MQ=255)                                                                                                                                   cgCTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAGCTGCTGCCGCTg    <  1:252889/140‑1 (MQ=255)                                                                                                                                     cTGGGAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAGCTGCTGCCGCTGtt  >  1:624072/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                        GCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTGAGCGGCGGGCTTTAATCACACGATGCGCTTTATCGTCGCTGGAAATGACCATAAAAGGCACCTGGAAATTTTGCTGATATTTATCATCATGGGCAAGGTATTGCACGTCTTTACCGCGCTCTTTAAAGGCCAGACCGTGGTCAGAAAAGTAAACCAGCGAGAAGCTGCTGCCGCTGTT  >  NC_000913/851099‑851370

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: ‑