Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,494,973	G→A	*92* (TAG→TAA)	ypdI →	putative lipoprotein involved in colanic acid biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,494,973	0	G	A	100.0%	207.3 / NA	64	*92* (TAG→TAA)	ypdI	putative lipoprotein involved in colanic acid biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/33);  total (31/33)

CACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGA  >  NC_000913/2494837‑2495107                                                                                                                                         |                                                                                                                                       cacacTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTCATTCATCCTGTGCGATAACg                                                                                                                                      >  2:659194/1‑139 (MQ=255)        tACAAAACTTCAGTCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTca                                                                                                                               >  2:334207/1‑139 (MQ=255)                ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       <  2:640722/139‑1 (MQ=255)                ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       >  2:1334004/1‑139 (MQ=255)                       aTCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCCCAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACa                                                                                                                <  2:330271/139‑1 (MQ=255)                          acaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATc                                                                                                             <  2:677479/139‑1 (MQ=255)                          acaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATc                                                                                                             <  1:781740/139‑1 (MQ=255)                                 aTTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGaa                                                                                                                  >  1:103030/1‑127 (MQ=255)                                 aTTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGaa                                                                                                                  <  2:103030/127‑1 (MQ=255)                                      gTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 <  2:1298007/139‑1 (MQ=255)                                          aaGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                                                            <  1:747107/108‑1 (MQ=255)                                          aaGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                                                            >  2:747107/1‑108 (MQ=255)                                            gATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       >  1:1220279/1‑111 (MQ=255)                                            gATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       <  2:1220279/111‑1 (MQ=255)                                                 gagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAacac                                                                                      <  1:1334004/139‑1 (MQ=255)                                                       aTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 >  1:1158276/1‑122 (MQ=255)                                                       aTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 <  2:1158276/122‑1 (MQ=255)                                                       aTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGccc                                                                                >  2:767714/1‑139 (MQ=255)                                                         gcgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCaa                                                                              <  2:271621/139‑1 (MQ=255)                                                         gcgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCaa                                                                              <  1:334207/139‑1 (MQ=255)                                                                    aaaCCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATcc                                                                   <  1:1033183/139‑1 (MQ=255)                                                                        cTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAc                                                               <  1:767714/139‑1 (MQ=255)                                                                          cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       >  2:1320016/1‑81 (MQ=255)                                                                          cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       <  1:1320016/81‑1 (MQ=255)                                                                          cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                         >  2:836799/1‑127 (MQ=255)                                                                          cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAccc                                                             >  2:435558/1‑139 (MQ=255)                                                                          cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGAccc                                                             >  1:1150949/1‑139 (MQ=255)                                                                          cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGAGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                         <  1:836799/127‑1 (MQ=255)                                                                             agaagaAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCatt                                                          >  2:1179185/1‑137 (MQ=255)                                                                                  aaaTGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTTATCCTGTGCGATAACGTc                                                                                                                                    <  2:748253/60‑1 (MQ=255)                                                                                  aaaTGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTc                                                                                                                                    >  1:748253/1‑60 (MQ=255)                                                                                       ccgccCCCTTTTCCGACAAAGGCGAGTCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGcc                                                >  1:1111175/1‑139 (MQ=255)                                                                                       ccgccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGcc                                                >  1:1284143/1‑139 (MQ=255)                                                                                         gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                         >  1:428824/1‑112 (MQ=255)                                                                                         gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                         <  2:428824/112‑1 (MQ=255)                                                                                           ccccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAacac                                                                                      >  2:793904/1‑97 (MQ=255)                                                                                           ccccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAacac                                                                                      <  1:793904/97‑1 (MQ=255)                                                                                               ttttCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGaa                                                                                                                  <  1:902678/65‑1 (MQ=255)                                                                                               ttttCCGACAAAGCCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGaa                                                                                                                  >  2:902678/1‑65 (MQ=255)                                                                                                         aaGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGc                                                                                                                         <  2:102025/48‑1 (MQ=38)                                                                                                         aaGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGc                                                                                                                         >  1:102025/1‑48 (MQ=38)                                                                                                           ggCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                                <  2:884397/71‑1 (MQ=255)                                                                                                           ggCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                                >  1:884397/1‑71 (MQ=255)                                                                                                             cGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGt                                                                                                                     <  1:666015/48‑1 (MQ=38)                                                                                                             cGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGt                                                                                                                     >  2:666015/1‑48 (MQ=38)                                                                                                                      ttATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       <  2:1062712/37‑1 (MQ=38)                                                                                                                      ttATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                       >  1:1062712/1‑37 (MQ=38)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 <  2:1110555/51‑1 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 >  1:1110555/1‑51 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATcc                                                                   >  1:794218/1‑81 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATcc                                                                   <  2:794218/81‑1 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                               <  1:72159/101‑1 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                               >  2:72159/1‑101 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAccc                                          <  1:676336/106‑1 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAccc                                          >  2:676336/1‑106 (MQ=255)                                                                                                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg         >  1:146352/1‑139 (MQ=255)                                                                                                                                 gtgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                         <  1:994166/72‑1 (MQ=255)                                                                                                                                 gtgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                         >  2:994166/1‑72 (MQ=255)                                                                                                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 <  1:76161/47‑1 (MQ=255)                                                                                                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                                 >  2:76161/1‑47 (MQ=255)                                                                                                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa     >  2:350120/1‑139 (MQ=255)                                                                                                                                   gCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCaa    <  1:1179185/139‑1 (MQ=255)                                                                                                                                    cGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg   <  2:1284143/139‑1 (MQ=255)                                                                                                                                     gATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGa  <  1:350120/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CACACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGA  >  NC_000913/2494837‑2495107

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