Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	852,261	A→G	Y113H (TAC→CAC)	opgE ←	OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	852,261	0	A	G	100.0%	175.4 / NA	48	Y113H (TAC→CAC)	opgE	OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/29);  total (19/29)

CGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAG  >  NC_000913/852145‑852390                                                                                                                     |                                                                                                                                  cGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTc                                                                                                             <  1:1467084/139‑1 (MQ=255) cGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTGTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAACCGACCAGGTGGCTGTACCAGGGAAAAATCgac                                                                                                             >  2:1299283/1‑137 (MQ=255)               ccAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  2:1478589/139‑1 (MQ=255)               ccAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  1:256397/139‑1 (MQ=255)               ccAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  1:433085/139‑1 (MQ=255)               ccAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  1:462147/139‑1 (MQ=255)                 aGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               >  2:392250/1‑137 (MQ=255)                 aGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  1:392250/137‑1 (MQ=255)                  gACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGActc                                                                                            >  1:40905/1‑139 (MQ=255)                     acaGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTc                                                                                         <  2:437577/139‑1 (MQ=255)                     acaGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTc                                                                                         <  2:952581/139‑1 (MQ=255)                     acaGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTc                                                                                         <  1:905054/139‑1 (MQ=255)                     acaGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTc                                                                                         >  1:591549/1‑139 (MQ=255)                       aGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCa                                                                                       >  2:122447/1‑139 (MQ=255)                        ggCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTc                                                                                                             <  2:281838/116‑1 (MQ=255)                        ggCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTc                                                                                                             >  1:281838/1‑116 (MQ=255)                        ggCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCCCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCat                                                                                      >  2:731207/1‑139 (MQ=255)                          cTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               >  1:508138/1‑128 (MQ=255)                          cTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  2:508138/128‑1 (MQ=255)                               gAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATtggtg                                                                               <  2:519460/139‑1 (MQ=255)                                             gcgcTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGc                                                                         >  2:60254/1‑131 (MQ=255)                                             gcgcTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGc                                                                         <  1:60254/131‑1 (MQ=255)                                              cgcTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGG‑AAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCa                                                               >  2:1169779/1‑139 (MQ=255)                                                    ttCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTa                                                          <  2:1160467/139‑1 (MQ=255)                                                                ttttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGG‑AAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  1:1169779/139‑1 (MQ=255)                                                                 tttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  2:40905/139‑1 (MQ=255)                                                                 tttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  2:591549/139‑1 (MQ=255)                                                                 tttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  2:815749/139‑1 (MQ=255)                                                                 tttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  1:1299283/139‑1 (MQ=255)                                                                 tttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  1:122447/139‑1 (MQ=255)                                                                 tttCCCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  1:731207/139‑1 (MQ=255)                                                                     gCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               >  2:1139292/1‑85 (MQ=255)                                                                     gCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGa                                                                                               <  1:1139292/85‑1 (MQ=255)                                                                                 ttACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTc                                                                                                 <  2:1340912/71‑1 (MQ=255)                                                                                 ttACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTc                                                                                                 >  1:1340912/1‑71 (MQ=255)                                                                                    cGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAGTGCCAGCa                                                               >  2:972004/1‑102 (MQ=255)                                                                                    cGCCGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCa                                                               <  1:972004/102‑1 (MQ=255)                                                                                       cGAGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAgg                       >  1:992818/1‑139 (MQ=255)                                                                                         aGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTAcc                                                                                                                           >  2:152112/1‑37 (MQ=38)                                                                                         aGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTAcc                                                                                                                           <  1:152112/37‑1 (MQ=38)                                                                                         aGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             <  1:585571/115‑1 (MQ=255)                                                                                         aGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcg                                             >  2:585571/1‑115 (MQ=255)                                                                                         aGCGCAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGa                     >  1:314497/1‑139 (MQ=255)                                                                                         aGCGCAAAAAT‑AATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCa                                                                  >  2:56450/1‑93 (MQ=255)                                                                                         aGCGCAAAAAT‑AATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCa                                                                  <  1:56450/93‑1 (MQ=255)                                                                                           cgcAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCata                                                                                     <  2:297810/73‑1 (MQ=255)                                                                                           cgcAAAAATAAATAAGCCGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCata                                                                                     >  1:297810/1‑73 (MQ=255)                                                                                                            cGACCAGGTGGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACCCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAg  <  2:314497/139‑1 (MQ=255)                                                                                                                     |                                                                                                                                  CGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAG  >  NC_000913/852145‑852390

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