Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	851,094	G→A	P502S (CCG→TCG)	opgE ←	OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	851,094	0	G	A	100.0%	209.9 / NA	61	P502S (CCG→TCG)	opgE	OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/36);  total (25/36)

GGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGG  >  NC_000913/850999‑851189                                                                                                |                                                                                                ggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCg                                                                                             <  2:346453/100‑1 (MQ=255)   cGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGc                                                                                           <  1:840843/100‑1 (MQ=255)    ggATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCt                                                                                          <  2:1148631/100‑1 (MQ=255)     ggtttttGTTTTTAAGGTTTCGTGTCGAAATTATCGGTTCTGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGc                                                                                           <  2:994890/96‑1 (MQ=255)         tttCTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTttctt                                                                                     <  2:120082/100‑1 (MQ=255)            gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                  <  2:225826/100‑1 (MQ=255)            gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                  <  1:1560054/100‑1 (MQ=255)            gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                  <  1:263750/100‑1 (MQ=255)             tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                 <  2:548888/100‑1 (MQ=255)             tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                 <  2:1322061/100‑1 (MQ=255)               tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGa                                                                               >  2:1278425/1‑100 (MQ=255)                ttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGat                                                                              <  1:790805/100‑1 (MQ=255)                            gTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAt                                                                  <  1:250759/100‑1 (MQ=255)                               gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                               <  2:466862/100‑1 (MQ=255)                               gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                               >  1:1501690/1‑100 (MQ=255)                               gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                               >  1:15804/1‑100 (MQ=255)                               gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCCGCTTTCTTCTTAGATATAAAAGGGTATTTg                                                               >  2:1457487/1‑100 (MQ=255)                                 aaaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGAt                                                             <  1:482153/100‑1 (MQ=255)                                  aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                            <  2:642655/100‑1 (MQ=255)                                     atCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTa                                                         <  2:898487/100‑1 (MQ=255)                                        ggTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAtt                                                      >  1:1086270/1‑100 (MQ=255)                                                 aTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCt                                             <  2:283325/100‑1 (MQ=255)                                                  tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                            >  1:635125/1‑100 (MQ=255)                                                  tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                            >  1:512041/1‑100 (MQ=255)                                                   ggTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCttt                                           <  1:1327683/100‑1 (MQ=255)                                                   ggTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCttt                                           <  1:1328825/100‑1 (MQ=255)                                                   ggTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                            <  2:1307070/99‑1 (MQ=255)                                                    gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                          >  2:272198/1‑100 (MQ=255)                                                    gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                          >  1:1373479/1‑100 (MQ=255)                                                    gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                          >  2:827393/1‑100 (MQ=255)                                                       gTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATc                                       <  1:658891/100‑1 (MQ=255)                                                        tAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcc                                      <  1:280115/100‑1 (MQ=255)                                                         aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATccc                                     >  2:1432994/1‑100 (MQ=255)                                                           tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATGCCCg                                   >  1:957787/1‑100 (MQ=255)                                                           tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                   >  1:610665/1‑100 (MQ=255)                                                           tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                   <  1:251372/100‑1 (MQ=255)                                                             cACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTc                                 >  2:870741/1‑100 (MQ=255)                                                             cACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTAGCTTTAATACCTGTc                                 >  2:1482913/1‑100 (MQ=255)                                                              aCCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTcc                                <  2:1127261/100‑1 (MQ=255)                                                                   cTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgg                           <  2:5743/100‑1 (MQ=255)                                                                   cTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgg                           <  1:318106/100‑1 (MQ=255)                                                                    tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                          <  2:797593/100‑1 (MQ=255)                                                                    tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                          <  1:1175193/100‑1 (MQ=255)                                                                    tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                          <  2:913496/100‑1 (MQ=255)                                                                    tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                          <  2:1298698/100‑1 (MQ=255)                                                                     gTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                         <  2:1571256/100‑1 (MQ=255)                                                                      tAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                        <  1:839612/100‑1 (MQ=255)                                                                       aaCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATACACGGGTATTTGATGTTAATTTCCTTCCCTTTAATCCCCATCCACTGGgaca                       >  2:510861/1‑98 (MQ=255)                                                                        aCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaa                      >  2:1391390/1‑100 (MQ=255)                                                                          tGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaa                     >  2:914781/1‑99 (MQ=255)                                                                            gAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAg                  <  1:960347/100‑1 (MQ=255)                                                                                 tGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaa             >  1:1388656/1‑100 (MQ=255)                                                                                  ggTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaag            <  2:535587/100‑1 (MQ=255)                                                                                  ggTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaag            <  1:140445/100‑1 (MQ=255)                                                                                       tGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTc       <  1:1380319/100‑1 (MQ=255)                                                                                       tGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTc       >  1:691896/1‑100 (MQ=255)                                                                                        gTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCa      >  2:1559071/1‑100 (MQ=255)                                                                                           tatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg   >  2:1380578/1‑100 (MQ=255)                                                                                            atCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg  >  1:600572/1‑100 (MQ=255)                                                                                            atCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg  >  2:388523/1‑100 (MQ=255)                                                                                            atCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg   >  2:1437862/1‑99 (MQ=255)                                                                                                |                                                                                                GGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGG  >  NC_000913/850999‑851189

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