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%	264.7 / NA	72	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 (36/36);  total (36/36)

TTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGAA  >  NC_000913/3320212‑3320466                                                                                                                          |                                                                                                                                      ttGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGc                                                                                                                      >  2:721703/1‑139 (MQ=255)         tCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACa                                                                                                              <  1:854048/139‑1 (MQ=255)              gCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGAcc                                                                                                         >  1:1002176/1‑139 (MQ=255)              gCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGAcc                                                                                                         >  2:1004857/1‑139 (MQ=255)                 aTCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAg                                                                                                      <  2:201691/139‑1 (MQ=255)                   ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                                    >  1:436586/1‑139 (MQ=255)                                         gATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               <  1:26809/122‑1 (MQ=255)                                         gATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               >  2:26809/1‑122 (MQ=255)                                         gATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                   <  2:10775/134‑1 (MQ=255)                                         gATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                   >  1:10775/1‑134 (MQ=255)                                              aCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTc                                                                                                                         <  1:415324/91‑1 (MQ=255)                                              aCATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTc                                                                                                                         >  2:415324/1‑91 (MQ=255)                                                 tCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg                                                                      >  1:538847/1‑139 (MQ=255)                                                   gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                                    >  1:208883/1‑139 (MQ=255)                                                   gggACATGATATGCCTACTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                                    <  2:1188674/139‑1 (MQ=255)                                                     gACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                  <  1:721703/139‑1 (MQ=255)                                                          gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCaa                                                             >  1:129614/1‑139 (MQ=255)                                                               gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               <  2:946205/100‑1 (MQ=255)                                                               gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               >  1:946205/1‑100 (MQ=255)                                                               gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                        <  2:436586/139‑1 (MQ=255)                                                               gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                        <  2:538847/139‑1 (MQ=255)                                                               gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                        <  2:1002176/139‑1 (MQ=255)                                                               gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                        <  1:45607/139‑1 (MQ=255)                                                                 ctgctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                      <  2:781110/139‑1 (MQ=255)                                                                  tgcttgctAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                        >  2:782836/1‑136 (MQ=255)                                                                  tgctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                        <  1:782836/136‑1 (MQ=255)                                                                   gctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCtgt                                                    >  2:1001796/1‑139 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               >  2:1269095/1‑95 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               <  1:1269095/95‑1 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGa                                                                                             <  2:741934/97‑1 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGa                                                                                             >  1:741934/1‑97 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                   <  1:9543/107‑1 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaata                                                                                   >  2:9543/1‑107 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg                                                                      >  2:38827/1‑120 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg                                                                      <  1:38827/120‑1 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCtgtg                                                   >  1:593382/1‑139 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCtgt                                                    <  2:292264/138‑1 (MQ=255)                                                                    ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCtgt                                                    >  1:292264/1‑138 (MQ=255)                                                                      ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAgg                                                                                                                                 <  1:1091818/59‑1 (MQ=255)                                                                      ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAgg                                                                                                                                 >  2:1091818/1‑59 (MQ=255)                                                                         tAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                  >  1:341270/1‑119 (MQ=255)                                                                         tAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                  <  2:341270/119‑1 (MQ=255)                                                                          aGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                  >  2:202862/1‑118 (MQ=255)                                                                          aGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCt                                                                  <  1:202862/118‑1 (MQ=255)                                                                           gCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               >  2:512832/1‑88 (MQ=255)                                                                           gCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               <  1:512832/88‑1 (MQ=255)                                                                            ccTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               <  1:451223/87‑1 (MQ=255)                                                                            ccTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               >  2:451223/1‑87 (MQ=255)                                                                                 cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTc                                                                                                                                     >  1:109993/1‑44 (MQ=255)                                                                                 cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTc                                                                                                                                     <  2:109993/44‑1 (MQ=255)                                                                                         aTGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAg                              >  2:1132920/1‑139 (MQ=255)                                                                                          tGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGc                             <  2:208883/139‑1 (MQ=255)                                                                                           gAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCa                            <  1:1004857/139‑1 (MQ=255)                                                                                            aaaaTAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAg                           <  2:129614/139‑1 (MQ=255)                                                                                               aTAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                      >  2:660054/1‑109 (MQ=255)                                                                                               aTAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                      <  1:660054/109‑1 (MQ=255)                                                                                                    cGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGt                                                                                                                          <  2:1305623/36‑1 (MQ=255)                                                                                                    cGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGt                                                                                                                          >  1:1305623/1‑36 (MQ=255)                                                                                                        cttcttcGTAGGGTCACGCTCCAGGCCGTGGTCgg                                                                                                                       >  2:16509/1‑35 (MQ=255)                                                                                                        cttcttcGTAGGGTCACGCTCCAGGCCGTGGTCgg                                                                                                                       <  1:16509/35‑1 (MQ=255)                                                                                                        cttcttcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCg               >  1:1002536/1‑139 (MQ=255)                                                                                                                 aGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                                    >  2:733850/1‑77 (MQ=255)                                                                                                                 aGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                                    <  1:733850/77‑1 (MQ=255)                                                                                                                 aGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGt      <  1:1132920/139‑1 (MQ=255)                                                                                                                   ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                      <  1:1201544/89‑1 (MQ=255)                                                                                                                   ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCt                                                      >  2:1201544/1‑89 (MQ=255)                                                                                                                   ggTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTg    >  2:404349/1‑139 (MQ=255)                                                                                                                     tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               <  2:816600/46‑1 (MQ=255)                                                                                                                     tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                               >  1:816600/1‑46 (MQ=255)                                                                                                                     tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGaa  >  2:570872/1‑139 (MQ=255)                                                                                                                      cACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTa                                                 <  2:485287/91‑1 (MQ=255)                                                                                                                      cACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTa                                                 >  1:485287/1‑91 (MQ=255)                                                                                                                          |                                                                                                                                      TTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTGCGGTTGAA  >  NC_000913/3320212‑3320466

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