Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,320,333	T→C	E427E (GAA→GAG)	yhbX ←	putative EptAB family phosphoethanolamine transferase, inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,320,333	0	T	C	100.0%	213.6 / NA	59	E427E (GAA→GAG)	yhbX	putative EptAB family phosphoethanolamine transferase, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/28);  total (31/28)

GATGTTGTTTTCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTG  >  NC_000913/3320198‑3320464                                                                                                                                        |                                                                                                                                    gATGTTGTTTTCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTcc                                                                                                                                  <  1:878359/139‑1 (MQ=255)                tgtgCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCTTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                                   <  1:922017/139‑1 (MQ=255)                            gCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGAcc                                                                                                       >  2:1068798/1‑139 (MQ=255)                                 ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                                  >  1:182538/1‑139 (MQ=255)                                 ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                                  >  1:950630/1‑139 (MQ=255)                                 ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                                  >  1:787430/1‑139 (MQ=255)                                 ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                                  >  1:584265/1‑139 (MQ=255)                                  cccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcg                                                                                                 <  1:1016023/139‑1 (MQ=255)                                  cccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcg                                                                                                 >  2:1016023/1‑139 (MQ=255)                                     aaGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcgg                                                                                              >  1:130013/1‑139 (MQ=255)                                                tataCCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaa                                                                                   <  2:1057721/139‑1 (MQ=255)                                                  taCCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                 <  2:1163509/139‑1 (MQ=255)                                                       gATAAACATCGGGACATGATATGTCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTc                                                                                                                       <  2:10418/96‑1 (MQ=255)                                                       gATAAACATCGGGACATGATATGTCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTc                                                                                                                       >  1:10418/1‑96 (MQ=255)                                                       gATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcgc                                                                                                <  1:641461/119‑1 (MQ=255)                                                       gATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcgc                                                                                                >  2:641461/1‑119 (MQ=255)                                                            aCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAAc                                                                       >  2:37926/1‑139 (MQ=255)                                                            aCATCGGGACATGATATGCATGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAAc                                                                       >  2:812145/1‑139 (MQ=255)                                                                cGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGa                                                                   >  2:860789/1‑139 (MQ=255)                                                                 gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                                  <  2:245726/139‑1 (MQ=255)                                                                  ggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAcc                                                                 <  2:584265/139‑1 (MQ=255)                                                                  ggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGACCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAcc                                                                 >  1:874666/1‑139 (MQ=255)                                                                   gACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                <  2:130013/139‑1 (MQ=255)                                                                        gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                <  1:924869/134‑1 (MQ=255)                                                                        gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                >  2:924869/1‑134 (MQ=255)                                                                            tGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTa                                                       <  2:950630/139‑1 (MQ=255)                                                                             gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             >  2:315766/1‑100 (MQ=255)                                                                             gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             <  1:315766/100‑1 (MQ=255)                                                                              ccTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATc                                                     >  2:257933/1‑139 (MQ=255)                                                                                  ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                >  1:1040005/1‑124 (MQ=255)                                                                                  ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                <  2:1040005/124‑1 (MQ=255)                                                                                   tgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             >  1:80143/1‑94 (MQ=255)                                                                                   tgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             <  2:80143/94‑1 (MQ=255)                                                                                    ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             >  1:35195/1‑93 (MQ=255)                                                                                    ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             <  2:35195/93‑1 (MQ=255)                                                                                    ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTa                                               >  1:949985/1‑139 (MQ=255)                                                                                        aGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTaatga                                           <  1:37926/139‑1 (MQ=255)                                                                                          ccTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTaatgaat                                         >  1:1174926/1‑139 (MQ=255)                                                                                             cccTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                      <  2:874666/139‑1 (MQ=255)                                                                                               cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGt                                    >  1:295380/1‑139 (MQ=255)                                                                                                     ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCAt                              <  2:949985/139‑1 (MQ=255)                                                                                                     ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCAt                              <  2:404072/139‑1 (MQ=255)                                                                                                      cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             <  2:989364/75‑1 (MQ=255)                                                                                                      cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                             >  1:989364/1‑75 (MQ=255)                                                                                                       aTGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAAc                                                                       >  2:395236/1‑96 (MQ=255)                                                                                                       aTGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAAc                                                                       <  1:395236/96‑1 (MQ=255)                                                                                                         gAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCa                          <  2:295380/139‑1 (MQ=255)                                                                                                         gAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCa                          <  2:182538/139‑1 (MQ=255)                                                                                                          aaaaTAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAg                         <  2:173643/139‑1 (MQ=255)                                                                                                              tAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                    >  1:459102/1‑108 (MQ=255)                                                                                                              tAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                    <  2:459102/108‑1 (MQ=255)                                                                                                                             gTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCg      >  1:833631/1‑139 (MQ=255)                                                                                                                               aGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGt    >  2:722006/1‑139 (MQ=255)                                                                                                                                 ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                      >  1:1197672/1‑87 (MQ=255)                                                                                                                                 ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                      <  2:1197672/87‑1 (MQ=255)                                                                                                                                 ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                      >  2:626032/1‑103 (MQ=255)                                                                                                                                 ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                      <  1:626032/103‑1 (MQ=255)                                                                                                                                 ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGTGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTg  >  2:764822/1‑139 (MQ=255)                                                                                                                                 ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTg  >  1:707042/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    GATGTTGTTTTCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTG  >  NC_000913/3320198‑3320464

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