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%	239.0 / NA	66	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/30);  total (36/30)

GTTGTTTTCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTG  >  NC_000913/3320201‑3320458                                                                                                                                     |                                                                                                                              gttgttTTCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAgg                                                                                                                         >  1:895786/1‑139 (MQ=255)                    tCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTTCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACa                                                                                                      >  2:477028/1‑139 (MQ=255)                         gCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGAcc                                                                                                 >  1:474117/1‑139 (MQ=255)                              ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                            >  1:6409/1‑139 (MQ=255)                              ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                            >  2:1060721/1‑139 (MQ=255)                              ccccAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgc                                                                                            >  1:629387/1‑139 (MQ=255)                                  aaGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGACGCgcgg                                                                                        >  1:857948/1‑139 (MQ=255)                                   aGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       <  1:708259/139‑1 (MQ=255)                                                     aTAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                             <  2:537501/99‑1 (MQ=255)                                                     aTAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                             >  1:537501/1‑99 (MQ=255)                                                           aTCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCt                                                               >  2:628055/1‑139 (MQ=255)                                                            tCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTg                                                              >  1:109775/1‑139 (MQ=255)                                                              gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGt                                                                                                                  >  2:666839/1‑85 (MQ=255)                                                              gggACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGt                                                                                                                  <  1:666839/85‑1 (MQ=255)                                                                   aTGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGc                                                       <  1:1095964/139‑1 (MQ=255)                                                                    tGATATGCCTGCTGGATAGCCTCCCTGCCTCTAGGAAAATAGACGTTCTTCTTCGTAGGGTCACGCCCCAGACCGTGGTCGGCAAAATACATGCCCGAGGCGCGGCGATATTTTAATAATTCAAAAACCTGACCTAGCa                                                      >  1:1090983/1‑139 (MQ=255)                                                                     gATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCaa                                                     >  2:6648/1‑139 (MQ=255)                                                                          gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       >  2:85074/1‑100 (MQ=255)                                                                          gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       <  1:85074/100‑1 (MQ=255)                                                                          gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                <  1:628055/139‑1 (MQ=255)                                                                          gCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                <  2:485544/139‑1 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCaaa                                                                                                           >  2:773207/1‑75 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCaaa                                                                                                           <  1:773207/75‑1 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       >  2:416618/1‑95 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       <  1:416618/95‑1 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                            <  2:1048786/122‑1 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGAc                                                            >  1:1048786/1‑122 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                <  2:965525/134‑1 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                >  1:965525/1‑134 (MQ=255)                                                                               ctgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCAGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCtgtg                                           >  2:338171/1‑139 (MQ=255)                                                                                tgGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATCCAAAAACCTGACCTAGCAAACTATCTgtgt                                          <  1:1060721/139‑1 (MQ=255)                                                                                 ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTa                                         <  2:895786/139‑1 (MQ=255)                                                                                 ggCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTa                                         <  1:6648/139‑1 (MQ=255)                                                                                  gCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTaa                                        >  1:944853/1‑139 (MQ=255)                                                                                    tAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCTGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       >  1:422377/1‑90 (MQ=255)                                                                                    tAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCTGTGGTCGGCAAAATACAAGACCGAGGCgcggc                                                                                       <  2:422377/90‑1 (MQ=255)                                                                                          cccTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                <  2:474117/139‑1 (MQ=255)                                                                                          cccTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGa                                <  2:944853/139‑1 (MQ=255)                                                                                           ccTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAg                               <  2:109775/139‑1 (MQ=255)                                                                                            cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                             >  2:518064/1‑60 (MQ=255)                                                                                            cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCa                                                                                                             <  1:518064/60‑1 (MQ=255)                                                                                            cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                <  2:421460/121‑1 (MQ=255)                                                                                            cTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTAt                                                >  1:421460/1‑121 (MQ=255)                                                                                              gCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTg                            >  2:339033/1‑139 (MQ=255)                                                                                                 tCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAg                                                        >  2:926806/1‑108 (MQ=255)                                                                                                 tCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAg                                                        <  1:926806/108‑1 (MQ=255)                                                                                                 tCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCa                         >  2:126050/1‑139 (MQ=255)                                                                                                  ccATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGAGATCTTTTAATAATTCAAAGACCTGACCTAGCAAACTATCTATGTAATGAATGGCGTTGTCAt                        <  2:1090983/139‑1 (MQ=255)                                                                                                   cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTa                                                 <  1:1111299/113‑1 (MQ=255)                                                                                                   cATGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTa                                                 >  2:1111299/1‑113 (MQ=255)                                                                                                    aTGAAAATAGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGACGCGCGGCGATCTCTTAATAATTCAAAAACCTGACCTAGCAAACTATCTACGTAATGAATGGAGTTGTCATAg                      >  2:301505/1‑139 (MQ=255)                                                                                                            aGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcg                                                                                           <  1:819412/62‑1 (MQ=255)                                                                                                            aGACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGgcg                                                                                           >  2:819412/1‑62 (MQ=255)                                                                                                             gACGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATcc             >  1:348246/1‑139 (MQ=255)                                                                                                              aCGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaat                                                                            <  2:886579/75‑1 (MQ=255)                                                                                                              aCGTTCTTCTTCGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaat                                                                            >  1:886579/1‑75 (MQ=255)                                                                                                                      cttcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTcc    >  1:631905/1‑139 (MQ=255)                                                                                                                        tcGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTg  <  1:126050/139‑1 (MQ=255)                                                                                                                         cGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTa                                                 >  1:105487/1‑91 (MQ=255)                                                                                                                         cGTAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTa                                                 <  2:105487/91‑1 (MQ=255)                                                                                                                           tAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       >  2:665495/1‑51 (MQ=255)                                                                                                                           tAGGGTCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCgcggc                                                                                       <  1:665495/51‑1 (MQ=255)                                                                                                                                tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGt                           <  2:586500/106‑1 (MQ=255)                                                                                                                                tCACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGt                           >  1:586500/1‑106 (MQ=255)                                                                                                                                 cACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaa                                                                             >  1:462315/1‑55 (MQ=255)                                                                                                                                 cACGCTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTtaa                                                                             <  2:462315/55‑1 (MQ=255)                                                                                                                                     |                                                                                                                              GTTGTTTTCCGTTGTGCGATCCACGCCATCCCCAAGAACGGGGCTATACCAGATAAACATCGGGACATGATATGCCTGCTGGCTAGCCTCCCTGCCTCCATGAAAATAGACGTTCTTCTTCGTAGGGTCACGTTCCAGGCCGTGGTCGGCAAAATACATGACCGAGGCGCGGCGATCTTTTAATAATTCAAAAACCTGACCTAGCAAACTATCTGTGTAATGAATGGAGTTGTCATAGCAGGCATCCTGATCGTCCTG  >  NC_000913/3320201‑3320458

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