Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,465,722	Δ1 bp	coding (422/933 nt)	yfdC →	putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,465,719	0	G	.	100.0%	273.8 / NA	57	coding (419/933 nt)	yfdC	putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (36/21);  total (36/21)

GCTGGAGAATCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTATGGGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGTGATGAAGAACACCCCCAGCGAGATGTT  >  NC_000913/2465576‑2465863                                                                                                                                                |                                                                                                                                                 gCTGGAGAATCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGcgtcg                                                                                                                                          <  1:139510/151‑1 (MQ=255)    ggAGAATCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGtg                                                                                                                                        <  1:1016853/150‑1 (MQ=255)    ggAGAATCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGtg                                                                                                                                        <  1:741637/150‑1 (MQ=255)      agaATCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGtgctg                                                                                                                                     >  2:361885/1‑151 (MQ=255)         aTCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGctgctg                                                                                                                                  >  2:645682/1‑151 (MQ=255)          tctcGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTgg                                                                                                                                 >  2:1062753/1‑151 (MQ=255)            tcGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGt                                                                                                                               >  2:820244/1‑151 (MQ=255)              ggTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTa                                                                                                                              >  2:592249/1‑150 (MQ=255)               gTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTa                                                                                                                              <  1:802160/149‑1 (MQ=255)               gTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAat                                                                                                                            >  2:57107/1‑151 (MQ=255)                 tataCCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAatat                                                                                                                          >  1:411277/1‑151 (MQ=255)                   taCCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATAtt                                                                                                                         >  2:249657/1‑150 (MQ=255)                    aCCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATtct                                                                                                                       <  1:1027856/151‑1 (MQ=255)                     ccTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATtctc                                                                                                                      <  1:433270/151‑1 (MQ=255)                        ttGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCg                                                                                                                     <  2:950654/149‑1 (MQ=255)                           gTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGAc                                                                                                                 >  2:548031/1‑150 (MQ=255)                            ttttATTACCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGAc                                                                                                                 >  1:931483/1‑149 (MQ=255)                               tattatCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGt                                                                                                             >  1:52137/1‑150 (MQ=255)                                 ttatCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTAtt                                                                                                          >  1:562448/1‑151 (MQ=255)                                      gTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCACCGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTgcggc                                                                                                     <  2:221683/151‑1 (MQ=255)                                      gTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTgcgg                                                                                                      >  1:96805/1‑150 (MQ=255)                                       tCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGCTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTgcggcg                                                                                                    >  2:1011976/1‑151 (MQ=255)                                            aTGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGc                                                                                               <  2:874632/151‑1 (MQ=255)                                                 ccGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCAttt                                                                                           >  2:226696/1‑150 (MQ=255)                                                 ccGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTg                                                                                          >  2:1048774/1‑151 (MQ=255)                                                 ccGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTg                                                                                          <  1:669272/151‑1 (MQ=255)                                                  cGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGa                                                                                         <  1:937414/151‑1 (MQ=255)                                                    ccAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAat                                                                                       >  2:952001/1‑151 (MQ=255)                                                            tattTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGcc                                                                                >  1:549381/1‑150 (MQ=255)                                                               ttACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTAt                                                                             <  1:450929/150‑1 (MQ=255)                                                                 aCCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCtt                                                                          <  1:171731/151‑1 (MQ=255)                                                                   cGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTc                                                                         >  1:785341/1‑150 (MQ=255)                                                                      aaaTACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATg                                                                     >  1:366279/1‑151 (MQ=255)                                                                       aaTACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATg                                                                     >  1:19803/1‑150 (MQ=255)                                                                           ccGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATgaagaa                                                                <  1:521310/151‑1 (MQ=255)                                                                              tGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACt                                                             >  2:614395/1‑151 (MQ=255)                                                                               gACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTc                                                            >  1:739137/1‑151 (MQ=255)                                                                                 cTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTcgc                                                          >  2:920761/1‑151 (MQ=255)                                                                                  tGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTcgcg                                                         >  1:370975/1‑151 (MQ=255)                                                                                     ggTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATg                                                      >  1:143078/1‑151 (MQ=255)                                                                                       tactacCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGc                                                     >  2:570583/1‑150 (MQ=255)                                                                                        actacCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCa                                                    <  1:565192/150‑1 (MQ=255)                                                                                           acCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTg                                                <  2:651070/151‑1 (MQ=255)                                                                                            cccGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGt                                               <  2:913218/151‑1 (MQ=255)                                                                                                   tGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATc                                        >  2:335111/1‑151 (MQ=255)                                                                                                       aaaaCCGACAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGc                                     >  1:535058/1‑150 (MQ=255)                                                                                                               cAATGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGAtgtg                            >  2:901364/1‑151 (MQ=255)                                                                                                                aaTGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGAtgtg                            >  1:1000702/1‑150 (MQ=255)                                                                                                                 aTGAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGtgat                          >  1:209574/1‑151 (MQ=255)                                                                                                                   gAGCAACGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGtgatg                         <  1:592147/150‑1 (MQ=255)                                                                                                                       aaCGTCGGCTTACTTATACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGTGATgaagaa                    <  2:916524/151‑1 (MQ=255)                                                                                                                               cttacttaTACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAGTCGGCATGGATGTGATGAAGAACAccccc             >  2:1060447/1‑150 (MQ=255)                                                                                                                                 tacttaTACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCGAAATCGGCATGGATGTGATGAAGAACACCCCCAg           >  1:187268/1‑150 (MQ=255)                                                                                                                                   cttaTACGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGTGATGAAGAACACCCCCAGCg         >  1:158502/1‑150 (MQ=255)                                                                                                                                      ataCGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTAAAAATCGGCATGGATGTGATGAAGAACACCCCCAGCGAGAt     <  2:438236/151‑1 (MQ=255)                                                                                                                                         cGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGTGATGAAGAACACCCCCAGCGAGATGtt  >  1:692860/1‑151 (MQ=255)                                                                                                                                         cGGTTAT‑GGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGTGATGAAGAACACCCCCAGCGAGATGtt  <  1:206416/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 GCTGGAGAATCTCGGTTATACCTTTGGTTTTATTATCGTCATTATGGCCCGCCAGCAATTATTTACCGAAAATACCGTGACTGCGGTACTACCCGTCATGCAAAAACCGACAATGAGCAACGTCGGCTTACTTATACGGTTATGGGGCGTCGTGCTGCTGGGTAATATTCTCGGGACAGGTATTGCGGCGTGGGCATTTGAATATATGCCTATCTTCAATGAAGAAACTCGCGATGCATTTGTCAAAATCGGCATGGATGTGATGAAGAACACCCCCAGCGAGATGTT  >  NC_000913/2465576‑2465863

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